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Abstract 

Seafloor  digital  imagery  and  bathymetric  data  are  used  to  evaluate  the  volcanic 
characteristics  of  the  85°E  segment  of  the  ultraslow  spreading  Gakkel  Ridge  (9  mm  yr'1). 
Imagery  reveals  that  ridges  and  volcanic  cones  in  the  axial  valley  are  covered  by 
numerous,  small- volume  lava  flows,  including  a  few  flows  fresh  enough  to  have 
potentially  erupted  during  the  1999  seismic  swarm  at  the  site.  The  morphology  and 
distribution  of  volcaniclastic  deposits  observed  on  the  seafloor  at  depths  of  ~3800  m, 
greater  than  the  critical  point  for  steam  generation,  are  consistent  with  having  formed  by 
explosive  discharge  of  magma  and  C02  from  source  vents. 

Microearthquakes  recorded  on  a  200  m  aperture  seismometer  network  deployed  on  the 
Trans-Atlantic  Geotraverse  active  mound,  a  seafloor  massive  sulfide  on  the  Mid- Atlantic 
Ridge  at  26°N,  are  used  to  image  subsurface  processes  at  the  hydrothermal  system.  Over 
nine-months,  32,078  local  microearthquakes  (ML  =  -1)  with  single-phase  arrivals  cluster 
on  the  southwest  flank  of  the  deposit  at  depths  <125  m.  Microearthquakes  characteristics 
are  consistent  with  reaction-driven  cracking  driven  by  anhydrite  deposition  in  the  shallow 
secondary  circulation  system.  Exit  fluid  temperatures  recorded  at  diffuse  vents  on  the 
mound  during  the  microearthquake  study  are  used  to  explore  linkages  between  seismicity 
and  venting. 
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Chapter  1:  Introduction 

The  mid-ocean  ridges,  the  Earth’s  longest  system  of  mountains,  mark  the  limbs  of  rising 
convection  cells  in  the  mantle,  where  rigid  tectonic  plates  move  apart  and  new  seafloor 
material  rises  to  fill  the  vacant  space.  The  generation  of  new  plate  material  (i.e.  oceanic 
lithosphere)  occurs  by  seafloor  spreading  along  the  active  mid-ocean  ridges.  Oceanic 
crust  is  formed  by  the  intrusion  of  basaltic  magma  from  below,  and  contact  with  seawater 
convection  of  seawater  through  the  crust  provides  an  efficient  cooling  mechanism. 
Convection  of  the  seawater  through  the  crust  can  result  in  hydrothermal  venting,  where 
hot  hydrothermal  fluids  are  formed  through  interactions  of  seawater  with  oceanic  crust, 
modifying  the  composition  of  the  crust,  affecting  ocean  chemistry,  building  metal-rich 
deposits,  and  providing  energy  for  deep  sea  biological  communities.  Both  volcanism  and 
hydrothermal  circulation  transfer  mass  and  energy  from  the  crust  and  mantle  to  the 
oceans.  Volcanic  processes  at  mid-ocean  ridges  (MORs)  generate  the  oceanic  crust 
comprising  more  than  60%  of  the  Earth’s  solid  surface.  Hydrothermal  circulation  through 
permeable  crust,  driven  by  heat  [Lister,  1972;  Strens  and  Cann,  1986],  contributes  to 
transport  of  heat  from  the  lithosphere  to  the  oceans,  which  is  estimated  to  account  for 
30%  of  the  predicted  global  oceanic  heat  flux  and  almost  25%  of  the  total  heat  flux  of  the 
Earth  [ Sclater  el  al.,  1980;  Stein  and  Stein,  1994]. 

Spreading  rates  along  the  global  MOR  system  range  from  ultrafast  (>160  mm  kyr  1  full 
rate)  to  ultraslow  (<12  mm  kyr  ').  Plate  spreading  rate  is  perhaps  the  dominant  parameter 
controlling  crustal  accretion  and  volcanic  processes  at  MORs,  with  faster  spreading 
ridges  exhibiting  a  higher  magma  supply  rate  and  a  hotter  thermal  structure  compared  to 
slower  spreading  ridges  [Macdonald  et  al.,  1982;  Smith  and  Cann,  1993;  Perfit  and 
Chadwick,  1998;  Edwards  et  al.,  2001;  Sinton  et  al.,  2002;  Cann  and  Smith,  2005; 
Cochran,  2008;  Searle  et  al.,  2010].  Fast-spreading  ridges  experience  frequent  (order  10- 
100  yrs),  small- volume  volcanic  eruptions  on  narrow  neovolcanic  zones  along  axis, 
believed  to  reflect  a  high,  continuous  magma  supply  rate  to  the  ridge  [e.g.,  Macdonald  et 
al.,  1982;  Perfit  and  Chadwick,  1998;  Sinton  et  al.,  2002].  The  resulting  ridge 
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morphology  is  an  axial  high  with  a  narrow  graben  along  the  ridge  crest  [Small,  1994; 
Fornari  et  al.,  2004;  Soule  et  al.,  2009] .  At  slow- spreading  ridges,  the  dominant 
morphologic  feature  is  a  rift  valley,  believed  to  reflect  lower  and  intermittent  magma 
supply  rates  and  increased  tectonism  [ Langmuir  et  al.,  1992;  Rubin  and  Sinton,  2007; 
Behn  and  Ito,  2008].  Volcanic  eruptions  at  slow  spreading  ridges  are  thought  to  be 
infrequent  (order  1000-10,000  yrs)  and  produce  discrete  point- source  volcanic 
constructions  across  the  rift  valley  floor  [ Smith  and  Cann,  1993;  Edwards  et  al.,  2001; 
Cann  and  Smith,  2005;  Cochran  et  al,  2008;  Searle  et  al.,  2010]. 

Melting  models  of  MORs  predict  a  dramatic  decrease  in  magma  production  and  crustal 
thickness  with  decrease  in  spreading  rate  [ Reid  and  Jackson,  1981;  Bown  and  White, 
1994].  Based  on  sparse  observations  it  has  also  been  suggested  that  eruption  volume 
should  vary  inversely  with  spreading  rate  with  the  largest  eruption  volumes  at  ultraslow 
spreading  ridges  [Perfit  and  Chadwick,  1998;  Sinton  et  al.,  2002] .  Sinton  et  al.  [2002] 
hypothesize  that  because  ephemeral  magma  chambers  at  ultraslow  spreading  ridges  lie 
beneath  thicker,  stronger  lithospheric  lids,  they  require  larger  pressures  for  eruption, 
leading  to  longer  repose  intervals  and  larger  volume  eruptions  (~10-100  km3).  The 
Gakkel  Ridge  figures  prominently  in  this  model  as  an  end-member  in  terms  of  spreading 
rate  and  eruption  volume,  and  the  hydrothermally  and  volcanically  active  85°E  segment 
is  particularly  well-suited  to  serve  as  a  study  site  for  observing  recent  volcanism  at 
ultraslow  spreading  rates. 

Melting  models  of  MORs  also  predict  that  a  decrease  in  crustal  magma  storage  with 
decrease  in  spreading  rate  should  result  in  a  corresponding  decrease  in  hydrothermal 
activity  [Baker,  et  al.,  1994].  Although  a  compilation  of  observations  at  ridge  segments 
of  different  spreading  rates  shows  a  robust  linear  trend  between  spreading  rate  and  spatial 
frequency  of  hydrothermal  venting  [Baker  et  al.,  2004],  investigations  of  vent  fields 
along  fast-,  medium  and  slow-spreading  ridges  have  shown  that  the  largest  hydrothermal 
deposits  tend  to  be  found  on  slower-spreading  ridges  [Hannington  et  al.,  1995;  Tivey, 
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2007].  Seafloor  massive  sulfide  deposits  found  at  slow-spreading  ridges,  like  the  TAG 
active  mound  on  the  MAR,  are  the  modern  equivalents  of  ancient  ophiolite-hosted  ore 
deposits,  such  as  those  mined  for  copper  and  gold  in  Cyprus,  Oman,  and  Newfoundland 
[Humphris  et  al.,  1995;  Hannington  etal.,  1998;  Humphris  and  Cann,  2000].  Many 
questions  remain  in  our  understanding  of  the  factors  that  control  focused  hydrothermal 
venting  and  deposit-building  on  the  seafloor  [c.f.  German  and  Von  Damm,  2004; 
Hannington  et  al.,  2005],  but  at  many  Atlantic  vent  fields,  including  the  TAG 
hydrothermal  field,  large  hydrothermal  deposits  are  associated  with  detachment  faults 
[deMartin  et  al.,  2007],  suggesting  that  detachment  faults  commonly  act  as  channels  for 
focusing  of  primary  hydrothermal  circulation  at  slow- spreading  ridges  [McCaig  et  al., 
2007] .  Genesis  of  a  seafloor  massive  sulfide  deposit  requires  a  magma  supply,  which  is 
thought  to  control  detachment  fault  formation  [Behn  and  Ito,  2008] ,  a  primary  system  of 
hydrothermal  circulation  that  extracts  heat  from  the  source,  and  a  secondary  circulation 
system  that  entrains  seawater  to  cool  the  shallow  deposit  and  fuel  mineral  precipitation 
[Sleep,  1991;  Tivey  etal.,  1995;  Lowell  and  Yao,  2002;  Lowell  et  al.,  2003]. 

The  purpose  of  this  thesis  is  to  better  understand  volcanic  and  hydrothermal  processes  in 
two  particular  mid-ocean  ridge  settings:  the  ultraslow-spreading  Gakkel  Ridge  at  85°E  in 
the  Arctic  Ocean  and  the  TAG  active  mound  in  the  TAG  hydrothermal  field  on  the  slow- 
spreading  Mid- Atlantic  Ridge  at  26°08’N.  The  work  presented  in  this  thesis  was 
designed  to: 

a)  identify  and  determine  volcanic  characteristics  of  the  ultraslow-spreading  Gakkel 
Ridge  at  85°E  on  the  scale  of  individual  eruptive  products  (Chapter  2), 

b)  describe  the  microearthquake  activity  local  to  the  TAG  active  mound  seafloor 
massive  sulfide  deposit  (Chapter  3), 

c)  better  understand  from  microearthquake  observations,  the  processes  of  secondary 
fluid  circulation  and  mineral  precipitation  within  the  hydrothermal  system 
(Chapter  3), 
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d)  assess  the  relationship  between  local  and  regional  seismicity  and  exit  fluid 
venting  temperatures  at  the  TAG  active  mound  (Chapter  4). 


1.1  Gakkel  Ridge,  85 °E  segment 

The  Gakkel  Ridge  is  among  the  slowest-spreading  sections  of  the  global  MOR  system. 
The  Gakkel  Ridge  and  the  Southwest  Indian  Ridge  fall  within  a  unique  classification  of 
ridges  termed  “ultraslow  spreading,”  which  are  characterized  by  alternating  magmatic 
and  amagmatic  segments,  and  a  lack  of  transform  faulting  [Cannat  et  al.,  2003;  Dick  et 
al.,  2003].  Limited  seafloor  mapping  studies  have  been  conducted  at  ultraslow  spreading 
ridges  [e.g.,  Dick  et  al.,  2003;  Standish  et  al.,  2008;  Michael  et  al.,  2003;  Edwards  et  al., 
2003]  and  even  fewer  with  observations  at  the  scale  of  individual  eruptions  owing  to  the 
logistical  and  technical  challenges  associated  with  working  in  the  rough  seas  of  the 
Southern  Ocean  (i.e.,  Southwest  Indian  Ridge)  and  under  the  perennial  ice  pack  of  the 
Arctic  Ocean  (i.e.,  Gakkel  Ridge),  where  a  large  fraction  of  ultraslow  spreading  ridges 
occur.  The  1800  km-long  Gakkel  Ridge  extends  from  near  Greenland  (12.8  mm  yr  1  full 
rate)  across  the  Eastern  Arctic  Basin  to  near  the  Siberian  margin  (6.5  mm  yr  1  full  rate) 
[Vogt  et  al.,  1979;  DeMets  et  al.,  1994;  Sella  et  al.,  2002] .  For  comparison,  the  ultraslow- 
spreading  Southwest  Indian  Ridge  spreads  at  rates  from  12  to  16  mm  yr1,  although 
obliquely  spreading  segments  may  experience  effective  spreading  rates  below  8  mm  yr 1 
[DeMets  et  al.,  1994;  Dick  et  al.,  2003]. 

In  1999  the  Global  Seismic  Network  detected  the  largest  MOR  earthquake  swarm  in 
recorded  history  at  85°E  [Muller  and  Jokat,  2000],  Earthquake  magnitudes  ranged  from 

3.2  to  5.8  and  some  events  had  non-double  couple  focal  mechanisms  similar  to  those  of 
known  volcanic  events  [Korger  and  Schlindwein,  2011].  The  swarm  lasted  >  7  months. 
The  SCICEX  expedition  of  the  Gakkel  Ridge  followed  the  detected  seismicity  with  a 
sidescan  sonar  survey  and  imaged  a  large  region  of  high  backscatter  intensity  interpreted 
as  unsedimented  lava  flows  [Edwards  et  al.,  2001].  The  2001  AMORE  expedition 
returned  to  the  study  area  to  conduct  dredging  and  water  column  surveys.  The  most 
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vigorous  hydrothermal  venting  detected  along  the  ridge  was  found  in  the  same  location 
(85°E)  where  a  1400  m  thick  megaplume  was  observed  [ Edmonds  et  al.,  2003]. 
Seismoacousitic  arrays  recorded  numerous  microearthquakes  (average  magnitude  = 

-0.33)  during  the  2001  AMORE  study  that  are  proposed  to  have  been  magmatic  in  origin 
[Schlindwein  et  al.,  2007;  Schlindwein  and  Riedel,  2010].  From  this  evidence,  it  was 
hypothesized  that  the  85°E  segment  is  a  hydrothermally  and  volcanically  active  area  that 
experienced  eruption  in  1999. 

In  Chapter  2,  we  use  samples  and  images  acquired  from  the  site  with  a  purpose-built 
camera  and  sampling  system  (CAMPER)  during  the  AGAVE  expedition  in  2007  to 
constrain  the  modes  of  volcanic  accretion  at  this  remote  site,  including  what  appears  to  be 
ubiquitous  explosive  activity  during  the  1999  earthquake  swarm.  Imagery  revealed  that 
major  portions  of  the  axial  valley  are  blanketed  with  volcaniclastic  deposits  interpreted  to 
result  from  magma  fragmentation  during  explosive  eruptions  [Sohn  et  al.,  2008],  We 
expand  on  evidence  for  magmatically-driven  explosive  eruptions  by  using  the  complete 
suite  of  dive  images  to  constrain  volcanic  processes  at  the  scale  of  individual  lava  flows. 
Pyroclastic  eruption  in  the  deep  ocean  is  a  topic  of  controversy  because  mid-ocean  ridge 
eruptions  are  generally  dominated  by  effusive  eruptions  of  basalt  with  low  volatile 
contents.  External  water  is  thought  to  be  unable  to  drive  explosive  eruptions,  especially  at 
these  depths  due  to  the  lack  of  boiling.  However  explosive  volcanic  deposits  have  been 
documented  at  similarly  deep  oceanic  spreading  centers  [Clague  et  al.,  2003,  2009]  and 
mechanisms  for  submarine  explosive  eruptions  have  been  evaluated  in  a  variety  of  arc, 
near-arc,  and  ocean  island  settings  [Wright,  1994;  Clague  et  al.,  2000,  2003,  2009; 
Chadwick  et  al.,  2008;  Schipper  et  al.,  2010;  Deardorff  et  al.,  2011;  Helo,  et  al.,  2011; 
Resing  et  al.,  2011].  Our  dive  imagery  demonstrates  that  the  seafloor  in  this  area  is 
covered  by  small  lava  flows  and  volcaniclastic  deposits  within  the  large  high-backscatter 
region  [Edwards  et  al.,  2003].  We  use  our  dive  images  to  develop  a  new  model  for 
volcanic  processes  at  the  85°E  segment,  and  discuss  the  implications  of  this  model  for 
ultraslow  spreading  volcanic  processes,  in  general. 
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In  the  future,  characterizing  the  individual  lava  flows  and  potential  products  of 
pyroclastic  eruption  will  help  to  answer  remaining  questions  about  volcanic  processes  on 
deep  and  ultraslow  spreading  ridges,  such  as: 

a)  How  are  volcaniclastic  materials  formed  and  dispersed  during  deep-sea  eruptions? 

b)  How  are  magmatic  volatiles  sequestered  and/or  how  is  magma  fragmented? 

c)  How  are  effusive  and  explosive  eruptions  related? 

d)  How  robust  is  the  relationship  between  eruption  volume  and  spreading  rate? 

These  questions  could  potentially  be  answered  by  more  direct  observations  of  the 
seafloor  on  deep  ridges.  Identifying  the  boundaries/extent  of  an  individual  volcaniclastic 
eruption  product,  in  particular,  may  add  needed  information  on  the  formation  and 
dispersal  of  clasts. 


1.2  The  TAG  active  mound,  Mid-Atlantic  Ridge  (26°08’N) 

The  Trans-Atlantic  Geotraverse  (TAG)  segment  on  the  slow  spreading  Mid- Atlantic 
Ridge  at  26°08’N  (half-spreading  rates  of  13  mm  year 1  to  the  east  and  11  mm  year  1  to 
the  west)  [McGregor  and  Rona,  1975;  McGregor  et  al.,  1977]  hosts  one  of  the  largest 
known  and  best  studied  hydrothermal  systems  on  the  seafloor  (e.g.,  Scott  et  al.,  1974a, 
1974b;  Rona  et  al.,  1975,  1986;  Campbell  et  al.,  1988;  Karson  and  Rona,  1990;  Herzig  et 
al.,  1998,  etc.)  The  TAG  active  mound,  a  seafloor  massive  sulfide  located  on  the 
extensional  hanging  wall  of  a  detachment  fault  at  3670  m  depth  [Tivey  et  al.,  2003; 
deMartin  et  al.,  2007],  has  been  the  subject  of  numerous  multidisciplinary  studies, 
including  segment-  and  mound-scale  multibeam  bathymetry  and  side-scan  sonar 
[Humphris  and  Kleinrock,  1996;  Kleinrock  and  Humphris,  1996;  White  et  al.,  1998; 
Humphris  and  Cann,  2000];  geophysical  surveys  [Evans,  1996;  Tivey  et  al.,  2003; 
Canales  et  al.,  2007;  deMartin  et  al.,  2007];  and  investigations  of  rocks  and  vent  fluids 
from  the  TAG  active  mound  [Humphris  et  al.,  1995;  Tivey  et  al.,  1995;  Sohn  et  al., 

2007a, b]. 
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Microearthquake,  magnetic,  and  seismic  refraction  studies  on  the  ridge  detected  a  zone  of 
demagnetization  and  low-velocity  attributed  to  the  extensional  environment  of  a  ~5  km 
long  shallow  dipping  (20°)  dome-shaped  detachment  fault  [ Tivey  et  al.,  2003;  Canales  et 
al.,  2007;  deMartin  et  al.,  2007]  hosting  the  TAG  hydrothermal  field,  a  zone  of  low- 
temperature  alteration  and  three  zones  of  high-temperature  alteration  with  an  assemblage 
of  active  and  relict  hydrothermal  mounds  up  to  50,000  years  old  [ Thompson  et  al.,  1985; 
Lalou  et  al.,  1990;  Rona  et  al.,  1993a,b].  The  seismogenic  zone  of  the  detachment  fault 
extends  to  7-8  km  below  the  seafloor  and  may  serve  as  a  pathway  for  deep  hydrothermal 
fluid  circulation  that  fuels  high  temperature  venting  at  the  active  TAG  mound  [ deMartin 
et  al.,  2007].  High  temperature  activity  is  currently  focused  at  the  TAG  active  mound, 
one  of  at  least  four  large  hydrothermal  deposits  in  the  field  [  White  et  al.,  1998], 

The  TAG  active  mound  is  a  200  m  diameter,  50  m  tall  massive  sulfide  deposit  that 
discharges  focused  hydrothennal  fluids  in  excess  of  360°C  from  a  Black  Smoker 
Complex  (B.S.C.)  of  sulfide  chimneys  at  the  top  of  the  mound  and  diffuse,  lower 
temperature  fluids  from  many  orifices  on  the  mound  surface  [ Campbell  et  al.,  1988; 
Humphris  et  al.,  1995;  Tivey  et  al.,  1995].  Within  and  beneath  the  active  mound,  buoyant 
high  temperature  fluids  are  focused  in  an  -80  m  wide  upflow  zone  to  the  Black  Smoker 
Complex  [Honnorez  et  al.,  1998].  Relatively  cooler  fluids  derived  from  entrained 
seawater  that  is  conductively  heated  and  mixed  with  black  smoker  fluids  are  diffusely 
vented  across  the  mound.  Time  scales  of  flow  reversals  in  exit  fluid  temperature  data  and 
measurements  of  low  conductive  heat  flow  are  indicative  of  seawater  entrainment  into  the 
mound  surface,  and  occur  on  complex  episodic  and  periodic  time  scales  [ Becker  et  al., 
1996;  Sohn,  2007a, b]. 

In  spite  of  the  large  body  of  knowledge  on  the  TAG  active  mound,  some  key  aspects  of 
the  hydrothermal  system  remain  unclear.  In  particular,  relatively  little  is  known  about  the 
physical  extent  and  geometry  of  hydraulic  pathways  and  subsurface  fluid  flow.  Defining 
the  zone  of  secondary  recharge  is  important  to  understanding  where  mineral  deposition, 
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especially  that  of  anhydrite,  occurs  within  the  mound.  There  is  strong  evidence  from  the 
TAG  active  mound  subsurface  that  faulting  is  necessary  to  maintain  shallow  hydraulic 
pathways,  and  studies  on  fast  spreading  mid-ocean  ridge  systems  have  shown  a  direct 
correlation  between  fluid  flow  and  local  seismicity  [Sohn  et  al.,  1998;  Johnson  et  al., 
2000;  Davis  et  al.,  2004] . 

In  Chapter  3,  we  investigate  hydraulic  phenomenon  by  deploying  a  dense  small  aperture 
array  of  ocean  bottom  seismometers  in  a  ring  on  the  periphery  of  the  TAG  active  mound 
seafloor  massive  sulfide.  We  detected  32,078  microearthquakes  in  a  9-month 
deployment.  We  model  the  microearthquakes  as  reaction-driven  fracturing  events  in 
response  to  the  deposition  of  anhydrite  in  the  secondary  circulation  system,  and  use  this 
model  to  constrain  the  rate  and  subsurface  location  of  anhydrite  deposition  at  the  active 
mound. 

In  Chapter  4,  we  use  statistical  methods  to  assess  the  impact  of  local  and  regional 
microearthquake  activity  on  exit-fluid  temperatures  at  the  TAG  active  mound.  Regional 
microearthquakes  occur  on  the  detachment  fault  hosting  the  TAG  hydrothermal  field,  as 
well  as  antithetic  normal  faults.  Local  microearthquakes  described  in  Chapter  3,  are 
located  in  the  shallow  subsurface  (<125  mbsf)  directly  on  the  periphery  of  the  TAG 
active  mound,  and  are  thought  to  be  intrinsic  to  hydrothermal  processes.  Comparison  of 
microearthquake  catalogs  with  concurrent  exit  fluid  temperature  venting  records  allows 
us  to  look  for  feedbacks  that  suggest  that  microearthquake  activity  affects  subsurface 
permeability  and  reflects  the  maintenance  of  hydraulic  pathy ways  in  the  hydrothermal 
system.  Results  of  initial  tests  suggest  that  more  refined  methods  may  be  successful  in 
defining  spatial  and  temporal  relationships  between  subsurface  cracking  and  hydraulic 
pathways  to  surface  vents. 
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In  the  future,  characterizing  the  precise  spatial  distribution  of  subsurface  mineralization 
and  the  secondary  recharge  zone  on  hydrothermal  mounds  will  help  to  answer  several 
remaining  questions  about  their  behavior,  such  as: 

a)  How  does  the  subsurface  permeability  of  a  seafloor  massive  sulfide  change  in 
time? 

b)  How  broad  is  the  source  and  depth  of  seawater  entrainment,  that  triggers  a  great 
deal  of  mineral  precipitation,  especially  that  of  anhydrite,  in  the  mound? 

These  questions  could  potentially  be  answered  by  using  a  more  dense,  small  aperture 
deployment  of  seismometers  on  seafloor  hydrothermal  systems.  Currently  few 
submarine  studies  of  hydrothermal  systems  exist.  Understanding  permeability  and  zones 
of  mineral  deposition  and  recharge  are  also  important  to  those  interested  in  the  subsurface 
biosphere  and  thermal  structure  of  the  mound. 

Examined  as  a  whole,  this  thesis  attempts  to  contribute  to  the  growing  body  of  work  on 
slow-  and  ultraslow-spreading  mid-ocean  ridges  and  the  vast  wealth  of  knowledge  and 
expertise  on  the  TAG  active  mound  and  lesser-studied  seafloor  massive  deposits.  This 
progress  towards  characterizing  the  volcanic  structure  of  the  remote  Gakkel  Ridge,  the 
ultraslow-spreading  end  member  of  our  global  mid-ocean  ridge  system,  sheds  light  on 
many  questions  about  deep-sea  explosive  volcanism.  Studies  at  the  TAG  active  mound 
make  progress  towards  the  geophysical  imaging  of  sulfide  deposits  and  a  better 
understanding  of  the  short-term  hydrothermal  processes  necessary  to  build  and  maintain 
seafloor  massive  sulfide  deposits. 
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Chapter  2:  Effusive  and  explosive  volcanism  on  the  ultraslow-spreading 
Gakkel  Ridge,  85°E* 


Abstract 

We  use  high-definition  seafloor  digital  imagery  and  multi-beam  bathymetric  data 
acquired  during  the  2007  Arctic  Gakkel  Vents  Expedition  (AGAVE)  to  evaluate  the 
volcanic  characteristics  of  the  85°E  segment  of  the  ultraslow  spreading  Gakkel  Ridge  (9 
mm  yr 1  full  rate).  Our  seafloor  imagery  reveals  that  the  axial  valley  is  covered  by 
numerous,  small-volume  (order  ~1000  nr )  lava  flows  displaying  a  range  of  ages  and 
morphologies  as  well  as  unconsolidated  volcanic lastic  deposits  with  thicknesses  up  to  10 
cm.  The  valley  floor  contains  two  prominent  volcanic  lineaments  made  up  of  axis- 
parallel  ridges  and  small,  cratered  volcanic  cones.  The  lava  flows  appear  to  have  erupted 
from  a  number  of  distinct  source  vents  within  the  ~12-15  km- wide  axial  valley.  Only  a 
few  of  these  flows  are  fresh  enough  to  have  potentially  erupted  during  the  1999  seismic 
swarm  at  this  site,  and  these  are  associated  with  the  Oden  and  Loke  volcanic  cones.  We 
model  the  widespread  volcaniclastic  deposits  we  observed  on  the  seafloor  as  having  been 
generated  by  the  explosive  discharge  of  C02  that  accumulated  in  (possibly  deep)  crustal 
melt  reservoirs.  The  energy  released  during  explosive  discharge,  combined  with  the 
buoyant  rise  of  hot  fluid,  lofted  fragmented  clasts  of  rapidly  cooling  magma  into  the 
water  column,  and  they  subsequently  settled  onto  the  seafloor  as  fall  deposits  surrounding 
the  source  vent. 


*  Reproduced  by  permission  of  American  Geophysical  Union.  Published  as:  Pontbriand, 
C.  W.,  S.  A.  Soule,  R.  A.  Sohn,  S.  E.  Humphris,  C.  Kunz,  H.  Singh,  K.  Nakamura,  M. 
Jakobsson,  and  T.  Shank  (2012),  Effusive  and  explosive  volcanism  on  the  ultraslow- 
spreading  Gakkel  Ridge,  85°E,  Geochemistry,  Geophysics,  Geosystems,  13(10),  1-22, 
doi:  1 0. 1 029/20 12GC004 187.  Copyright  2012  American  Geophysical  Union. 
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2007  Arctic  Gakkel  Vents  Expedition  (AGAVE)  to  evaluate  the  volcanic  characteristics  of  the  85  E  seg¬ 
ment  of  the  ultraslow  spreading  Gakkel  Ridge  (9  mm  yr  1  full  rate).  Our  seafloor  imagery  reveals  that 
the  axial  valley  is  covered  by  numerous,  small-volume  (order  ~1000  m1)  lava  flows  displaying  a  range 
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1.  Introduction 

[2]  Volcanic  processes  at  mid-ocean  ridges  (MORs) 
generate  the  oceanic  crust  comprising  more  than 
60%  of  the  Earth’s  solid  surface.  Plate  spreading 
rale  is  perhaps  the  dominant  parameter  controlling 
crustal  accretion  and  volcanic  processes  at  MORs, 
with  faster  spreading  ridges  receiving  a  higher 
magma  supply  rale  and  exhibiting  a  holler  thermal 
structure  compared  to  slower  spreading  ridges 
[Macdonald.  1982;  Smith  and  Catin.  1992;  Perfit 
and  Chadwick.  1998;  Edwards  ct  al..  2001;  Sin  ton 
ct  al.,  2002;  Conn  and  Smith.  2005;  Cochran, 
2008;  Searle  et  al..  20101.  Although  the  first-order 
effects  of  spreading  rale  on  ridge  morphology  (e.g., 
axial  rise  versus  rift  valley)  and  thermal  structure  are 
well-established  [Simon  and  Doric h.  1992;  Smith 
and  Cann ,  1993;  Small,  1994;  Fornari  et  al., 
2004;  Cann  and  Smith.  2005:  Rubin  and  Sinton. 
2007;  Helm  and  Ito,  2008;  Soule  et  al.,  2009],  the 
implications  for  volcanic  eruption  processes  (e.g., 
volume,  spatial  extent,  frequency,  effusion  rate,  and 
lava  composition)  are  not  frilly  understood.  Many 
seafloor-mapping  studies  have  constrained  die 
nature  of  volcanic  processes  across  the  global 
MOR  system  [e.g.,  Ballard  and  van  Andel.  1977; 
Bonatti  and  Harrison.  1988;  Perfit  and  Chadwick, 
1998;  Simon  et  al.,  2002;  Baker  et  al.,  2004;  Soule 
et  al,  2005.  2007],  however,  volcanic  processes  at 
the  most  slowly  diverging  plate  boundaries  arc 
poorly  constrained.  Limited  seafloor-mapping 
studies  have  been  conducted  at  ultraslow  spreading 
ridges  [e.g..  Edwards  el  al..  2001:  Dick  et  al..  2003; 
Edwards  and  Coakley,  2003:  Michael  et  al .  2003; 
Standish  el  al.,  2008]  and  even  fewer  with  obser¬ 
vations  at  die  scale  of  individual  eruptions  owing  to 
the  logistical  and  technical  challenges  associated 
with  working  in  the  rough  seas  of  the  Southern 
Ocean  (i.e„  Southwest  Indian  Ridge)  and  under  the 
perennial  ice  pack  of  die  Arctic  Ocean  (i.c.,  Gakkel 
Ridge),  where  a  large  fraction  of  ultraslow-spreoding 
ridges  occur. 

[sj  The  Gakkel  Ridge  exhibits  full  spreading  rates 
that  systematically  decrease  from  12.8  mm  yr  1  at 
its  western  end  (near  Greenland)  to  6.5  mm  yr  1  at 
its  eastern  end  (the  Siberian  margin)  fl rogt  et  al.. 


1979:  DeAlets  et  al..  1994;  Sella  el  al..  20021. 
The  latter  arc  the  slowest  MOR  spreading  rates 
found  on  Earth.  Melfrng  models  of  MORs  predict  a 
sharp  decrease  in  magma  production  and  crustal 
thickness  when  full  spreading  rates  are  lower  tliau 
~20  mm  yr  1  [ Reid  and  Jackson,  1981;  Bown  and 
mile,  1994],  suggesting  that  magma  supply  should 
be  reduced  along  the  entire  1800  km  length  of  tire 
ridge.  Bathymetric  surveys  and  dredging  studies 
conducted  during  the  AMORE  expedition  [Dick 
et  al..  2003;  Michael  et  al..  2003],  along  with 
geophysical  surveys  [Coakley  and  Cochran.  1998; 
Jokat  et  al..  2003 ].  have  demonstrated  that  the 
igneous  crust  does  become  'vanishingly  thin'  in 
some  places,  presumably  because  relatively  low  and 
intermittent  magma  supply  rates  lead  to  avoleanie 
extension  and  exhumation  of  ultra-mafic  lithologies 
[Smith  and  Cann.  1993;  Sinton  and  Detrich,  1992; 
Cannat.  1996;  Dick  et  al.  2003;  Cann  and  Smith. 
2005;  Rubin  and  Sinton,  2007;  Bchn  and  Ito,  2008], 
However,  the  ridge  exhibits  significant  along-axis 
variations  in  magma  supply,  and  some  segments 
appear  to  be  magmatically  robust  with  crustal 
thicknesses  up  to  2  km  greater  than  at  nearby  a  vol¬ 
canic  segments  [Jokat  and  Schmidt- Attrsch.  2007], 
These  results  indicate  that  the  relationship  between 
spreading  rale  and  magma  supply  is  not  systematic 
on  the  scale  of  individual  ridge  segments,  and  that 
volcanic  processes  are  highly  heterogeneous  at 
ultraslow  spreading  rates. 

[-]  Spreading  rates  along  the  global  MOR  system 
range  from  ultrafasl  (>160  mm  yr  *.  full  rale)  to 
ultraslow  (<12  mm  yr  *).  Fast  spreading  ridges  are 
believed  to  experience  fr  equent  (order  10-100  years), 
small-volume,  volcanic  eruptions  aloug  a  narrow, 
axial  neovolcanic  zone  that  reflect  a  relatively  liigh 
magma  supply  rate  [e.g..  Macdonald.  1982;  Perfit 
and  Chadwick.  1998;  Sinton  et  al..  2002;  Soule  et 
al..  2007],  Morphologically,  fast  spreading  ridges 
exhibit  an  axial  high  with  a  narrow,  scmi-conlinuous, 
graben  [Small.  1994;  Fornari  et  al.,  2004;  Soule 
et  al.,  2009],  By  contrast,  at  slow-spreading  ridges 
the  dominant  morphologic  feature  is  a  riti  valley, 
which  is  believed  to  reflect  relatively  low  and  inter¬ 
mittent  magma  supply  rates  combined  with  increased 
tectonic  crustal  extension  [Smith  and  Cann.  1993; 
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Simon  ami  Detrich,  1992;  Cmm  am I  Smith,  2005; 
Rubin  ami  Sinton,  2007;  Bella  anti  Ito,  2008].  Vol¬ 
canic  eruptions  at  slow  spreading  ridges  are  thought 
to  be  infrequent  (order  1000-10,000  years),  and  to 
produce  discrete,  point-source  volcanic  constructions 
across  the  rift  valley  floor  [Smith  anti  Conn,  1993; 
Cann  and  Smith ,  2005;  Cochran,  2008;  Searle  el  al.. 
2010], 

|s|  Given  the  rilt  valley  morphology  and  intermit¬ 
tent  and  likely  reduced  magma  supply  for  volcanic 
processes  al  ullraslow-spreading  ridges,  magma 
chambers  are  thought  to  lie  beneath  thicker,  stron¬ 
ger,  lithospheric  lids,  and  therefore  require  larger 
magma  pressures  to  induce  emption.  Thus  it  has 
been  hypothesized  that  ultraslow  spreading  ridges 
experience  large  (1-10  km')  and  infrequent  (order 
1 0.000- 1 00,000  years )  volcanic  eruptions  and.  more 
generally,  that  eruption  volume  scales  inversely  with 
spreading  rate  [Perfit  and  Chadwick.  1998;  Sinton 
et  al.,  2002],  Recently  the  deep-submergenee  tech¬ 
nology  to  allow  for  direct  visual  observations  of  the 
Gakkel  Ridge  and  detailed  mapping  of  lava  flows 
beneath  the  perennial  pack  ice  has  made  it  possible  to 
test  these  hypotheses.  In  2007  the  AGAVE  expedi¬ 
tion  utilized  a  novel  vehicle  (CAMPER)  to  acquire 
high-definition  digital  imagery  and  precisely  located 
rock  samples  from  the  Gakkel  Ridge  at  the  7CE  and 
85  E  sites  beneath  the  pack  ice.  Initial  results  from 
the  85  E  site  demonstrated  that  large  portions  of  the 
seafloor  in  this  area  are  blanketed  with  pyroclastic 
material  from  explosive  eruptions  [Si dm  et  al., 
2008],  We  expand  on  this  result  by  using  the  com¬ 
plete  suite  of  dive  images  to  constrain  volcanic 
processes  at  the  scale  of  individual  lava  flows.  Our 
dive  imagery  demonstrates  that  the  seafloor  in  this 
area  is  covered  by  a  large  number  of  small  lava 
flows,  as  opposed  to  a  single,  large-volume  flow, 
which  is  at  odds  the  hypothesis  that  eruption  volume 
varies  inversely  with  spreading  rate  [Perfit  and 
Chadwick.  1998;  Sinton  et  al..  2002],  We  use  our 
dive  images  to  develop  a  new  model  for  volcanic 
processes  at  the  85'  E  segment,  and  discuss  the 
implications  of  this  model  for  ullraslow-spreading 
volcanic  processes,  in  general. 

2.  Gakkel  Ridge  and  the  85 "E  Site 

|f,|  Dredging  studies  along  the  Gakkel  Ridge  have 
recovered  a  mix  of  ultramafic  and  mafic  lithologies 
reflecting  spatial  and  temporal  variations  in  magma 
supply  [Dick  et  al.,  2003;  Michael  et  al.,  2003]. 
Dredging  from  the  AMORE  expedition  in  2001  led 
to  the  division  of  the  Gakkel  Ridge  into  three 


distinct  tectono-magmatic  provinces:  the  Western 
Volcanic  Zone  (WVZ)  from  7  W  to  3  E,  the 
Sparsely  Magmatic  Zone  (SMZ)  from  3  to  29- E. 
and  the  Eastern  Volcanic  Zone  (EVZ)  from  29  to 
96  E  [Michael  et  al.,  2003]  (Figure  la).  The 
AMORE  expedition  dredged  almost  exclusively 
basalt  lithologies  from  within  the  axial  valley  of  the 
Western  Volcanic  Zone,  and  aeromagnetic  studies 
revealed  a  continuous  magnetic  anomaly  along  the 
rift  valley  [Jokat  and  Schmid t-Aursch,  2007],  con¬ 
sistent  with  sustained  volcanic  accretion  for  the  past 
at  least  25  Myr  [Jokat  et  al.,  2003].  By  contrast, 
dredges  acquired  from  lire  Sparsely  Magmatic  Zone 
contained  a  mixture  of  basalt  and  peridotile  litholo¬ 
gies,  indicating  that  in  some  places  the  magma  sup¬ 
ply  is  insufficient  to  sustain  volcanic  accretion, 
leading  to  the  exhumation  of  peridotite  lithologies  by 
extensional  faulting.  This  hypothesis  is  consistent  with 
large  (rcc-waicr  gravity  anomalies  (85-150  niGal) 
observed  over  the  ridge  axis  in  this  zone  [Coaklev 
and  Cochran.  1998],  Along  the  Eastern  Volcanic 
Zone  the  ridge  appears  to  consist  of  at  least  seven 
widely  spaced  (~30  km)  volcanic  centers.  The 
average  depth  of  the  axial  valley  of  the  Gakkel 
Ridge  (4600  m  to  5100  m)  is  deeper  than  those  of 
the  Mid-Atlantic  Ridge  (3500-4000  m)  and  South 
West  Indian  Ridge  (3090-4730  m)  (Figure  2)  [Le 
Douaran  and  Francheteau,  1981;  Sempere  et  al.. 
1990.  1993;  Neumann  and  Forsyth,  1993:  Cannat 
et  al.,  1999],  At  the  easternmost  end  of  the  ridge 
the  axial  valley  becomes  buried  in  sediments  from 
the  Siberian  margin,  and  gravity  data  suggest  the 
crust  becomes  ‘vanishingly’  thin  [Coaklev  and 
Cochran,  1998]. 

|i|  The  85  E  segment  of  the  Gakkel  Ridge  is  of 
particular  interest  because  in  1999  it  hosted  the 
largest  earthquake  sw'ann  delected  to-date  by  the 
Global  Seismic  Network  on  the  spreading  axis  of  a 
MOR  [Midler  and  Jokat,  2000;  Tolstoy  et  al.,  2001], 
Tw'o  hundred  and  fifty  two  earthquakes  with  mag¬ 
nitudes  ranging  from  3.2  to  5.8  were  detected  during 
a  seven-month  swarm,  and  many  events  during  the 
first  few  months  exhibited  non-double  couple  focal 
mechanism  solutions,  which  arc  believed  to  indicate 
volcanic  activity  [Muller  and  Jokat,  2000;  Tolstoy 
et  al..  2001;  Korger  and  Schlindwein,  2012].  Geo¬ 
physical  surveys  of  the  85  E  segment  from  a  nuclear 
submarine  (SCICEX  expedition)  following  the  onset 
of  the  swarm  (summer  of  1999)  imaged  a  large 
region  of  high  acoustic  backscatter  covering  the 
entire  axial  valley  (Figure  3),  which  was  interpreted 
to  represent  multiple  unsedimented,  voluminous  lava 
flow's  [Edwards  et  al.,  2001].  The  AMORE 
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Figure  1.  Study  Area  at  85°E,  Gakkel  Ridge,  (a)  The  Gakkel  Ridge  extends  1 800  km  across  the  Arctic  Basin.  A  star 
denotes  our  study  area  within  the  Eastern  Volcanic  Zone  (EVZ).  (b)  A  bathymetric  chart  of  the  axial  valley  of  the  Gakkel 
Ridge  at  85°E.  Depths  range  from  2543  m  at  the  rift  valley  walls  to  435 1  m  on  the  axial  valley  floor.  Axial  ridges  and 
volcanic  cones  lie  on  two  axis-parallel  lineations  of  focused  volcanism  in  the  center  of  the  rift  valley.  Black  tracks 
denote  the  paths  of  near-seafloor  CAMPER  vehicle  deployments.  Circles  denote  the  locations  of  samples  of  volcani- 
clastic  deposits  and  the  Figure  6  mosaic. 


expedition,  which  returned  to  the  85°  E  site  in  2001, 
found  evidence  for  a  hydrothermal  ‘megaplume’  in 
the  water  column  over  the  axial  valley  [ Edmonds  et 
al.,  2003]  and  seismometers  deployed  on  ice  floes 


over  the  site  recorded  explosive  reports  inteipreted 
to  represent  ongoing  volcanic  activity  [Schlindwein 
el  al,  2007;  Schlindwein  and  Riedel,  2010].  All  of 
this  evidence  led  to  the  hypothesis  that  the  85°  E 


34 


P0NIBR1AND  ETAL.;  GAKKEL  R1DGU  VOLCANISM.  85  U 


1 0 . 1029/20 1 2GC004 1 87 


Geochemistry  , 
Geophysics  f  < J 
Geosystems  r 


-  EAST  PACIFIC  PISE  13'N 

-  MID-AT ANTIC  RIDGE  2 7'N 

-  GAKKEL  RIDGE  85‘E 
SOUTHWEST  INDIAN  RIDGE  65‘E 

Figure  2.  Profiles  of  the  fast  spreading  (Hast  Pacific 
Rise),  slow  spiralling  (Mid-Atlantic  Ridge),  and  ulffiislow- 
spreading  ridges  (Southwest  htdtan  Ridge,  Gakkel  Ridge) 
show  how  axial  motphology  and  depth  vary  w  ith  sprcad- 
ing  rate. 


segment  has  been  undergoing  a  phase  of  active 
volcanism  wherein  it  has  produced  frequent  and 
voluminous  lava  flows  and  vigorous,  sustained, 
hydrothermal  discharge  [Muller  and  Jokal,  2000; 
Edwards  et  al,  2001;  Tolstoy  el  al.  2001;  Edmonds 
el  al ,  2003;  Maker  et  al. ,  2004] 

[s]  The  AGAVH  expedition  returned  to  the  85" K 
site  in  2007  with  a  towed  camera  and  sampling 
system  (CAMPER)  capable  of  returning  high- 
definition  video  and  still  imagery  of  die  seafloor 
beneath  die  ice  cap.  This  imagery  revealed  that 
major  portions  of  die  axial  valley  are  blanketed 
with  volcaniclastic  deposits  interpreted  to  result 
from  magma  fragmentation  during  explosive  erup¬ 
tions  [Sohn  et  al..  2008],  The  fragmentation  mech¬ 
anism  for  the  volcaniclasts  is  enigmatic  because: 
(a)  mid-ocean  ridge  basalts  typically  exhibit  vola¬ 
tile  contents  well  below  what  is  necessary  for 
magmatic  fragmentation  [Wright,  1994;  Clague 
el  al.  2000;  Head  and  Wilson.  2003],  (b)  hydro¬ 
static  pressures  on  the  axial  valley  floor  at  an 
average  depdi  of  3870  m  are  well  above  die  critical 
point  for  steam  formation,  which  appears  to  rule 
out  hydrovolcanic  fragmentation  (i.c..  interaction 
with  seawater)  [Dries ner  and  Geiger,  2007],  and 
(c)  fragmentation  due  to  cooling  (quench  granula¬ 
tion)  is  inconsistent  w'idi  the  seafloor  distribution 
and  moqihology  of  die  clasts.  These  considerations 
led  Solw  et  al.  [2008]  to  propose  that  die  volcani¬ 
clastic  material  was  fragmented  by  magmatic 
volatiles  (e.g.,  C02)  sequestered  in  large  crustal 
melt  bodies.  Melt  inclusions  trapped  in  volcani¬ 
clastic  glasses  recorded  high  (TT  concentrations 
(167-1596  ppm)  consistent  with  entrapment  of 


COj  saturated  magma  at  depths  >9  km  in  die  crust 
[S/imv  el  al,  2010]  Similar  interpretations  have 
also  been  suggested  for  volcaniclastic  deposits  at 
other  deep  oceanic  spreading  centers  [Clague  el  al, 
2003,  2009]  and  diese  are  among  a  growing  num¬ 
ber  of  observations  attributing  fragmentation  to 
magmatic  volatiles  in  submarine  volcanic  systems 
[Wright,  1994;  Clague  et  al,  2000,  2003,’ 2009; 
Chadwick  et  al,  2008;  Schipper  and  White.  2010; 
Iieardorff  et  aL,  201 1;  Hein,  et  al. ,  201 1;  Resing 
el  al..  201 1], 

3.  Data  and  Methods 
3.1.  Bathymetry 

[«]  In  2007  we  conducted  a  bathymetric  survey  of 
die  axial  valley  and  northern  rift  valley  wall  of  the 
85'F.  volcanic  center  using  I/B  Oden 's  Kongsberg 
EM  120  1°  x  1°  multibeam  sonar.  Badiymetric 
data  were  collected  while  die  ship  drifted  in  the 
pack  ice  (with  engines  silent)  along  overlapping 
tracks  at  the  85*0  segment.  The  resulting  low- 
noise.  high-density  multibeam  soundings  allowed 
us  to  construct  a  30  m  resolution  bathymetric  grid 
of  the  study  area  (Figure  I  b).  Ihis  new  badiymetric 
grid  was  merged  with  die  International  Badiymetric 
Chart  of  the  Arctic  Ocean  (IBCAO)  Version  2.0 
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Figure  3.  SCICEX  1999  acoustic  reflectivity  map. 
A  SCICEX  1999  acoustic  reflectivity  map  that  encom¬ 
passes  our  study  area  (yellow  box)  shows  an  —30  km 
long  region  of  high  reflectivity  [Ettwanls  et  al. ,  2001], 
Our  study  area,  including  Jessica’s  Hill  (JH),  Duquc's 
Hill  (1)H ),  Oden  Volcano,  and  Ixikc  Volcano  fall  w  ithin 
the  zone  of  high  reflectivity. 
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[Jakobsson  et  al. ,  2008].  IBCAO  2.0  was  in  turn 
compiled  using  all  available  bathymetric  data  in  the 
Arctic  Ocean  at  the  time  of  compilation,  including 
soundings  collected  in  the  Gakkcl  Ridge  area  dur¬ 
ing  the  SCICEX  1999  and  A  MORI:  cruises.  These 
data  are  available  through  Marine  Geoscience  Data 
System,  the  NOAA  IBCAO  website,  and  al  http:// 
oden.geo.su.se/AGAVE2007.php  [Jakohsson  elal., 
2008;  Ryan  et  al..  2009], 

3.2.  Seafloor  Imagery,  Seawater 
Measurements,  and  Rock  Samples 

|  in)  Near-bottom  imaging  surveys  were  conducted 
with  the  deep-lowed  CAMPER  vehicle,  a  towed, 
fiber-optic  vehicle  that  operates  at  altitudes  of  1 .5  to 
3  m  above  the  seafloor.  The  CAMPER  vehicle  col¬ 
lected  digital  video  and  still  imagery,  measurements 
of  near-bottom  water  properties,  and  geologic  and 
biological  samples  on  thirteen  deployments  at  the 
85  E  site  ( Figure  1  b).  Because  the  vehicle  is  operated 
from  an  idle  ship  within  the  pack  ice,  its  trajectory 
over  the  bottom  is  controlled  by  the  direction  and 
flow  of  the  local  ice  pack,  with  typical  drift  speeds 
oT  0.15  to  0.3  knots.  The  position  of  CAMPER  on 
the  seafloor  was  determined  from  ship  GPS  navi¬ 
gation,  collected  for  ship  heading,  and  assuming  an 
essentially  vertical  wire  from  the  ship  to  the  towed 
vehicle  [cf.  Stranne  et  a!.,  2010]. 

1 1 1 )  The  multipurpose  vehicle  was  outfitted  with 
five  standard  and  high-definilion  video  cameras  that 
streamed  real-time  video  to  the  ship  via  a  fiber-optic 
cable.  Conduetivity-Temperalure-Deplh  (CTD)  and 
an  Eh  sensor  for  meastuing  the  oxidation-reduction 
potential  of  seawater  were  also  mounted  on  the 
vehicle  in  order  to  identify  the  presence  of  reduced 
(i.e..  fresh)  hydrothermal  fluids.  The  Eh  values 
presented  in  this  paper  are  the  measured  Pt  electrode 
voltage  against  the  Ag/AgCl  reference  electrode  in 
saturated  K.CI  solution.  Eh  measurements  can  be 
converted  to  the  redox  potential  of  seawater,  i.e., 
against  the  standard  or  normal  hydrogen  electrode 
(SHE  or  NIIE)  by  applying  a  temperature  correc¬ 
tion,  but  we  do  not  make  the  conversion  because 
the  electrodes  on  the  moving  CAMPER  platform 
never  reach  equilibrium  with  seawater  [Nakamura 
et  al..  2000],  Measured  thermal  anomalies  repre¬ 
sent  potential  temperature  anomalies  derived  from 
absolute  temperature,  pressure  and  salinity  mea¬ 
surements  adjusted  for  an  adiabatic  displacement  to 
a  reference  pressure  of  3000  dliars  [Bryden.  1973; 
Fofonojf  and  Millard ,  1983;  Baker  et  al..  1994]. 
Approximately  4  kg  of  unconsolidated  volcani- 
clastic  deposits.  Sample  12-1  and  Sample  13-1, 


were  sampled  from  two  sites  on  the  axial  valley 
seafloor  at  4000  m  and  4116  m  water  depths 
respectively  using  a  suction  sampler  on  CAMPER 
(Figure  I  b).  Larger  rock  samples  were  also  grabbed 
from  7  locations  [Shaw  et  al..  2010]. 

3.3.  Volcanic  Deposits  at  die  85°E  Site 

[  is]  We  logged  seafloor  volcanic,  tectonic,  and 
biological  features  in  approximately  55  h  of  I  ID 
seafloor  footage  over  17.5  linear  kilometers  of  sea¬ 
floor  co-registered  with  high-resolution  bathymetric 
data  (~5  km  x  10  km  area.  Figure  lb).  Observa¬ 
tions  were  made  in  seven  primary  categories,  each 
with  multiple  subcategories;  lava  surface  morphol¬ 
ogy,  coloration,  ornamentation,  brokcn/intact  lava 
surfaces,  volcaniclastic  deposit  thickness,  microbial 
mat  coverage,  and  prescncc/absenec  of  macro¬ 
biology  (Figures  S1-S7  in  the  auxiliary  material).1 
These  observations  allowed  us  to  identify  lava  How 
contacts  and  How  directions,  relative  eruption  ages, 
possible  source  vents,  and  relationships  between 
volcanic  processes,  hydrothermal  discharge,  and 
biological  colonization. 

|n|  We  categorized  lava  surface  morphologies  as 
pillow,  lobate.  or  sheet  flow  [e.g..  Ballard  et  al., 
1975,  1979;  Lonsdale,  1977a,  1977b;  Bonatti  and 
Harrison,  1988;  Fundis  et  al.,  2010],  which  arc 
thought  to  rcllcct  increasing  volume  flux  [Fink  and 
Griffiths,  1992;  Gregg  and  Fink.  1995].  In  addition, 
we  logged  the  orientation  of  lineations,  folds,  or 
pillow/lobalc  elongation,  which  can  be  related  to  the 
direction  of  lava  transport  [e.g.,  Chadwick  et  a!., 
2001;  Garry  et  al.,  2006],  We  logged  lava  colora¬ 
tions  (black/glassy  Versus  orange/oxidized),  and  the 
presence  or  lack  of  ornamentation  (delicate  fingers 
or  buds),  as  indicators  of  relative  age.  We  did  not 
detect  any  pelagic  sediment,  which  is  commonly 
used  as  an  indicator  of  relative  lava  age  in  mid- 
ocean  ridge  environments  [e.g.,  Cann  and  Smith, 
2005],  in  our  study  area.  Low  biological  produc¬ 
tivity  beneath  the  sea  ice  limits  sedimentation  rates 
but  pelagic  sediment  has  been  observed  in  the 
region  [Vogt  el  al..  1979;  Thiede  et  al..  1990], 
Recent  volcanism  may  overprint  pelagic  sediments. 

[u]  We  categorized  volcaniclastic  deposit  cover¬ 
age.  which  was  present  in  all  imagery,  as  light, 
medium,  or  thick.  ‘Light'  coverage  represents  lava 
surfaces  that  are  visible  through  a  light  "dusting”  of 
volcaniclastic  material.  'Medium'  coverage  repre¬ 
sents  surfaces  that  are  not  visible  through  the 


’Auxiliary  materials  ate  available  in  the  HTML.  doi;lO.I029i 
21)1 2CK'004 187. 
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volcaniclaslic  material,  with  a  deposit  thickness  of  a 
few  centimeters.  ‘Thick'  coverage  represents  sur¬ 
faces  where  volcaniclasts  pile  on  lava  surfaces  with 
a  thickness  of  up  to  ~10  cm,  and  fill  in  void  spaces 
between  seafloor  features  with  tens  of  centimeters 
of  relief.  Potential  controls  on  submarine  volcani¬ 
claslic  deposit  thickness  include  the  age  of  the 
underlying  surface  and  the  proximity  to  the  source 
regions  (e.g.,  Barreyre  et  al.,  2011],  Volcaniclaslic 
material  was  sampled  from  two  sites  (Figure  lb),  as 
described  in  section  4.2.2. 

|  is]  We  logged  the  presence  or  absence  of  biolog¬ 
ical  organisms  and  whether  the  observed  organisms 
are  motile  or  sessile.  Sessile  organisms,  which  arc 
fixed  to  the  seafloor,  included  sponges  and  ane¬ 
mones.  The  presence  or  absence  of  microbial  mats 
was  also  logged.  Motile  organisms,  which  are  free 
to  move  to  available  food  sources,  include  fish,  an 
octopus,  shrimp,  and  amphipods.  The  presence  of 
sessile  organisms  provides  some  indication  of  rel¬ 
ative  seafloor  (i.e.,  eruption)  age  because  it  is  set  to 
zero  when  a  surface  is  repaved  during  an  eruption, 
and  grows  over  time  as  an  area  is  colonized  [Shank 
et  al.,  1998].  The  presence  of  microbial  mats  can  be 
used  as  a  proxy  for  identifying  zones  of  active 
hydrothermal  discharge  because  these  organisms 
derive  metabolic  energy  by  oxidizing  hydrothermal 
fluids  [Helntke  et  al.,  2007], 

4.  Results 

4.1.  Seafloor  Geology 

|  if.]  The  85  h  segment  of  the  Gakkel  Ridge  has  a 
12-15  km  wide  and  1.8  km  deep  rill  valley  mor¬ 
phology  with  an  average  valley  floor  depth  of 
3870  m,  and  a  maximum  depth  of  4351  m.  Linear 
volcanic  ridges  and  cones  (many  of  which  have 
large  central  craters)  are  the  predominant  construc¬ 
tional  volcanic  features  inside  the  axial  valley  and 
rise  several  hundred  meters  above  the  valley  floor. 
These  features  arc  concentrated  along  two  axis- 
parallel  lineaments  located  in  the  central  portion 
(5  km)  of  the  axial  valley,  with  the  notable  excep¬ 
tion  of  a  large  (1250  m  diameter)  cone  located  on  an 
apparent  fault  terrace  about  halfway  up  the  northern 
rill  valley  wall  (Figure  lb).  Only  12  faults  and  fis¬ 
sures  are  observed  along  20  linear  km  of  track 
(Figure  7).  No  peridotite  lithologies  are  observed  on 
the  seafloor. 

[ it]  Tlte  near-bottom  CAMPER  transects  provide 
coverage  over  portions  of  the  two  intravalley  vol¬ 
canic  lineaments,  including  Oden  and  Loke  volca¬ 
nic  cones  and  Jessica’s  Hill  and  Duque's  Hill 


volcanic  ridges.  Oden  and  Loke  arc  located  on  the 
southern  lineament;  Jessica’s  Hill  and  Duque’s  Hill 
are  located  on  the  northern  lineament  (the  linea¬ 
ments  are  offset  by  2-3  km,  Figure  lb).  The  lava 
morphology  (lobate  vs  sheet),  coloration  (orange/1 
altered  vs  glassy),  and  volcaniclaslic  deposit  thick¬ 
ness  observed  along  the  CAMPER  tracklines  are 
shown  in  Figure  4.  Talus  and  pillow  lavas  arc  the 
dominant  deposit  types  observed  in  the  study  area. 
Pillow  lavas  are  observed  along  44%  of  the  track¬ 
line.  Some  lobate  and  sheet  flows  are  observed  on 
the  volcanoes  (14%  of  the  trackline)  and  ridgclines 
(<1%)  (Figures  S1-S7  in  the  auxiliary  material). 
Observations  of  glassy  (i.e..  young)  lavas  are 
restricted  to  the  Oden  and  Loke  Volcanoes,  where 
higher  effusion  rate  sheet  flows  are  also  found.  The 
linear  ridges  appear  to  be  older  volcanic  features 
consisting  largely  of  oxidized  and  broken  pillow 
basalts.  Seafloor  talus  blocks  cover  the  steep  slopes 
and  the  bases  of  fault  scarps.  Glassy,  unconsoli¬ 
dated,  volcaniclastic  material  was  observed  on  all 
surfaces  regardless  of  apparent  age  or  morphology 
(96%  of  trackline).  Volcaniclastic  deposits  are 
thickest  on  the  Duque’s  Hill  and  Jessica’s  Hill  lin¬ 
ear  ridges,  thinnest  on  the  Oden  and  Loke  volcanic 
cones,  and  are  present  as  a  light  ‘dusting’  on  the 
youngest  (i.e.,  glassy)  lavas.  Microbial  mats  arc 
observed  over  5.5%  of  the  trackline. 

4.1.1.  Oden  Volcanic  Cone 

|is|  Oden  is  one  of  several  volcanic  cones 
(including  Loke  and  Thor)  that  constitute  the  SE 
end  ol'the  southern  volcanic  lineament  (Figure  lb). 
Three  of  the  CAMPER  deployments  pass  over  the 
Oden  volcanic  cone,  and  a  fourth  passes  just  to  the 
east  (Figure  lb),  providing  some  of  the  densest 
track  line  coverage  in  our  survey  area.  The  flat- 
topped  cone  is  1400  nr  in  diameter  and  170  m  tall, 
with  a  600  m  wide  by  70  m  deep  central  crater 
(Figure  5).  Our  deployments  pass  over  the  edges  of 
the  crater,  and  orange-yellow  microbial  mats  blan¬ 
ket  much  of  the  seafloor  in  this  area  (Figure  5c). 
Where  visible,  the  seafloor  wtthin  the  crater  is 
covered  with  talus  blocks  and  a  thin  (<1  cm)  coat¬ 
ing  of  glassy  volcaniclastic  material  (Figure  5d). 

(i«|  On  the  crater  rim  we  observe  mostly  pillow 
lavas  and  lobate  flows,  including  obvious  contacts 
between  older  flows  and  fresh,  glassy  surfaces 
(Figures  5e  and  5f ).  Some  of the  lavas  with  orange, 
altered,  surface  colorations  have  retained  delicate 
ornamentations,  suggesting  that  coloration  does  not 
necessarily  correlate  with  eruption  age.  We  also 
observed  a  glassy,  ropy  sheet  flow  with  a  flow 


37 


juMS-  Geochemistry  _  ,  \ 
W  '  Geophysics  tf 
Geosystems 


PONTQR1AND  I  T  AL.:  GAKKEL  RIDGE  VOLCANISM.  85  E 


1 0. 1 029/20 12GC004 187 


A.  LOBATE  (•)  AND  SHEET  FLOWS  (♦)  C.  VOLCANICLASTIC  DEPOSIT  THICKNESS 


85'39' 


85'36’ 


85'00'  85'30' 


Figure  4.  Dislribulions  oflavas  and  volcaniclasts.  (a)  A  bathymetric  map  shows  that  lobate  flows  (dots)  are  observed 
on  all  four  seafloor  features  examined  in  digital  imagery.  Lobate  flows  (e.g..  Figure  40  are  found  at  the  ridge  crests  of 
Jessica's  Hill  and  Duque’s  Hill  axial  volcanic  ridges.  They  arc  also  abundant  on  Oden  Volcano  and  the  northern  apron 
of  Loke  Volcano.  Sheet  flows  (diamonds;  e.g..  Figure  4e),  products  of  more  effusive  eruptions,  are  also  abundant  on 
the  volcanic  cones,  (b)  Oxidized  or  weathered  pillow  flows  (dots)  arc  observed  on  all  seafloor  features,  as  well  as  the 
seafloor  in  the  depression  between  the  two  axial  valley  trends.  Glassy,  fresh,  and  highly  ornamented  lavas  (diamonds; 
e.g..  Figure  4d)  are  only  observed  on  the  steep  slopes  and  summits  of  Oden  and  Loke  Volcanos  and  represent  the 
youngest  eruptive  products  in  the  study  area,  (c)  Volcaniclastic  deposits  of  variable  thickness  are  observed  throughout 
the  study,  (d  g)  Coverage  from  a  light  dusting  to  a  thick  cover  (~  1 0  cm)  that  buries  lava  surfaces  and  fills  in  crevices 
between  rocks  on  the  seafloor.  Duque’s  Hill  has  the  thickest  volcaniclastic  cover.  “Fresh"  lavas  on  Oden  and  Loke 
Volcanoes  are  lightly  dusted  with  volcaniclastic  material. 


orientation  radial  to  the  center  of  the  volcano  on  the 
crater  rim. 

l»J  The  outer  walls  of  the  cone  consist  largely  of 
pillow  lubes  elongated  in  the  downslope  direction, 
lobate  flows,  and  talus.  In  some  places  the  pillow 
tubes  are  glassy  and  ornamented,  in  other  places  they 
have  orange  surfaces.  The  talus  slope  on  the  NE  wall 
is  interrupted  by  an  orange-stained  lobate  flow  (with 
a  thin  microbial  coating).  The  abrupt  nature  of  this 
contact  and  the  lack  of  talus  on  the  lobate  flow  indi¬ 
cates  that  the  flow  is  a  relatively  young  feature  that 
erupted  from  within  the  talus  field. 


[:i  |  The  marginal  regions  around  Oden  exhibit  a  mix 
of  orange-stained  flows  and  fresh,  glassy  lavas.  The 
4000  m  deep  bathymetric  saddle  between  Oden  and 
Duque's  Hill  contains  weathered,  unornamented, 
pillow  lavas  and  talus  breccia,  with  no  evidence  of 
recent  volcanic  activity.  On  the  4100  m  deep  plat¬ 
form  between  Oden  and  Tit  or  we  observe  pillow 
lavas  and  an  ornamented,  orange-stained,  lobate 
flow.  Glassy,  elongate,  ornamented  pillow  lavas 
were  observed  at  the  base  of  the  Thor  volcanic  cone. 
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Figure  5.  (a)  Oden  Volcano  Map  view  and  (b)  profile  of  Oden  Volcano — a  1600  m  diameter  volcano  with  ~200  m 
of  relief  and  an  80  m  deep  x  700  m  wide  circular  crater,  (c)  Talus  and  microbial  mats  in  the  crater,  (d)  Talus  ejecta  on 
the  flank  of  the  volcano  covered  in  a  moderate  thickness  of  large  pyroclasts.  (e)  Fresh,  glassy  lobate  flow  with  evi¬ 
dence  of  hydrothermal  alteration  at  edges,  (f)  Weathered,  elongate  pillow  flow  with  moderate  pyroclastic  coverage. 


4.1.2.  Loke  Volcanic  Cone 

[22 1  The  Loke  volcanic  cone  is  1200  m  in  diameter 
with  150  m  of  vertical  relief  from  a  basal  depth  of 
more  than  4000  m.  One  CAMPER  track  line  passed 
over  the  SE  side  of  the  structure,  and  two  more 
surveyed  the  area  immediately  NE  of  the  cone.  The 
southeastern  wall  of  the  cone  appears  to  have  been 
breached  by  a  lava  flow  and/or  collapsed  such  that 
the  volcanic  edifice  is  horscshoc-shapcd  (Figure  lb). 
The  bathymetric  “apron”  that  extends  outwards  from 
the  volcano  is  comprised  of  multiple,  low-relief,  lava 


flows  with  obvious  contacts  that  extend  radially  from 
the  base  of  the  cone.  At  the  edge  of  the  apron,  over 
500  m  from  the  slopes  of  Loke,  we  found  black, 
glassy  pillows.  Within  the  apron  we  observed 
a  >50-m-wide  lava  channel  that  displayed  typical 
morphologies  of  lineated  sheets  at  the  channel  center 
surrounded  by  chaotic,  jumbled  textures  and  bor¬ 
dered  by  lobates  (Figure  6)  [Soule  et  al.,  2005;  Garry 
et  al.,  2006],  Sheet  flows  had  hackly,  chaotic,  line¬ 
ated,  and  ropy  surface  textures.  On  the  volcanic  cone 
some  glassy  lava  surfaces  (Figure  4d)  were  overlain 
by  microbial  mats. 
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Figure  6.  Sheet  flow  on  Lokc's  apron  Mosaic  of  a  sheet  flow  on  Lokc’s  apron  (see  Figure  lb).  Altered  lavas  covered 
in  a  moderate  thickness  of  volcaniclastic  deposits  show  a  transition  from  chaotic  sheet  flow  to  lincatcd  sheet  flow, 
boarded  by  lobates.  The  white  arrow  in  Figure  fib  denotes  the  approximate  orientation  of  flow. 


4.1.3.  Jessica’s  Hill  and  Duque’s  Hill 
Linear  Ridges 

[23]  Jessica's  Hill  is  5.5  km  long, 1.8  km  wide,  and 
0.4  km  high  with  a  peak  at  3650  mbsl;  Duque’s  Hill 
is  2.5  km  long,  1.2  km  wide,  and  0.125  km  high 
with  a  peak  at  3820  mbsl.  Both  of  these  features  are 
located  on  the  northern  volcanic  lineament  in  the 
rift  valley,  and  appear  to  be  volcanic  constructional 
features  similar  to  axial  volcanic  ridges  described 
on  the  slow-spreading  Mid-Atlantic  Ridge  [ Ballard 
and  van  Andel,  1977;  Ramberg  and  van  Andel, 
1977;  Smith  and  Cann.  1990,  1992;  Searle  et  a!., 
2010].  The  ridges  consist  primarily  of  orange- 
stained.  broken,  pillow  basalts  and  talus,  and  we 
did  not  observe  any  evidence  of  fresh  eruptive 
products  (i.e.,  glassy  surfaces)  on  these  features. 
Some  pillows  were  ornamented,  but  all  had  dull  or 
oxidized  surfaces.  The  steep  slopes  of  the  ridges  are 
characterized  by  elongate  pillow  lavas,  and  steep, 
9-1 0  m  high  scarps  consisting  of  fractured  pillows 
were  observed  in  six  places.  Both  Jessica’s  Hill  and 
Duque’s  Hill  are  covered  with  moderate  to  thick 
layers  of  volcaniclastic  deposits,  with  the  thickest 
coverage  on  Duque’s  Hill  where  the  volcaniclasts 
completely  bury  the  lavas  and  talus  in  many  places 
(Figure  4g). 

4.2.  Volcaniclastic  Deposits 

4.2.1.  Deposit  Thickness 

[24 1  Our  qualitative  evaluation  of  volcaniclastic  layer 
thickness  varied  from  a  light  dusting  (Figure  4d)  to 
up  to  ~10  cm  thick  cover  (Figure  4g)  along  the 
seafloor  tracks  (Figure  4c).  Volcaniclastic  deposits 
were  thickest  on  the  linear  ridges  and  in  the  depres¬ 
sion  between  the  two  volcanic  lineaments,  and  thin¬ 
nest  on  the  glassy,  unaltered  lava  flows  identified  on 


Oden  and  Loke  Volcanoes.  Abrupt  transitions  in 
volcaniclastic  deposit  thickness  were  observed  on  the 
seafloor  on  the  apron  of  Loke  and  on  the  flanks  ofthe 
volcanoes,  and  oflen  corresponded  with  other  indi¬ 
cators  of  flow  contacts  between  lavas  of  different 
relative  ages  (motphology,  coloration).  Volcaniclasts 
were  present  on  lava  surfaces  of  all  ages  and  on  non- 
volcanic  surfaces  (c.g.,  talus  slopes)  suggesting  that 
they  were  not  derived  from  quench  granulation  of  lava 
flows  [e.g.,  Bariza  et  a!.,  1984],  In  addition,  volcani¬ 
clasts  were  piled  on  the  tops  of  pillows  high  on  rid¬ 
gelines  of  Duque's  Hill  and  Jessica's  Hill,  indicating 
deposit  by  fall  rather  than  flow.  Although  difficult  to 
quantify,  there  appeared  to  be  a  tendency  for  the  clast 
layer  to  thicken  downslope  on  the  flanks  ofthe  linear 
ridges,  suggesting  some  post-depositional  downslope 
movement  of  material. 

1 25 1  In  general,  there  is  a  correlation  between  the 
apparent  (relative)  age  of  the  lava  flows  and  the 
thickness  of  the  volcaniclastic  deposits,  with  the 
freshest-looking  lavas  being  covered  with  a  light 
dusting  of  material  and  the  oldest-looking  lava 
surfaces  being  covered  with  the  thickest  deposits. 
This  suggests  that  the  volcaniclasts  have  been 
deposited  progressively  over  time,  and  that  they  are 
buried  by  younger  lava  flows.  In  most  places  it  was 
not  possible  to  discern  any  systematic  spatial  var¬ 
iations  in  clast  grain  size,  but  on  the  NE  side  of 
Oden  Volcano  the  clast  grain  size  appeared  to 
increase  toward  the  crater,  and  the  largest  clasts 
observed  during  our  surveys  were  observed  on  the 
crater  rim,  suggesting  that  the  crater  may  have  been 
a  source  vent  for  volcaniclastic  material.  Within  the 
crater,  a  light  dusting  of  volcaniclastic  material  sat 
on,  or  was  embedded  into,  the  microbial  mats,  and 
could  be  observed  on  top  of  talus  blocks  protruding 
from  the  mats  as  well  (Figure  5c). 
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Figure  7.  SH.M  images  of  voleaniclast  morphologies,  (a-c)  Limu  o  Pele.  or  bubble-wall  fragments  and  tluidal 
morphologies  (d-f )  Platy  clasts  with  crystalline  inclusions  and  elongate  vesicles,  (g— i)  Blocky  dense  fragments  dom¬ 
inate  the  samples. 


4.2.2.  Clast  Morphologies 

1 |  We  analyzed  volcaniclaslic  material  from  two 
suction  samples  acquired  during  our  CAMPER 
surveys.  Sample  12-1  was  collected  from  the  sad¬ 
dle  between  Jessica’s  Hill  and  Duque’s  Hill  AVRs 
at  a  depth  of  4000  m,  and  Sample  13-1  was  col¬ 
lected  from  thick  deposits  at  the  base  of  Duque’s 
Hill  at  a  depth  of  4116  m  (Figure  lb).  Grain  size 
analysis  revealed  clasts  from  fine  volcanic  ash 
(<0.149  mm)  to  greater  than  4  mm  in  length.  Mean 
grain  size  was  consistent  with  volcaniclaslic 
deposits  from  the  slow-spreading  Mid-Atlantic 
ridge  [Eissen  el  al.,  2003],  Details  of  grain-size 
sorting  are  described  by  Barrexre  et  al  [201 1], 

[27 1  We  examined  a  representative  sample  of 
~  1 200  volcaniclasts  by  stereo  microscope  and  SEM 
for  morphology  and  other  characteristics  including 
apparent  age  (juvenile  or  non-juvenile),  material 
(glass,  biogenic,  or  silt),  and  morphology  (blocky, 
platy  or  tabular,  or  sheet-like).  Both  samples  consist 
almost  entirely  of  juvenile  clasts  of  glassy  basalt. 
Slightly  vesicular  blocky  fragments  dominate  the 
deposits  (Figures  7g— 7i;  see  also  Figures  S8  and  S9 


in  the  auxiliary  material).  A  moderate  fraction  of 
fragments  (<22  wt.%  and  <12  wt.%  for  Samples 
12-1  and  13-1,  respectively)  are  limu  o  pele  bubble 
wall  fragments  [Mattox  and  Mangan,  1997; 
Matcher  and  White,  2001;  Clague  et  al. ,  2009], 
Some  fragments  preserve  elongate  vesicles  or  lubes, 
also  representative  of  rapid  quenching.  Most  vesi¬ 
cles  imaged  are  not  interconnected,  indicating  that 
the  magma  had  limited  permeability,  which  could 
facilitate  the  development  of  gas  slugs  [Namiki  and 
Manga,  2008;  Polacci  et  al.,  2008;  Schipper  and 
White,  2010], 

[is|  Bubble-wall  fragments  arc  flat,  curved,  or  com¬ 
plexly  folded  (Figures  7a-7c).  We  also  observed 
tubelike  structures  (Figure  7d),  triple  junctions,  elon¬ 
gated  vesicles,  and  surface  textures  preserved  by  rapid 
cooling  (Figure  7a)  [Clague  et  al.,  2009],  Surfaces  of 
limu  o  Pele  fragments  range  from  smooth  to  bumpy, 
the  latter  due  to  small  inclusions  of  crystals  and  bub¬ 
bles.  No  Pele’s  hair  or  rod  shapes  are  evident  in  the 
samples.  Dense,  blocky  volcaniclasts  commonly  dis¬ 
play  conchoidal  fractures  indicative  of  fracture  at 
temperatures  lower  than  the  glass-transition  tempera¬ 
ture  [Maicher  et  al,  2000].  There  is  a  very  small 
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Figure  8.  Near-bottom  seawater  anomalies,  (a)  A  map  shows  locations  of  near-seafloor  seawater  temperature  and 
Eh  anomalies  and  faults  and  fissures.  CTD  data  were  not  collected  on  Jessica’s  Hill.  The  yellow  star  denotes  a  visual 
observation  of  cloudy  fluid  90  in  above  the  seafloor  of  Lokc’s  apron,  (b)  Concurrent  measurements  of  potential  tem¬ 
perature  and  Eh  records  during  deployment  DCS-1 5  on  Oden  show  temperature  and  Eh  anomalies  (diamonds  and  dots, 
respectively)  detected  within  tens  of  meters  of  one  another  on  the  edifice,  (c)  Cloudy  fluid  emanates  from  beneath  an 
overhanging  ledge  on  Oden  Volcano  where  temperature  and  Eh  anomalies  are  detected.  Yellow  microbial  mats  coat 
the  rocks. 


component  of  crystalline  and  altered  basalt  rock  in  the 
samples  but  it  is  not  clear  whether  these  arc  from  the 
volcaniclastic  deposits  or  originated  front  underlying 
intact  lava  flows.  Free  olivine  containing  melt  inclu¬ 
sions  [5/taw  el  al.,  2010]  and  a  small  component  of 
altered  clasts  with  orange  colorations,  which  may  be 
palagonite,  were  also  observed. 

4.3.  Near-Bottom  Water 
Column  Anomalies 

|2»|  As  described  in  section  3.2.  the  CAMPER  vehi¬ 
cle  was  equipped  with  sensors  to  detect  near-bottom 
water  column  anomalies  in  temperature  and  Eh, 
which  can  be  used  to  delect  regions  of  hydrothermal 
discharge  [Baker  el  al.,  1994;  Nakamura  el  al.,  2000], 

|.iu|  Near-bottom  (1-3  m  above  the  seafloor)  tem¬ 
perature  anomalies  (SB  of  0.005-0.0 1 9°C)  were 
observed  over  Duque’s  Hill.  Oden  Volcano,  and 
Loke  Volcano  (Figure  8;  see  also  Table  SI  in  the 
auxiliary  material).  Twenty-three  Eh  anomalies  were 
detected  on  the  flank  of  Oden  Volcano.  No  anoma¬ 
lies  were  detected  on  the  seafloor  between  the  vol¬ 
canic  lineaments. 

[ji]  Many  of  the  Eh  anomalies  measured  on  Oden 
Volcano  were  concurrent  with  observations  of 
dense  microbial  mats  (Figure  9).  One  or  two  near¬ 


bottom  Eh  anomalies  were  detected  over  Duque’s 
Hill,  Jessica’s  Hill,  and  the  apron  of  Loke  Volcano. 
Oden  and  the  apron  of  Loke  arc  the  only  places 
where  temperature  and  Eh  anomalies  were  co¬ 
located.  Within  the  crater  of  Oden  Volcano,  the  tips 
of  orange-stained  talus  blocks  protruded  through 
near-continuous  microbial  mat  (Figure  5c)  and 
seawater  measurements  identified  temperature  and 
Eh  anomalies  within  20  m  of  one  another 
(Figure  8b).  On  the  flanks  of  the  volcano,  in  one 
collapse  pit,  a  pool  of  cloudy  water  was  observed 
(the  only  place  in  our  imagery  where  hydrothermal 
fluids  could  be  observed  discharging  from  the  sea¬ 
floor)  simultaneously  with  a  strong  Eh  anomaly 
(Figure  8c)  and  observations  of  microbial  mats 
(Figure  9).  Loke  Volcano  also  hosted  extensive 
microbial  mats.  On  the  apron  of  Loke  Volcano, 
temperature  and  Eh  anomalies  were  detected  in  the 
vicinity  of  a  fissure  and  cloudy  seawater  was 
observed  at  90  ni  above  the  seafloor  during  ascent  of 
the  CAMPER  vehicle.  On  Duque's  I  fill  no  microbial 
mats  were  observed,  but  near-seafloor  temperature 
and  Eh  anomalies  were  measured  at  several  areas 
along  the  steep  slopes.  We  observe  Eh  anomalies  on 
Jessica’s  Hill  and  south  of  Oden  in  areas  where 
microbial  mats  arc  not  observed,  suggesting  that  Eh 
anomalies  have  a  hydrothermal  source  and  are  not 
generated  by  the  microbial  mats. 
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Figure  9.  Locations  of  sessile  biology,  (a)  A  bathymetric  map  shows  locations  of  microbial  mats  and  sessile  mega¬ 
fauna.  Locations  of  (b)  microbial  mats  correspond  with  areas  of  diffuse  hydrothermal  venting  and  observations  of  the 
youngest  eruptive  products  on  Oden  and  Loke  Volcanoes.  Sessile  macrofauna  such  as  (c)  sponges  and  anemones  are 
abundant  on  the  relatively  older  seafloor  of  the  axial  volcanic  ridges. 


[32]  Temperature  and  Eli  anomalies  (detected  by 
CTD  ‘tow-yos’)  were  ubiquitous  higher  in  the  water 
column  over  wide  portions  of  the  rifl  valley 
| Stranne  et  al,  2010],  and  modeling  of  these  data 
indicated  that  there  are  up  to  six  distinct  sites  of 
hydrothermal  discharge  on  the  seafloor  within  the 
axial  valley.  Inferred  discharge  sites  from  water 
column  data  agree  with  sites  of  our  near-bottom 
anomalies,  and  near-seafloor  measurements  show 
that  water  anomalies  are  detected  over  100's  of 
meters  along  the  volcanic  cones  and  ridges. 

5.  Discussion 

[jj]  We  assign  relative  ages  of  flow  units  using  the 
abundance  of  volcanic  glass,  the  degree  ofalteration 
of  lava  surfaces,  and  the  sediment  cover  |e.g.,  van 
Antlel  anil  Ballard ,  1979J.  These  subjective  criteria 
are  not  easily  quantified  and  are  affected  by  factors 
other  than  age  alone.  For  example  staining  and 
palagonitization  are  useful  indicators  of  age,  but  the 
presence  of  diffuse  hydrothermal  venting  and 
microbial  mats  may  also  alter  lava  surface  colora¬ 
tion.  Although  volcaniclastic  deposition  is  likely  not 
continuous  in  time  or  space,  we  use  abrupt  transi¬ 
tions  in  volcaniclastic  thickness  as  a  useful  tool  to 
support  relative  age  inferences  from  other  observa¬ 
tions.  The  thickness  of  a  submarine  volcaniclastic 
deposit  depends  on  the  age  of  the  lava  flow  on 


which  it  rests,  proximity  to  the  source  vent,  and  to  a 
lesser  degree,  post-depositional  transport  |e.g., 
Barrevre  et  al.,  2011],  The  relative  abundance  oT 
sessile  biology  also  provides  an  indication  of  sea¬ 
floor  age  because  colonization  begins  at  the  time  or 
lava  flow  emplacement  and  appears  to  progress 
predictably  with  time  [Shank  et  al.,  1998J. 

5.1.  Recent  Eruptions  at  85°E 

[34]  Correlated  lava  surface  morphology  transitions 
and  relative  age  indicators  (e.g.,  staining,  presence/ 
absence  of  fine  ornamentation,  volcaniclastic  deposit 
thickness)  suggest  that  the  survey  area  has  experi¬ 
enced  numerous  eruptions  front  volcanic  vents  scat¬ 
tered  throughout  the  axial  valley  (Figure  10).  Lava 
surfaces  on  the  Oden  volcanic  cone  vary  from  glassy, 
sometimes  ornamented  surfaces  to  orange-stained 
(Figures  5e  and  5f).  Although  the  extent  of  surface 
coloration  may  indicate  eruption  age,  we  note  that  the 
presence  of  hydrothermal  anomalies  in  the  area  may 
expedite  the  lava  surface  alteration  process,  making 
pillows  appeal’  “prematurely”  orange  [Stroncik  and 
Schmineke,  2002],  Alternatively,  the  microbial  mats, 
which  appear  similar  to  mats  composed  of  Fe- 
oxidizing  bacteria  on  Loihi  Volcano  |e.g.,  Emerson 
and  Moyer,  2002;  Edwards  et  a!.,  2011]  may  be  the 
primary  cause  of  coloration  [e.g.,  Standigel  and 
Clague,  20 1 OJ.  As  such,  we  are  cautious  when 
interpreting  flow  contacts  and  require  variations  in 
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Figure  10.  Interpretive  map  of  volcanic  features,  small-volume  effusive  lava  flows,  pyroclastic  deposits,  and  possi¬ 
ble  volcanic  vents  extrapolated  from  near-seafloor  observations  and  bathymetry. 


coloration  to  be  accompanied  by  variations  in  mor¬ 
phology,  degree  of  ornamentation,  and  volcani- 
clastic  thickness. 

1 35 )  To  better  constrain  relative  ages  of  eruptions  on 
Oden  Volcano,  we  evaluate  the  distribution  of 


volcaniclastic  material  on  lava  surfaces.  Although 
volcaniclastic  deposition  is  likely  not  continuous  in 
time  or  space,  differences  in  volcaniclastic  thickness 
are  a  useful  tool  to  support  relative  age  inferences 
from  other  observations.  For  example,  orange- 
stained  lava  pillows  on  Oden’s  eastern  rim  arc 
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covered  in  a  moderate  thickness  of  voleaniclaslic 
deposits,  while  the  unstained,  glassy,  ornamented 
pillows  on  Oden’s  southern  slopes  have  only  a  light 
dusting  of  voleaniclaslic  deposits.  Oden’s  northern 
rim  is  also  covered  in  ornamented  pillows  with  a  light 
dusting  of  voleaniclaslic  deposits,  hut  lava  surfaces 
on  the  northern  rim  are  orange-stained,  possibly  due 
to  rapid  alteration  by  interaction  with  diffuse  hydro- 
thermal  fluids  detected  in  the  area.  Abrupt  changes  in 
voleaniclaslic  deposit  thickness  must  represent  con¬ 
tacts  between  flow  units  rather  than  spatial  vanations 
in  distance  from  an  explosive  source,  which  would 
result  in  more  gradual  changes  in  deposit  thickness. 
Thus  Oden  is  constructed  of  products  of  at  least  two 
eruptions  of  distinct  ages,  one  being  very  recent  and 
possibly  after,  or  coincident  with,  the  emplacement 
of  voleaniclaslic  deposits. 

| m |  A  distinct  contact  was  also  observed  between 
the  relatively  older,  moderately  volcaniclast -coveted 
lava  surfaces  on  Oden  and  lightly  dusted,  glassy, 
ornamented  pillows  oriented  downslope  on  Thor 
Volcano,  which  look  similarly  fresh  to  glassy  flows 
on  Oden.  On  Loke  Volcano  fresh,  glassy  lava  sur¬ 
faces  in  conjunction  with  the  light  distribution  of 
volcaniclasls  suggests  that  the  elongate  pillow  lavas 
that  mantle  the  flanks  of  the  volcanic  cone  erupted 
more  recently  and  overlay  the  eruptive  units  on 
the  apron.  The  eruptive  apron  is  composed  of  high 
elTusion  rale  sheet  flows,  is  pervasively  stained,  and 
has  accumulated  more  voleaniclaslic  material  on  its 
surface.  Orientations  of  sheet  flows  suggest  that  lavas 
were  transported  radially  from  the  central  crater  of 
the  cone, 

| sr)  We  infer  from  the  coincidence  of  glassy,  elon¬ 
gate  pillow  lavas  with  ornamental  buds,  and  a  thin 
dusting  of  volcanidastic  deposits  that  Oden,  Thor, 
and  Loke  Volcanoes  are  sites  of  the  most  recent 
eruptions.  The  extent  of  individual  eruptive  units  are 
poorly  constrained,  but  the  identification  of  a  few 
flow  contacts  suggest  aerial  extents  on  the  order  of 
~33%  or  less  of  the  cone  areas  (~3  km2).  Given 
their  thin  volcanidastic  cover  and  fresh,  glassy 
surfaces,  these  young  Hows  on  the  volcanoes  were 
emplaced  slightly  before  or  at  the  time  of  the  most 
recent  volcanidastic  eruption  that  is  presumed  to 
have  occurred  in  1999  [Muller  and  Jokat ,  2000; 
Tolstoy  et  at.,  2001],  Multiple,  small-volume,  erup¬ 
tions  on  these  volcanoes  thus  characterize  the  most 
recent  effusive  volcanic  activity  at  85c  E. 

|j*|  Duque’s  Hill  and  Jessica’s  Hill  are  axial  volca¬ 
nic  ridges  (AVRs)  constructed  from  multiple  effusive 
eruptions  that  appear  to  have  occutred  earl  ier  than  the 
recent  eruptions  on  the  southern  lineament.  Both 


ridges  arc  dominated  by  talus  slopes,  oxidized  pil¬ 
lows,  fractured  seafloor,  abundant  sessile  macro¬ 
fauna,  and  a  lack  of  microbial  mats.  The  absence  of 
sessile  macrolauna  on  Oden  and  Loke  Volcanoes 
and  their  presence  on  Jessica  and  Duque’s  Hills 
suggest  that  eruptions  on  the  AVRs  may  be  older 
than  any  of  the  volcanoes  on  the  southern  lineament. 
Despite  their  apparent  older  age,  few  faults  and  fis¬ 
sures  are  observed,  although  the  presence  of  steep 
talus  slopes  provides  indirect  evidence  of  faulting. 
Jessica’s  Hill  has  some  lobate  Hows  at  the  ridgeline, 
indicating  moderate  ilux  eruptions.  However,  pillow 
eruptions  dominate  the  lava  morphology  of  both 
ridges.  Both  ridges  are  covered  in  a  moderate  to 
heavy  thickness  of  volcanidastic  deposits.  This  is 
consistent  with  older  eruption  ages  and  hence  longer 
periods  of  volcamclast  accumulation  or  proximity  to 
a  source. 

|  w|  A  relatively  younger  age  for  southern  linea¬ 
ment  volcanoes  is  supported  by  the  abundance  of 
near-seafloor  temperature  and  Eh  anomalies,  and 
the  absence  of  sessile  macrofauna  (Figure  9). 
Oden  and  Loke  are  the  sources  of  numerous  near¬ 
bottom  seawater  temperature  and  Eli  anomalies 
detected  by  the  CAMPER  vehicle,  indicating  wide¬ 
spread  diffuse  hydrothermal  discharge.  Two  plume 
horizons  were  also  detected  over  the  axial  valley 
during  the  2001  AMORE  cruises,  one  at  ~2500  m 
and  one  at  ~3300  m  (AO  >  0.7  C),  and  were  the 
most  extensive  and  voluminous  of  any  on  the  Gakkel 
Ridge  [ Edmonds  et  al.,  2003;  Baker  et  ah,  2004] 
Water  column  plume  experiments  during  the 
AGAVE  cruise  also  identified  hydrothermal  plumes 
confined  within  the  axial  valley  walls  at  3000-3800  m 
that  were  generated  by  vents  on  the  valley  floor. 
Given  their  low-temperature  (T  <  10°C),  these 
hydrothermal  fluids  did  not  disperse  great  distances 
despite  the  weakly  stratified  water  column  and  were 
traced  to  six  spatially  distinct  sites,  including 
plumes  above  Oden  and  Loke  [Strcuwe  et  al,  ,2010]. 
Near-seafloor  temperature  and  Eh  anomalies  iden¬ 
tified  distinct  locations  of  diffuse,  low-temperature 
hydrothermal  fluids  over  100’s  of  meters  on  Oden 
and  the  apron  of  Loke.  The  volcanoes  also  host 
extensive  microbial  mats  that  are  supported  by 
hydrothermal  fluid  (low  [Helmke  et  al.,  2007].  The 
microbial  mats  observed  at  the  85"E  site  are  similar 
in  appearance  (i.e„  orange,  gelatinous  texture)  to 
those  produced  by  Fe-oxidizing  bacteria  found  at 
Loihi  Seamount,  Hawaii  [e.g.,  Emerson  and  Moyer, 
2002;  Edwards  el  al.  201 1  ].  The  presence  of  diffuse 
hydrothermal  discharge  along  both  axis-parallel 
volcanic  lineations  suggests  that  regional  tectonics 
make  the  seafloor  permeable  to  hydrothermal 
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circulation  to  a  hot  rock  body  underlying  the  valley 
floor  and  arc  key  in  controlling  fluid  discharge  sites. 

I  lowevcr.  the  greater  abundance  of  hydrothermal 
anomalies  on  the  southern  lineament  suggests  that 
the  vigor  of  discharge  is  related  to  the  relative  age  of 
volcanic  activity. 

5.2.  Volcanic  Processes  at  85 r  E 

|  jo  |  Our  evaluation  of  seafloor  imagery  from  the 
2007  AClAVb  cruise,  which  is  located  within  the  area 
of  high  intensity  backscatler  identified  by  Edwards 
el  at.  [2001]  (Figure  3).  indicates  that  the  region  has 
experienced  multiple,  small-volume  eruptions  from 
discrete  source  vents  over  time.  Our  results  agree  with 
Edw  ards  et  al.  [2001]  interpretation  of  the  area  of 
high  intensity  backscalter  as  a  region  formed  by 
multiple  eruptions,  which  they  based  on  observations 
of  faults  crosscutting  the  area.  The  sedimentation  rate 
in  the  Gakkel  Ridge  is  estimated  to  be  1-3  cm  kyr  1 
[Bacbiian  et  al..  2004]  and  penetration  into  pelagic 
sediment  by  the  12  kllz  SCICHX  side  scan  sonar 
could  easily  exceed  1  m  without  dramatic  signal 
attenuation  [Mitchell.  1993].  thus  the  region  of  high 
backscalter  at  85E  could  be  produced  by  erup¬ 
tions  >1  million  years  in  age.  All  glasses  sampled 
on  Oden  and  Lokc  [cf.  Shaw  et  al..  2010]  and  vol- 
caniclastic  samples  (Figure  la)  are  age-constrained 
by  U-scrics  disequilibria  (i.e„  “'’Ra/’^'Th  -  3.07  to 
3.65  ±  3%)  to  less  than  8000  years,  and  one  volca- 
niclastic  sample  (Sample  13-1)  is  age-constrained 
(i.e.,  210Pb/22bRa  —  0.89  dt  9%)  to  less  than  100  years 
old  [Elkins.  2009],  Thus,  we  conclude  that  the  large 
region  of  high  intensity  backscalter  from  SC1CEX 
1999  represents  a  zone  of  seafloor  (up  to  I  Myr  in 
age)  with  thin  pelagic  sediment  coverage  and  only  a 
small  area  of  fresh  lavas  on  the  southern  lineament 
may  be  as  recent  as  the  1999  seismic  swarms.  Our 
new  observations  of  small-volume  eruptions  at  the 
Gakkel  Ridge  challenge  the  hypothesis  that  eruption 
volume  varies  inversely  with  spreading  rate  across  the 
global  MOR  system  [Pet fit  and  Chadwick.  1998; 
Sintun  etui.  2002], 

5.3.  Evidence  of  Deep  Explosive  Eruptions 

|ji|  Unconsolidated  volcanictaslic  deposits — 
dominated  by  fine-grained,  glassy  clasts  but  also 
containing  bubble-wall  fragments — blanket  the 
axial  valley  with  variable  thicknesses  up  to  10  cm 
over  an  area  ~5  km  x  10  km.  Although  many 
formation  mechanisms  are  postulated  for  submarine 
basaltic  volcaniclasls  (e.g.,  fuel-coolant  interaction, 
quench  granulation,  hydromagmalic  processes)  the 
widespread  distribution  of  deposits  within  the  axial 


valley  at  depths  greater  than  3.6  km  (e.g.,  at  pres¬ 
sures  greater  than  the  critical  point  of  seawater) 
point  to  formation  by  pyroclastic  explosion  [Biittner 
et  al..  2002;  Zimanowski  and  Biittner.  2003].  In 
particular,  the  abundance  of  the  fresh  volcaniclastic 
deposits  on  lava  surfaces  of  varying  ages  argues 
against  quench  granulation  from  intact  lava  flows  as 
a  primary  mechanism  of  formation.  In  addition,  the 
presence  of  fresh  volcaniclasts  in  the  shallowest 
portions  of  the  study  area,  where  they  are  locally 
piled  on  lop  of  pillows  and  on  top  of  ridgelines, 
suggests  that  the  volcaniclasts  arc  deposited  by 
falling  through  the  water  column  rather  than  flow 
along  the  seafloor.  This  suggests  that  volcaniclasts 
arc  lofted  into  the  water  column  at  the  eruptive  vent. 
Similar  clasts  were  recovered  in  box  cores  at  the 
60' E  volcanic  center  on  the  Gakkel  Ridge  and  may 
indicate  similar  explosive  eruptions  along  the  EVZ, 
though  in  that  case  clasts  were  inteiprcted  as  spall¬ 
ation  from  lava  flow  surfaces  [Miihe  et  al..  1997], 

1 4z |  Volcaniclastic  deposit  thicknesses  correlate  with 
relative  apparent  ages  of  lava  surfaces,  with  the 
thinnest  deposits  on  near-zero-age  lavas  on  Oden  and 
Lokc  Volcanoes.  On  more  oxidized,  relatively  older 
flows  of  Oden  and  Lokc  we  observe  thicker  cover¬ 
age.  Abrupt  changes  in  volcaniclastic  deposit  thick¬ 
ness  at  eruptive  unit  boundaries  indicate  thickness 
variations  related  to  multiple  episodes  of  explosive 
eruption  rather  than  distance  from  a  source,  which 
would  result  in  more  gradual  deposit  variation. 
Deposits  are  thickest  on  the  axial  ridges,  especially 
Duque’s  Hill,  which  we  interpret  to  be  the  oldest 
surfaces  in  the  study  area.  Overall,  there  is  no  clear 
gradient  in  volcaniclastic  deposit  thickness  as  would 
be  expected  with  increasing  proximity  to  an  explo¬ 
sive  source.  If  a  single  explosive  source  were  to  be 
invoked,  it  would  have  to  be  within  or  north  of  the 
two  AVRs.  However,  the  volcaniclastic  deposits  all 
appear  to  be  much  fresher  than  the  lava  flows  on  the 
AVRs,  and  the  source  of  the  youngest  eruptions  in 
the  study  arc  south  of  the  AVRs  on  the  southern 
lineament.  From  settling  velocity  experiments  and 
modeling,  Barreyre  et  al.  [201 1  ]  found  that  a  buoy¬ 
ant  plume  of  hydrothermal  fluid  during  an  eruption 
(i.e.,  hydrothermal  mcgaplumc)  could  disperse  clasts 
~400  m  al  most  from  a  source  vent  on  the  Gakkel 
Ridge  (assuming  nominal  current  speeds  of  ~1  cm/s 
[Stranne  et  al.,  2010]).  Thus,  multiple  source  vents 
distributed  within  the  axial  valley  are  required  to 
disperse  clasts  along  -- 1 0  km  of  seafloor  and  to 
create  abrupt  boundaries  in  deposit  thickness.  We 
further  note  that  the  cratered  morphology  and  pres¬ 
ence  of  unweathered  talus  on  the  summit  and  flanks 
of  Oden  Volcano  arc  consistent  with  a  Vulcanian 
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explosion  mechanism,  where  the  catastrophic  expan¬ 
sion  of  magmatic  volatiles  occurs  in  the  shallow 
subsurface,  thereby  generating  a  central,  talus-filled 
crater  and  ejecting  talus  blocks  onto  the  summit  rim 
and  flanks  of  the  cone. 

5.4.  A  Model  of  Effusive  and  Explosive 
Eruption  at  85°  E 

[«|  The  depth  of  the  axial  valley  of  our  study  area 
at  3600-4200  mbsl  is  deeper  than  the  maximum 
depth  at  which  steam  explosion  can  occur.  The 
physical  characteristics  and  widespread  dispersal  of 
the  pyroclastic  deposits  lead  us  to  suspect  a  highly 
energetic  eruption  driven  by  magmatic  volatiles. 
Given  the  typical  MORB  geochemistry  [Shaw  et  al., 
2010]  and  low  vesicularily  of  the  samples,  we  pre¬ 
sume  that  the  C02  driving  the  eruption  must  have 
accumulated  at  depth,  largely  decoupled  from  the 
erupting  lava  [Head  and  Wilson,  2003].  Thus  mag¬ 
matic  volatiles,  particularly  C02  because  it  has  the 
smallest  solubility  of  the  major  gases,  represent  the 
most  likely  driving  force  for  explosivity  [Bottinga 
and  J avow  1990;  Head  and  Wilson,  2003;  White 
et  al,  2003], 

[ 44 1  C02  exsolution  can  begin  somewhere 
between  45-20  km  below  sea  level  in  upward 
moving  MORB  depending  on  the  initial  concen¬ 
tration  in  the  melt  [Bottinga  and  Javoy,  1990], 
Although  Shaw  et  al.  [2010]  suggest  that  initial 
concentrations  may  noL  have  been  much  higher  than 
the  maximum  C02  that  they  measured  ( 1 596  ppm), 
even  small  supersaturations  at  depths  >1 0  km  would 
allow  for  some  free  gas  to  cxsolvc.  We  propose  that 
over  time  a  significant  volume  of  cxsolved  volatiles 
may  accumulate  in  crustal  magma  chambers.  We 
assess  the  feasibility  of  such  a  model  by  examining 
the  timescales  of  melting  and  magnitude  of  melt 
accumulation  necessary  to  generate  the  required 
amounts  ofC02  (see  Appendix  A).  By  modeling  the 
amount  of  magma  produced  due  to  passive  mantle 
upwcliing  [e.g..  Forsyth,  1993]  and  calculating  the 
supcrsaturalion  in  CO,  at  a  given  storage  depth 
within  the  crust,  we  estimate  that  realistic  timescales 
for  gas  accumulation  would  be  on  the  order  of 
1500-2000  years.  As  the  expected  repose  interval 
for  eruptions  at  ultraslow  spreading  ridges  is 
1 0.000  to  1 00.000  years  [e.g.,  Perfu  and  Chadwick. 
1 998],  this  timescale  for  melt  accumulation  suggests 
storage  of  COi  in  crustal  reservoirs  is  a  viable 
mechanism. 

1 45]  Gas  bubbles  can  form  a  foam  layer  in  a  reser¬ 
voir  with  a  gas  volume  fraction  that  can  be  stable 
over  90%  [Vergniolle  and  Jaupart ,  1990].  The 


sequestered  volatiles  may  then  be  released  during 
an  eruptive  episode,  leading  to  explosive  discharge 
spreading  pyroclastic  material  over  seafloor  areas 
proximal  to  the  source  vent  [Head  and  Wilson, 
2003;  Sohn  et  al,  2008].  We  propose  that  gas- 
deplclcd  magma,  rising  more  slowly  in  an  eruption 
conduit,  may  also  erupt  during  these  explosive 
episodes.  A  slowly  rising  magma  would  erupt  low- 
vesicularity,  effusive  flows  in  equilibrium  with 
C02  at  seafloor  depths,  similar  to  those  observed  al 
the  85' t  segment  [e.g.,  Shaw  et  al,  2010],  This 
model  is  consistent  with  our  observations  of  wide¬ 
spread  pyroclastic  deposits  and  avcsicular.  effusive 
lavas,  and  this  type  of  eruption  has  been  directly 
observed  at  W.  Mata  Volcano  in  the  Lau  Basin 
[Re sing  el  al,  2011]. 

6.  Conclusions 

(•«.]  Recent  volcanic  activity  at  the  85  E  segment  of 
the  ultraslow  spreading  Gakkel  Ridge  consists  of 
small  volume  effusive  and  explosive  eruptions 
sourced  from  multiple  volcanic  features  within  the 
axial  valley.  These  eruptions  have  produced  AVRs 
and  volcanic  cones,  with  many  of  the  latter  exhi¬ 
biting  distinctive  central  craters.  Similar  volcanic 
cones  and  AVRs  arc  observed  in  the  axial  valley  at 
the  very  slow  spreading  Southwest  Indian  Ridge 
[Mendel  and  Sauter,  1997;  Cannot  et  al,  1999; 
Sauter  et  al,  2001]  as  well  as  the  slow  spreading 
Mid-Atlantic  Ridge  [Smith  and  Conn,  1990],  All  of 
the  “fresh”  lavas  observed  in  our  surveys  were 
found  on  the  southernmost  volcanic  lineament  al  the 
Oden,  Thor,  and  Loke  volcanic  cones,  which  also 
coincides  with  the  region  where  near-bottom  Eh  and 
thermal  anomalies  were  observed  in  the  water  col¬ 
umn.  and  where  microbial  mats  were  observed  on 
the  seafloor.  These  observations  demonstrate  that 
the  1999  seismic  swarm  may  have  been  associated 
with  a  series  of  small  volume  eruptions  featuring 
both  explosive  and  effusive  discharge. 

[47]  We  model  the  widespread  volcaniclastic  depos¬ 
its  observed  on  the  seafloor  at  the  85“  E  segment  as 
having  been  generated  by  the  explosive  discharge  of 
C02  that  accumulated  in  (possibly  deep)  crustal 
melt  reservoirs  over  a  period  of  time  on  the  order 
of  the  eruption  interval  for  an  ultraslow-sprcading 
ndge  (~  10,000  vr).  We  envision  that  COi  ‘bubbles’ 
escaped  from  crustal  melt  reservoirs  during  eruptive 
events,  rising  through  a  more  slowly  ascending 
magma  column  before  explosively  discharging  at.  or 
in  some  cases  just  beneath,  the  seafloor.  The  energy 
released  during  explosive  discharge,  combined  with 
the  buoyant  rise  of  hot  fluid,  lofts  fragmented  clasts  of 
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rapidly  cooling  magma  into  the  water  column,  and 
they  subsequently  settle  onto  the  seafloor  as  fall 
deposits  in  the  region  surrounding  the  source  vent. 
The  volcaniclasis  themselves  arc  avcsicular  with  dis¬ 
solved  CO?  in  equilibrium  with  seafloor  piessures 
because  magmas  (unlike  the  gas  bubbles)  rise  slowly 
in  the  eruption  conduit,  and  therefore  continually 
cxsolve  excess  CO?  along  their  vertical  path  to  the 
seafloor  from  the  crustal  melt  reservoir,  hxcess  gas 
exsolved  Irom  magmas  during  a  slow  rise  to  the  sea¬ 
floor  may  be  incorporated  into  gas  ‘bubbles'  escaping 
from  the  melt  reservoir,  and  thereby  enhance  explo¬ 
sive  activity. 

Appendix  A:  Accumulation  of  C02 
in  a  Crustal  Magma  Chamber 

1 48 1  The  physical  characteristics  of  the  volcaniclas- 
tic  deposits  lead  us  to  suspect  a  highly  energetic 
eruption  driven  by  magmatic  volatiles.  Given  the 
typical  MORB  geochemistry  and  low  vesicularity  of 
(he  samples,  we  presume  that  the  CO?  driving  the 
eruption  must  have  accumulated  at  depth,  largely 
decoupled  from  the  erupting  lava.  We  assess  the 
feasibility  of  such  a  model  by  examining  the  time¬ 
scales  of  melting  and  magnitude  of  melt  accumula¬ 
tion  necessary  to  generate  the  required  amounts  of 
COj.  Although  the  volume  of  gas  required  to  drive 
an  explosive  eruption  in  this  environment  is  not 
known,  a  conservative  estimate  is  that  the  gas  vol¬ 
ume  is  3  times  that  of  the  magma  volume  emplaced 
at  the  seafloor  [e.g.,  Sohn  et  ai,  2008;  Head  and 
Wilson ,  2003]. 

|4«|  A  minimum  time  to  accumulate  the  CO?  is 
determined  by  assuming  that  CO?  behaves  per¬ 
fectly  incompatibly  and  that  all  CO?  extracted  from 
the  mantle  is  available  to  contribute  to  an  explosive 
eruption.  The  latter  assumption  is  likely  unreason¬ 
able,  but  places  a  lower  bound  on  the  timescale  of 
melt  generation.  If  all  mantle  that  passes  through 
the  melting  zone  is  stripped  of  its  carbon  by  newly 
formed  magmas,  we  can  estimate  the  CO?  accu¬ 
mulation  by  assuming  a  mantle  carbon  concentra¬ 
tion  and  the  length  of  ridge  over  which  melts  are 
focused.  The  volume  flux  of  mantle  (£?,„)  through 
the  melting  region  is  related  to  spreading  rate  by. 

Q.  -  (L,L,)(2U0/ir)  (At) 

in  the  case  of  a  uniform  viscosity  mantle  [Reid  and 
Jackson.  1981].  The  melting  column  width  (Lx)  is 
derived  from  the  height  of  the  melting  column, 
which  we  set  to  40  km  (Figure  S10  in  the  auxiliary 
material).  The  along-axis  length  of  ridge  from 


which  melts  accumulate  (/., )  is  set  to  30  km,  con¬ 
sistent  with  the  spacing  of  volcanic  centers  in  the 
tastern  Volcanic  Zone.  The  full  spreading  rate  ( U„) 
is  9  mm  yr  ‘.  With  a  mantle  density  of 
3500  kg  m  \  a  mantle  CO?  concentration  of 
134  ppnt  [Shaw  et  ai..  2010],  and  assuming  ideal 
gas  behavior,  the  time  necessary  to  extract  suffi¬ 
cient  C02  front  the  mantle  is  1130  years  for  a 
106  nv  eruption  (lime  scales  linearly  with  eruption 
volume  in  this  calculation). 

[so]  A  ntoie  realistic  estimate  of  the  timescale  of  melt 
accumulation  can  be  detemtined  by  modeling  the 
amount  of  magma  produced  due  to  passive  mantle 
upwclling  and  calculating  the  supersaturation  in  CO? 
at  a  given  storage  depth  within  the  crust.  Here  we  rely 
on  the  formulation  from  Forsyth  [1993],  where  melt 
production  per  unit  length  of  ridge  per  unit  time  (M ) 
is  related  to  spreading  rate  by; 


where  B  is  the  mantle  melt  productivity  (fraction  of 
melt  per  km  of  uplift  above  the  solidus)  and  Za  is  the 
thickness  of  the  melting  column  (depth  of  solidus 
measured  from  the  base  or  the  crust).  The  amount  of 
melt  produced  over  the  region  (Figure  S10  in  the 
auxiliary  material)  as  a  function  of  time  is  combined 
with  the  CO?  contributed  Horn  the  mantle  over  the 
same  time  interval  and  then  the  magma  is  stored  at  a 
depth  of  7.5  km  depth  in  the  crust,  consistent  with 
data  from  melt  inclusions  from  this  area  [5/iair  etal.. 
2010].  At  that  depth  CO?  solubility  is  ~  1 200  ppm 
(0.0012  wt%)  given  observed  water  concentrations 
of  0.25  wl%  {Shaw  et  ul.,  2010).  The  annual  melt 
production  wx  calculate  of  1.87  x  1  if  kg  yr  1  and 
the  annual  CO?  contribution  from  the  mantle  of 
6.1 8  x  10h  kg  yr  1  results  in  CO?  concentrations  in 
the  melt  of  ~3300  ppm.  Thus,  nearly  two-thirds  of 
the  CO?  extracted  from  the  mantle  could  be 
stored  as  a  free  gas  in  a  reasonably  placed  crustal 
magma  chamber.  As  such,  more  realistic  time- 
scales  for  gas  accumulation  w'ould  be  on  the  order  of 
1500-2000  years.  As  the  expected  repose  interval 
for  eruptions  at  ullraslow-spreading  ridges  is 
10,000  to  100,000  years  [e.g..  Pet  fit  and  Chadwick , 
1998],  this  timescale  for  melt  accumulation  sug¬ 
gests  storage  of  CO?  in  crustal  reservoirs  is  a  v  iable 
mechanism. 
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Chapter  3:  Reaction-driven  microearthquakes  at  the  Trans-Atlantic 
Geotraverse  (TAG)  active  mound,  Mid-Atlantic  Ridge,  26°08’N 


Abstract 

A  9-month  deployment  of  five  ocean  bottom  seismometers  (OBS)  on  the  periphery  of  the 
Trans-Atlantic  Geotraverse  (TAG)  active  mound  detected  32,078  very  small,  local 
microearthquakes  (average  ML  =  -0.95)  on  at  least  four  instruments  from  June  2003  to 
February  2004.  Events  occurred  at  a  constant  rate  of  up  to  ~243  events  per  day  for 
several  months  prior  to  an  instrument  failure  and  are  characterized  by  durations  of  <1 
second  and  single-phase  P-wave  arrivals  (i.e.,  no  5-arrivals).  Accurate  locations  were 
determined  for  6,207  of  the  earthquakes,  with  hypocenters  clustered  on  the  south  and 
west  flanks  of  the  deposit  at  depths  of  <  125  m  below  the  seafloor.  The  characteristics  of 
the  microearthquakes  and  their  location  within  the  horizon  of  mineral  deposition 
observed  in  drill  cores  are  consistent  with  a  ‘frost  heave’  fracturing  mechanism  driven  by 
anhydrite  deposition  in  the  shallow  secondary  recharge  circulation  system  of  the  TAG 
active  mound. 
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1.  Introduction 


Seafloor  massive  sulfide  deposits  are  the  modern  equivalents  of  ancient  ophiolite-hosted 
ore  deposits,  such  as  those  mined  for  copper  and  gold  in  Cyprus,  Oman,  and 
Newfoundland  [Humphris  et  al.,  1995;  Hannington  et  al.,  1998].  They  were  first 
recognized  during  explorations  of  deep-sea  hydrothermal  vent  systems,  where 
interactions  between  seawater  and  oceanic  crust  build  massive  sulfide  deposits  on  the 
seafloor  [Scott  et  al.,  1974a,  Scott  et  al.,  1974b;  Rona  et  al.,  1975;  Rona  et  al.,  1986; 
Campbell  et  al.,  1988;  Karson  and  Rona,  1990;  Herzig  et  al.,  1998].  The  TAG 
hydrothermal  field  on  the  seafloor  of  the  Mid- Atlantic  Ridge  at  26°N  (spreading  rates  13 
mm  year  1  to  the  east  and  1 1  mm  year  1  to  the  west)  [McGregor  and  Rona,  1975; 
McGregor  et  al.,  1977]  consists  of  a  zone  of  low-temperature  deposits  and  alteration  and 
three  zones  of  high-temperature  deposits  and  alteration  spanning  ages  up  to  100,000 
years  old  [Rona  et  al.,  1975,  1986,  1993;  Thompson  et  al.,  1985;  Karson  and  Rona, 

1990], 

Genesis  of  seafloor  massive  sulfide  deposit  requires  a  heat  source,  and  sometimes 
includes  a  secondary  recharge  system  of  seawater  entrainment,  a  primary  system  of 
hydrothermal  circulation,  and  chemical  reaction  with  host  rock.  The  TAG  hydrothermal 
field  sits  in  the  extensional  environment  of  a  ~5  km  long  shallow  dipping  (20°)  dome¬ 
shaped  detachment  fault  [Tivey  et  al.,  1993;  Canales  et  al.,  2007;  deMartin  et  al.,  2007] 
(Figure  1A).  The  detachment  fault  hanging  wall  is  characterized  by  intersecting  faults 
and  fissures  thought  to  structurally  control  the  focusing  of  hydrothermal  activity  to  the 
seafloor  [Zonenshain  et  al.,  1989;  Rona,  1990;  Karson  and  Rona,  1990;  Kleinrock  and 
Humphris,  1996].  The  seismogenic  zone  of  the  detachment  fault  extends  to  7-8  km  below 
the  seafloor  and  may  serve  as  a  pathway  for  primary  hydrothermal  fluid  circulation  that 
fuels  current  high  temperature  venting  [deMartin  et  al.,  2007] .  High  temperature  activity 
is  currently  focused  at  the  TAG  active  mound,  one  of  at  least  four  large  hydrothermal 
deposits  in  the  field  [White  et  al.,  1998]. 
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Investigations  of  the  TAG  active  mound  have  revealed  a  zoned  subsurface  stratigraphy 
resulting  from  low  and  high-temperature  reactions  and  an  abundance  of  anhydrite  within 
the  deposit  built  by  complex  interactions  between  hydrothermal  fluid,  entrained  seawater, 
and  stockwork  [Herzig  et  al.,  1998] .  Conceptually,  the  mound  is  built  both  by  a  focused 
primary  fluid  circulation  system  and  by  a  secondary  circulation  system  of  conductively  or 
convectively  heated  seawater,  which  is  entrained  into  the  deposit  [c.f.  Humphris  and 
Tivey,  2000] .  The  mound  has  developed  as  a  breccia  pile  through  periods  of  venting 
activity  and  quiescence  over  at  least  20,000  years  [Rona,  1974;  Lalou  et  al.,  1990],  and 
exhibits  complex  short-term  episodic  and  periodic  behaviors  of  venting  [Sohn,  2007 a,b\. 
In  spite  of  the  large  body  of  knowledge  on  the  TAG  active  mound,  some  key  aspects  of 
the  hydrothermal  system  remain  unclear.  In  particular,  relatively  little  is  know  about  the 
physical  extent  and  geometry  of  hydraulic  pathways  and  subsurface  fluid  flow.  Defining 
the  zone  of  secondary  recharge  is  important  for  understanding  where  mineral  deposition, 
especially  that  of  anhydrite,  occurs  within  the  mound  and  the  rate  at  which  deposition 
occurs. 

Hydrothermal  activity  intrinsically  generates  many  microearthquakes  [e.g.,  Kieffer,  1984; 
Ross  et  al.,  1996;  Bianco  et  al.,  2004;  Foulger  et  al.,  2004,  etc.].  Studies  of  hydrothermal 
systems  have  been  done  on  land,  but  due  to  the  difficulty  of  submarine  observation,  most 
oceanic  microearthquake  experiments  have  been  on  the  scale  of  volcanic  centers  or  ridge 
segments.  Marine  studies  on  fast  spreading  mid-ocean  ridge  systems  have  shown  a  direct 
correlation  between  fluid  flow  and  local  seismicity  [Sohn  et  al.,  1998;  Johnson  et  al., 
2000;  Davis  et  al.,  2004],  however  most  studies  represent  feedbacks  between  tectonic 
faulting  and  venting.  We  know  that  changes  in  the  flow  within  the  subsurface  at  the 
TAG  active  mound  are  manifest  at  its  surface  by  changes  in  its  morphology  and  in  the 
distribution  of  venting  on  very  short  time  scales  [Kleinrock  and  Humphris,  1996;  Sohn, 
2007a,b].  This  is  the  first  time  that  a  seismic  network  small  enough  to  detect  and  locate 
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microearthquakes  intrinsic  to  hydrothermal  processes  has  been  deployed  at  a  seafloor 
massive  sulfide  deposit. 


Analysis  of  small  microearthquakes  intrinsic  to  hydrothermal  processes  in  the  shallow 
crust  requires  a  dense  on-bottom  seismic  array  with  short  inter-element  spacings,  in  this 
case  of  98  m  to  203  m.  To  investigate  hydraulic  phenomenon,  we  deployed  a  dense  small 
aperture  array  of  OBSs  in  a  ring  on  the  periphery  of  the  TAG  active  mound.  We  detected 
32,078  microearthquakes  in  a  9-month  deployment,  and  located  6,207  events  with 
standard  errors  of  20  ms  or  less.  The  microearthquake  hypocenters  collapse  on  the 
southern  and  western  extents  of  the  sulfide  deposit  in  the  top  125  m  of  crust.  We  use  our 
results  to  investigate  natural  hydraulic  fracturing  processes  and  the  maintenance  of  fluid 
pathways  beneath  the  active  hydrothermal  system. 


2.  The  TAG  Active  Mound 

The  TAG  active  mound  is  located  within  the  rift  valley,  approximately  3  km  west  of  the 
detachment  fault  termination  and  2.4  km  east  of  the  ridge  axis  [Tivey  et  al.,  2003] .  Within 
the  hanging  wall  block,  orientations  of  ~40%  of  fissures  deviate  significantly  (  >  45°) 
from  those  predicted  by  the  regional  tectonic  stress  field,  and  small  extensional  fissures 
(several  tens  of  meters  in  length  and  0.7  m  wide  on  average)  crosscut  100,000-year-old 
crust  (Figure  2B)  [ Kleinrock  and  Humphris ,  1996;  Bohnensteihl  and  Kleinrock,  2000]. 
First  order  estimates  indicate  that  these  fissures  are  restricted  to  <  500  m  depth  with  a 
mean  value  of  70  m  [ Bohnenstiehl  and  Kleinrock,  2000] .  These  results  indicate  that  the 
deposit  is  hosted  on  a  highly  permeable  crustal  block  with  a  stress  field  that  has  been 
modified  by  the  detachment  faulting  process. 

The  TAG  active  mound  is  a  ~200  m  diameter,  ~50  m  tall  massive  sulfide  deposit  that 
discharges  hydrothermal  fluids  in  excess  of  360°C  with  a  total  heat  flux  of  ~1  GW  [Rona, 
1984,  Wickers,  2005].  On  the  basis  of  geochronological  studies,  high  temperature  venting 
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at  the  TAG  active  mound  is  thought  to  have  begun  20,000  -  50,000  years  ago,  with 
renewed  pulses  of  activity  every  5,000-6,000  years.  These  studies  suggest  that  the  most 
recent  pulse  of  activity  began  ~90  years  ago  [Lalou  et  al.,  1990,  1993,  1998;  You  and 
Bickle,  1998],  The  active  mound  is  believed  to  experience  periods  of  quiescence  during 
which  anhydrite  dissolves  in  response  to  cooling  temperatures,  leading  to  mass  wasting 
and  extensive  reworking  of  the  sulfides  [ Lalou  et  al.,  1990;  Humphris  and  Cann,  2000]. 

High-temperature  fluid  discharge  is  focused  in  an  ~80  m  wide  upflow  zone  to  the  Black 
Smoker  Complex  (BSC),  which  constitutes  a  set  of  sulfide  chimneys  at  the  top  of  the 
mound  [ Honnorez  et  al.,  1998].  Lower  temperature  fluids,  presumably  resulting  from 
mixing  of  high-temperature  fluids  with  cold  seawater,  discharge  from  a  multitude  of 
locations  on  the  mound  [ Campbell  et  al.,  1988;  Edmond  et  al.,  1995;  Humphris  et  al., 
1995;  Tivey  et  al.,  1995;  Edmonds  et  al.,  1998].  Exit-fluid  temperature  time-series  data 
acquired  from  the  active  mound  demonstrate  that  discharge  and  flow  reversals  in  exit 
fluid  temperature  data,  indicative  of  seawater  entrainment  into  the  mound  surface,  occur 
on  complex  episodic  and  periodic  time  scales  [Sohn,  2007a, b]. 


3.  Experiment 

An  ocean  bottom  seismometer  (OBS)  network  was  deployed  at  the  TAG  hydrothermal 
field  from  June,  2003  to  February,  2004  as  part  of  the  Seismicity  and  Fluid  Flow  of  the 
TAG  Hydrothermal  Mound  (STAG)  Experiment  (Figure  2A).  The  OBS  network 
consisted  of  three  rings  of  instruments  centered  on  the  active  hydrothermal  mound  with 
apertures  of  200  m,  2  km,  and  8  km.  Here  we  report  results  from  the  smallest  aperture 
sub-network,  which  was  deployed  around  the  periphery  of  the  TAG  active  mound  in  an 
attempt  to  record  very  small  microearthquakes  associated  with  fluid  flow  in  the  deposit 
(Figure  2A).  The  sub-network  consisted  of  five  WHOI  D2  OBSs,  which  are  four- 
component  (3  geophone  +  1  hydrophone  channel),  short-period  (4.5  Hz  natural 
frequency)  instruments.  The  seismometers  recorded  continuous  data  at  100  Hz  for  the 
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nine-month  duration  of  the  deployment,  except  for  STIR-12,  which  failed  77  days  into 
the  study  for  unknown  reasons.  The  instruments  were  deployed  by  wire  (rather  than  free- 
fall)  from  the  R/V  Knorr,  and  were  positioned  around  the  perimeter  of  the  active  mound 
using  long  baseline  navigation  techniques.  This  was  the  first  deployment  of  a  small- 
aperture  (inter-element  spacings  of  98  -  203  m)  OBS  network  on  a  deep-sea  hydrothermal 
field. 


4.  Data  Analysis 

We  present  analyses  for  microearthquakes  that  were  detected  only  on  the  200  m  aperture 
sub-network  in  order  to  focus  on  shallow  seismic  activity  beneath  the  hydrothermal 
mound.  The  OBS  network  also  detected  a  large  number  of  local  and  regional  tectonic 
events  [de Martin  et  al.,  2007],  but  these  events  are  excluded  from  our  analyses.  In 
addition,  we  analyze  events  that  are  only  detected  by  at  least  four  instruments  on  the 
small  aperture  sub-network.  We  picked  arrivals  for  32,078  microearthquakes  using  a 
short-term  average  to  long-term  average  algorithm,  and  then  refined  these  picks  using 
spectral  methods  [e.g.  Song  et  al.,  2010] .  Linear  clock  drift  corrections  were  applied  to 
each  station,  and  these  corrections  were  further  refined  to  account  for  non-linear  clock 
drift  by  cross-correlating  long-wavelength  (>  500  m)  microseism  signals  [Webb,  1992; 
Webb  and  Crawford,  1999]  propagating  across  the  small-aperture  sub-network.  Finally, 
station  corrections  were  applied  to  each  record  by  subtracting  the  mean  arrival  time 
residual  for  all  events. 

Hypocenters  were  estimated  using  a  stochastic  grid  search  routine  [Sohn  et  al.,  1998] 
minimizing  root  mean  square  (rms)  travel-time  residuals  between  the  phase  picks  and 
raytracing  predictions  over  a  1 .2  km  cubic  grid  with  40  m  node  spacing.  We  used  the 
average  seismic  (P-wave)  velocity  vs  depth  profile  for  the  TAG  region  from  Canales  et 
al.,  [2007]  as  a  starting  point  for  developing  a  velocity  model  for  the  raytracer.  Because 
the  seismic  velocities  in  the  shallowmost  crust  (i.e.,  extrusive  lithologies  in  seismic  layer 
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2A)  are  not  well-constrained  in  the  Canales  et  al.,  [2007]  model  we  systematically 
perturbed  the  shallow  velocities  and  selected  the  model  that  minimized  hypocentral  errors 
for  a  subset  of  1000  randomly  selected  events.  Our  final  model  has  a  P-wave  velocity  of 
2.5  km  s'1  at  the  seafloor  that  increases  linearly  with  depth  to  a  velocity  of  ~4  km  s  '  at  1 
km  below  the  seafloor,  and  it  is  draped  over  the  local  bathymetry,  including  the  active 
mound  [White  et  al.,  1998;  Roman  and  Singh,  2005].  The  velocity  model  does  not 
account  for  the  different  velocity  structure  of  the  sulfide  deposit  in  the  oceanic  crust, 
which  may  add  error  to  hypocenter  estimates.  Microearthquake  moments  and  local 
magnitudes  were  calculated  using  standard  techniques  [e.g.,  Deichmann,  2006]. 


5.  Results 

We  located  32,078  events  detected  during  the  9-month  OBS  deployment.  Seismicity  rates 
during  the  deployment  were  essentially  constant,  although  the  rate  decreased  from  243  to 
128  microearthquakes  per  day  when  STIR- 12  failed  77  days  into  the  experiment,  and 
then  decreased  again  to  97  events  per  day  when  whale  calls  [e.g.,  McDonald  et  al.,  1995] 
in  the  study  area  masked  signals  from  the  smallest  earthquakes  (Figure  3).  The  space- 
time  microearthquake  patterns  do  not  exhibit  mainshock-aftershock  sequences,  swarm¬ 
like  behavior,  nor  evidence  for  migration.  The  spectrum  of  the  microearthquake  catalog 
(formed  by  treating  the  catalog  as  a  point  process,  e.g.,  Brillinger,  [1974])  does  not 
exhibit  strong  periodicities  and  is  not  coherent  with  concurrently  measured  deep-sea  tide 
gauge  data  (i.e.,  ocean  tides)  nor  predicted  solid  earth  tides  (Figure  4). 

Seismic  moments  for  the  microearthquakes,  estimated  from  average  event  spectral  levels 
between  5  -  50  Hz  [ Pearson ,  1982],  fall  between  10l3-1016  dyne-cm,  corresponding  to 
local  magnitudes  of -1.4  <  ML  <0.5,  with  an  average  of  -0.95  (Figure  5).  The  total 
moment  release  for  these  events  over  the  9-month  deployment  is  1018  6  dyne-cm,  which  is 
equivalent  to  a  single  event  with  a  moment  of  ~  1.7.  For  comparison,  the 
microearthquakes  associated  with  detachment  faulting  and  tectonic  extension  observed  at 
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the  TAG  segment  over  the  same  time  period  with  the  large  aperture  OBS  network  had 
local  magnitudes  ranging  from  1  <  ML  <4  [deMartin  et  al.,  2007]. 


5.1  Microearthquake  event  characteristics 

Microearthquakes  are  characterized  by  impulsive,  short-period  (15-30  Hz),  P-wave 
arrivals  with  durations  of  <1  s  (Figure  6).  Event  spectra  contain  peaks  related  to  the 
frequency  content  of  the  source  mechanism,  as  well  as  peaks  at  individual  stations 
corresponding  to  coupling  frequencies.  5-arrivals  are  not  generally  observed  in  the  data 
records,  indicating  that  the  events  are  generated  by  volume  changes  rather  than  shear 
failure  on  a  fault  surface.  The  small  size  of  the  earthquakes  and  poor  signal-to-noise  ratio 
makes  it  difficult  to  accurately  determine  the  polarity  of  first  arrivals.  This,  combined 
with  the  small  number  of  stations  available  for  analysis,  precludes  the  estimation  of 
reliable  focal  plane  solutions  for  the  microearthquakes.  Furthermore,  many  of  the  events 
are  demonstrably  within  the  seismic  near-field  (i.e.,  less  than  one  wavelength  from  the 
stations),  which  complicates  waveform  analysis  and  interpretation.  We  cross-correlated 
waveforms  for  5000  randomly  selected  events  with  the  next  100  events  in  time  and  100 
additional  random  events  from  the  catalog  to  investigate  waveform  similarities  that  could 
be  exploited  for  analytical  purposes.  We  found  that  less  than  1%  of  the  catalog  is 
correlated  above  a  threshold  of  r  >  0.7,  such  that  it  was  not  possible  to  use  relative 
relocation  nor  stacking  methods  to  refine  our  hypocentral  and  arrival  polarity  estimates. 
These  results  argue  against  repeating  sources,  although  low  signal-to-noise  ratios  and 
near- field  propagation  effects  may  also  contribute  to  the  lack  of  waveform  similarity. 


5.2  Hypocenter  locations 

We  estimated  hypocenters  for  32,078  events,  but  roughly  half  of  these  estimates  have  rms 
errors  >200  ms,  and  are  therefore  not  suitable  for  detailed  analyses  and  interpretation. 
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This  is  to  be  expected  considering  the  small  number  of  stations  (4-5  OBSs)  available  for 
analysis  and  the  lack  of  5-arrivals.  Most  of  the  hypocenters,  which  have  four  degrees  of 
freedom  (t0,  x,  y,  z),  were  estimated  with  four  P-wave  arrival  times,  leading  to  large 
uncertainties,  particularly  for  events  well  outside  the  network  aperture.  We  addressed  this 
issue  by  focusing  our  analyses  on  the  6,207  events  with  hypocenter  rms  errors  <  20  ms 
(Figure  7).  Hypocenters  for  these  events  (Figure  8)  have  average  l-o  confidence  limits  of 
214  m,  202  m,  and  187  m  in  the  north-south,  east-west,  and  vertical  directions, 
respectively.  The  most  intense  activity  occurs  to  the  south  and  west  of  the  active  mound 
within  a  narrow  depth  interval  down  to  ~200  mbsf.  The  velocity  model  used  in 
microearthquake  location  is  a  simple  one-dimensional  velocity  model  draped  over  a 
three-dimensional  grid.  Microearthquake  hypocenters  often  fall  outside  of  the  network. 

It  is  possible  that  many  waveforms  travel  paths  through  the  sulfide  deposit,  which  is 
likely  to  have  a  very  different  velocity  structure  than  the  assumed  model.  If  the  different 
crustal  structure  beneath  the  sulfide  deposit,  which  is  a  breccia,  has  lower  seismic 
velocities,  calculated  locations  of  epicenters  would  locate  closer  to  the  hydrothermal 
mound  than  they  are  shown  in  Figure  8 . 

The  low  number  of  stations  available  for  analysis  combined  with  the  low  signal-to-noise 
ratio  for  many  events  produces  considerable  scatter  in  the  hypocenters.  Relative 
relocation  methods  based  on  waveform  cross-correlation  cannot  be  used  to  reduce  scatter 
because  of  the  lack  of  waveform  similarity,  but  the  collapsing  method  [Jones  and 
Stewart,  1997]  offers  an  appealing  alternative.  The  collapsing  method  seeks  the  simplest 
structure  that  can  be  fit  to  the  hypocenters  by  shifting  event  locations  towards  the 
geometric  locus  of  the  3-dimensional  error  volumes.  In  other  words,  event  hypocenters 
are  moved  within  a  confidence  interval  in  a  direction  toward  the  center  of  mass  of  the 
events  within  the  confidence  interval  (i.e.  error  ellipsoid).  The  only  variable  in  the 
method  is  the  significance  level  of  the  confidence  ellipsoid.  Collapsing  method  does  not 
suggest  that  closely-collapsed  events  have  similar  source  mechanisms  or  wave  paths. 
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Before  applying  the  method  we  divided  the  events  into  two  separate  catalogues 
corresponding  to  before,  and  after,  the  failure  of  STIR- 12  so  that  the  results  are  not 
influenced  by  the  change  in  network  geometry  that  occurred  at  that  time.  After  applying 
the  collapsing  method,  the  epicenters  delineate  a  narrow  (~  50  m  wide)  zone  of  seismicity 
that  wraps  around  the  southwest  comer  of  the  active  mound  (Figure  9).  Many  of  the 
events  cluster  in  a  narrow  depth  interval  at  ~  125  mbsf,  with  essentially  no  activity  above 
this  horizon,  and  more  diffuse  activity  below.  Some  north-south  striping  is  evident  in 
hypocenters  to  the  south  of  the  mound,  but  this  is  most  likely  a  result  of  the  north-south 
orientation  of  the  error  ellipses  in  this  area  rather  than  geologic  structure.  Comparing  the 
results  of  the  collapsing  method  (Figure  9)  to  a  density  map  of  uncollapsed  epicenters 
(Figure  10),  we  see  that  the  method  simply  serves  to  highlight  the  densest  seismic  zone. 


6.  Discussion 

Deploying  a  small- aperture  (~200  m)  network  of  five  OBSs  around  the  perimeter  of  the 
TAG  active  mound  allowed  us  to  detect  and  locate  large  numbers  (~32,000)  of  very  small 
(Ml  ~  -1)  microearthquakes  within  the  zone  of  mineral  deposition  and/or  alteration  in  the 
shallow  crust  (depth  <  200  m).  These  events  are  characterized  by  short-duration  (<  Is), 
impulsive  P- waves  with  little  to  no  evidence  for  accompanying  5-arrivals,  indicating  a 
non-double-couple  (i.e.,  volume  change)  as  opposed  to  a  shear  failure  source  mechanism. 
Flow-induced  microearthquakes  generated  by  volume  change  mechanisms  have  been 
observed  in  geothermal  areas  [e.g.,  Ross  et  al.,  1996;  Miller  et  al.,  1998;  Foulger  et  al., 
2004]  and  hydraulic  fracturing  experiments  [e.g.,  Aki  et  al.,  1982;  Cuenot  et  al.,  2006; 
Julian  et  al.,  2007;  Stteny  et  al.,  2009],  but  this  is  the  first  time  that  a  network  small 
enough  to  detect  and  locate  them  has  been  deployed  at  a  deep-sea  hydrothermal  field. 

The  dissimilarity  between  seismograms  from  different  stations  for  the  same  event,  despite 
the  small  network  aperture,  is  striking.  This  is  partly  attributable  to  coupling  effects  -  for 
example  the  spectral  plots  in  Figure  6  show  that  STIR- 13  has  a  coupling  resonance  at  ~15 
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Hz  -  but  most  of  the  waveform  dissimilarity  is  likely  caused  by  scattering.  Dissimilar 
microearthquake  waveforms  have  also  been  observed  on  small- aperture  networks 
deployed  at  the  Old  Faithful  geyser  in  Yellowstone  National  Park  [e.g.,  Kiejfer,  1984; 
Kedar  and  Kanamori,  1998],  suggesting  that  hydrothermal  systems  are  associated  with  a 
high  degree  of  structural  heterogeneity  in  the  shallow  crust.  This  is  consistent  with  the 
deposition  of  hydrothermal  minerals  in  a  complex  fracture  network,  and  mound  building, 
in  general,  and  it  suggests  that  a  sufficiently  dense  seismic  network  may  have  the 
potential  to  image  these  subsurface  features. 

Certain  spectral  peaks  are  observed  at  all  stations  for  individual  events,  and  these  are 
likely  associated  with  source  processes.  For  the  events  shown  in  Figure  6,  for  example, 
there  are  common  spectral  peaks  in  the  20  -  40  Hz  range,  and  if  the  source  mechanisms 
were  known  it  might  be  possible  to  use  these  peaks  to  constrain  the  physical  properties 
and  length  scales  of  the  source  region.  Unfortunately  we  do  not  have  enough  information 
to  conduct  this  type  of  analysis  with  the  data  from  this  experiment,  but  here  again  we  can 
see  the  potential  value  of  deploying  a  dense  seismic  network  on  a  deep-sea  hydrothermal 
field. 

Most  of  the  events  we  located  are  within,  or  near  the  limit  of,  the  seismic  near-field  of  the 
network.  Seismic  wavelengths  for  the  P-arrivals  are  on  the  order  of  50-150  m,  which  is 
comparable  to  the  hypocentral  distance  for  many  events  to  some  seismometers,  such  that 
the  far-field  assumption  implicit  in  our  analyses  may  not  always  be  valid.  Near-field 
propagation  intertwines  with  both  P  and  S  arrivals,  thus  contaminating  their  observed 
polarizations  [Aki  and  Richards,  2002;  Lokmer  and  Bean,  2010],  and  thus  could 
potentially  play  a  role  in  the  waveform  dissimilarity  that  we  observed.  The  near-field 
term  decays  as  1/r2  to  a  distance  of  half  the  P-wavelength  (~25-75  m)  from  the  source  and 
1/r3  beyond  that  with  distance  (r)  from  the  source  [ Lokmer  and  Bean,  2010] .  Thus  near¬ 
field  disturbance  cannot  explain  the  lack  of  S-arrivals  in  our  data,  because  most  event 
hypocenters  lie  outside  of  the  network  and  are  more  than  a  wavelength  from  at  least  some 
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stations  in  the  network.  Filtering  may  be  used  to  suppress  near-field  effects  in 
homogeneous  media,  but  in  heterogeneous  media  it  can  only  be  accounted  for  with  full  3- 
dimensional  modeling  or  full  waveform  inversion  of  the  seismic  wavefield  [Lokmer  and 
Bean,  2010;  Li  et  al.,  201 1] .  A  detailed  3-dimensional  seismic  velocity  model  for  the 
TAG  active  mound  will  therefore  be  required  to  accurately  model  the  effects  of  near- field 
propagation  on  the  seismograms. 


6.1  Potential  source  mechanisms 

Here  we  assess  the  compatibility  of  the  observed  microearthquakes  with  potential  source 
mechanisms  to  identify  plausible  seismogenic  scenarios.  We  do  not  have  enough  data  to 
formally  constrain  the  focal  mechanisms,  but  the  lack  of  5-arrivals  can  be  used  to  rule  out 
any  source  mechanism  involving  shear  failure  across  a  fault  surface.  Thus,  although  the 
hanging  wall  hosting  the  active  mound  is  clearly  tectonized  and  extended  [e.g.,  Kleinrock 
and  Humphris,  1996],  tectonic  faulting  is  inconsistent  with  the  data.  In  addition,  we  do 
not  observe  any  space-time  correlation  between  the  events  observed  at  the  active 
hydrothermal  mound  with  the  larger-scale  seismicity  associated  with  tectonic  extension 
on  the  active  detachment  fault  [deMartin  et  al.,  2007]. 

Thermal  contraction  is  a  process  than  can  generate  very  small  microearthquakes,  but  it 
does  not  appear  to  provide  a  plausible  mechanism  for  our  data.  Thermal  contraction  in  a 
deep-sea  hydrothermal  system  occurs  in  places  where  hot  rock  is  being  cooled  by  cold 
seawater,  and  this  may  generate  small  microearthquakes  [e.g.,  Sohn  et  al.,  2004] . 
However,  the  water-rock  reaction  zone  where  hot  rock  is  cooled  by  hydrothermal  fluids 
at  TAG  is  somewhere  at  or  below  the  detachment  footwall  interface  (i.e.,  depths  >1.5  km 
below  the  seafloor),  not  in  the  upper  200  m  of  crust.  The  only  plausible  way  to  generate 
sustained  rates  of  thermal  contraction  events  in  the  shallow  crust  is  by  cooling  a 
magmatic  intrusion.  For  example,  very  small  (ML  ~  -1)  microearthquakes  with  purely  P- 
wave  arrivals  were  observed  by  seismometers  on  the  Kilauea  Iki  lava  lake  as  it  cooled  in 
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1976  [Chouet,  1979],  and  these  events  bear  some  resemblance  to  our  data.  However,  the 
presence  of  a  magma  body  in  the  shallow  crust  is  precluded  by  the  results  of  an  active 
source  seismic  refraction  study  conducted  during  the  STAG  experiment  [Canales  et  al., 
2007] .  We  note  that  heat  transfer  in  the  shallow  crust  is  a  feature  of  hydrothermal 
circulation  models  for  TAG,  both  through  conductive  heating  from  the  ‘pipe’  that  focuses 
hot  fluids  and  convective  mixing  between  hot  hydrothermal  fluids  and  cold  pore  fluids 
[Teagle  et  al.,  1998;  Lowell  et  al.,  2003].  This  heat  transfer,  which  drives  a  secondary 
circulation  system  in  the  shallow  crust,  does  not  cool  the  crust  and  therefore  does  not 
trigger  thermal  contraction  events. 

Hydrothermal  processes  can  generate  non-double-couple  events  [e.g.,  Ross  et  al.,  1996; 
Miller  et  al.,  1998;  Foulger  et  al.,  2004],  and  the  dense  clustering  of  activity  immediately 
beneath  the  mound  is  consistent  with  a  hydrothermal  mechanism.  Three  potential  source 
mechanisms  associated  with  hydrothermal  processes  have  been  identified:  hydraulic 
fracturing,  bubbles/vibrating  fluids,  and  reaction-driven  cracking.  We  consider  each  of 
these  in  turn.  Hydraulic  fracturing  occurs  when  the  fluid  pressure  somewhere  in  the 
circulation  system  exceeds  the  strength  of  the  local  host  rock.  Fluid  pressures  in  the  TAG 
circulation  system  are  not  well-constrained  by  observations,  but  the  maximum  fluid 
pressures  in  a  high-temperature  deep-sea  hydrothermal  system  are  constrained  to  be  less 
than  or  equal  to  the  buoyancy  pressure,  which  is  the  difference  between  the  pressure  at 
the  base  of  the  recharge  and  discharge  limbs  of  the  convection  cell  [Cann  and  Strens, 
1989].  The  buoyancy  pressure  ( AP )  created  by  a  convection  cell  with  a  depth  of  z,  down¬ 
going  seawater  at  2°C  and  density  psw ,  and  a  discharge  limb  with  hot  fluid  with 
density  pht  at  depth  z  is: 

AP  =  gz(psw  -  pht )  (1) 

where  g  is  gravitational  acceleration.  If  we  consider  a  secondary  circulation  system  up  to 
200  m  deep  (i.e.  the  base  of  the  microearthquake  seismogenic  zone)  with  2°C  seawater 
(Psw  —  1028  kg  m3  )  recharging  the  downgoing  limb  and  200°C  to  400°C  fluids  with 
densities  (pw)  of  874  and  579  kg  m 3,  respectively,  at  the  base  of  the  discharge  limb 
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[Bischoff  and Rosenbciuer,  1985;  Seyfried,  1987],  buoyancy  pressure  generated  is  0.30  to 
0.88  MPa.  Assuming  that  fracturing  occurs  when  buoyancy  pressures  reach  lithostatic, 
which  is  ~43  MPa  with  overburden  at  200  mbsf,  it  is  not  feasible  to  hydrofracture  the 
host  rock  with  such  small  buoyancy  pressures.  We  can  rule  out  hydrofracturing  as  a 
seismic  source  mechanism  in  the  secondary  recharge  system.  In  the  primary  circulation 
system  where  fluids  may  penetrate  along  the  detachment  fault  with  temperatures  of  700- 
900°C  and  fluid  densities  of  300-850  kg  m 3  [Palliser  and  McKibbin,  1998;  Driesner, 
2007],  buoyancy  pressures  generated  are  12-50  MPa.  Again,  buoyancy  pressure  is  much 
less  than  lithostatic  pressure  (>200  MPa)  at  these  depths,  and  hydrofracturing  is  not  a 
likely  source  mechanism  in  the  primary  circulation  system. 

Collapsing  bubbles  and  vibrating  fluids  are  potential  sources  of  hydrothermal 
microearthquakes  and/or  noise  [ Kieffer ,  1984;  Kedar  et  al.,  1996 ;  Kedar  et  al,  1998; 
Vandeulmeulebrouck  et  al.,  2000],  but  they  are  typically  observed  in  subaerial  systems 
where  there  are  large  density  changes  associated  with  the  phase  transition  from  liquid  to 
gas  [Kieffer  and  Delany,  1979].  The  TAG  active  mound  is  located  at  a  depth  of  3670  m 
where  ambient  pressures  are  super-critical  (critical  pressure  in  the  ocean  corresponds  to 
depths  of  ~3000  m),  and  the  volume  change  associated  with  the  transition  from  a  fluid  to 
a  super-critical  vapor  is  ~500  times  smaller  at  these  conditions  compared  to  a  subaerial 
geothermal  field  [e.g.,  Driesner,  2007].  In  addition,  the  frequency  content  of  a  bubble 
collapse  event  is  controlled  to  first-order  by  the  size  (radius)  of  the  bubble,  the  fluid 
density,  and  the  ambient  pressure  [ Rayleigh ,  1917].  For  idealized  Rayleigh  collapse  at  the 
ambient  pressure  of  the  TAG  mound  a  bubble  radius  on  the  order  of  10  m  would  be 
required  to  generate  the  20-40  Hz  signals  we  observed,  which  is  clearly  unrealistic. 
Finally,  none  of  the  TAG  exit  fluid  samples  acquired/analyzed  to-date  are  saturated  in 
any  dissolved  gases  (e.g.,  C02,  H2)  [Charlou  et  al.,  1996].  Vibrating  fluid  mechanics  also 
generate  seismic  harmonics  that  are  absent  from  our  microearthquake  data.  These  pieces 
of  evidence  argue  strongly  against  bubble  collapse  and  vibrating  fluid  mechanisms  for  the 
microearthquakes  we  observed. 
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6.2  Reaction-driven  cracking 


Reaction-driven  cracking  of  the  veins  and  breccias  in  response  to  mineral  deposition 
appears  to  be  the  most  plausible  source  mechanism  for  the  microearthquakes  we 
observed.  The  solid  volume  of  the  mineral  deposit  is  continually  increasing  as 
hydrothermal  minerals  are  deposited  in  flow  veins,  and  this  process  (also  known  as  ‘frost 
heave’  and  ‘salt  weathering’  in  the  literature)  can  generate  large  stresses  that  fracture  the 
host  matrix  [Walder  and  Hallet,  1985;  Scherer,  1999;  2004;  Steiger,  2005a, b].  Cores 
recovered  from  the  TAG  active  mound  during  ODP  Leg  158  revealed  breccias  with  an 
abundance  of  anhydrite  veins,  many  of  which  contain  textural  evidence  for  fracturing. 
Anhydrite  is  a  retrograde  soluble  mineral  that  precipitates  from  seawater  when 
temperatures  rise  above  ~150°C  through  conductive  heating  or  convective  mixing  with 
hydrothermal  fluid,  and  dissolves  at  cooler  temperatures  [Edmond  et  al.,  1995;  Chiba  et 
al.,  1998;  Teagle  et  al.,  1998;  Lowell  and  Yao,  2002].  The  anhydrite  veins  fill  complex, 
multistage  fractures  up  to  45  cm  in  width  [Humphris  et  al.,  1996],  with  textures  typical  of 
open-space  filling.  Individual  growth  bands  reflect  sequential  precipitation  within  cavities 
and  repeated  opening  of  fractures.  Some  anhydrite  veins  contain  angular  clasts  of 
competent,  low-porosity  wallrock,  indicative  of  fracturing  (Figure  11).  In  samples  with 
multiple  veins,  the  host  rock  is  commonly  shattered  and  healed  by  late  anhydrite 
[Humphris  et  al.,  1996].  Textures  of  late  anhydrite  veining  show  that  cementation  and 
fracturing  is  a  continual  process  during  hydrothermally  active  periods.  Anhydrite  was 
recovered  on  broken  clast  surfaces  and  in  veins  and  vugs  throughout  the  stockworks 
zone,  indicating  that  seawater  penetrates  to  depths  of  at  least  120  mbsf  [Humphris , 

Herzig  et  al.,  1996], 

The  ODP  cores  thus  demonstrate  that  anhydrite  precipitation  and  fracturing  of  the  basalt 
matrix  occur  concurrently  within  the  depth  horizon  of  our  microearthquakes.  In  the  cores 
it  was  not  possible  to  distinguish  between  a  hydraulic  vs.  reaction-driven  cracking 
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mechanism,  but  as  we  showed  earlier  hydraulic  pressures  are  not  likely  to  reach 
fracturing  thresholds  in  the  secondary  circulation  system  where  anhydrite  is  deposited. 
Reaction-driven  cracking  is  therefore  the  most  probable  fracturing  mechanism 
accompanying  anhydrite  deposition,  and  here  we  assess  the  compatibility  of  this 
mechanism  with  our  microearthquake  data. 

Data  from  drilled  cores  show  that  the  anhydrite  within  the  mound  precipitated  from 
mixtures  of  hydrothermal  fluid  and  seawater  heated  to  >350°C  [ Mills  et  al.,  1998;  Tivey 
et  al.,  1998]  in  particular  the  analysis  of  fluid  inclusions,  show  that  anhydrite  precipitates 
in  the  secondary  circulation  system  when  seawater  is  conductively  heated  to  108-383  °C. 
Anhydrite  precipitates  from  a  fluid  flowing  through  a  permeable  network  (i.e.  with 
deposition  in  fractures),  which  results  in  expansion  of  the  solid  volume  filling  the 
fracture.  Volume  expansion  increases  stresses  on  the  fracture  surface,  and  if  the  stresses 
are  large  enough  a  cracking  event  is  triggered  that  opens  the  fracture  and  generates  elastic 
waves  (i.e.,  seismic  energy)  [e.g.,  Grossi  et  al.,  1997].  Cracking  maintains  permeability 
and  exposes  fresh  mineral  surfaces  in  the  fracture  network,  most  likely  anhydrite, 
creating  a  positive  feedback  with  fluid  flow  and  mineral  precipitation. 

To  model  the  stresses  generated  by  the  precipitation  of  anhydrite  at  the  TAG  active 

mound  we  consider  mode  I  tensile  cracking  of  a  penny-shaped  fracture  (Figure  1 1).  We 

wish  to  estimate  the  fracture  size  required  to  generate  the  ‘average’  microearthquake 

observed  in  our  study  (M,  =  -0.95,  Ma  =  3.7xl014  dyne-cm).  Consider  a  crack  opening 

under  tension  in  plane  so  that  only  the  strain  component  Ae33  is  nonzero  (Figure  12A). 

The  moment  density  tensor  in  an  isotropic  medium  is 

/AAe  0  0  \ 

m  —  j  0  AAe  0  j  (2) 

\0  0  {A  +  2\£)AeJ 

where  A  and  \i  are  the  Lame  parameters  of  the  medium  (Table  1).  The  moment  tensor  is 
found  by  integrating  over  the  source  volume,  V, 
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Thus,  solving  equation  (5)  for  strain,  we  substitute  the  stress-strain  relationship  into  the 
scalar  moment  equation  to  determine  what  source  volumes  are  required  to  generate  a 
microearthquake  typical  of  those  observed  at  the  TAG  active  mound  (i.e.  M0  =  3.7xl014 
dyne-cm)  for  a  given  stress  drop.  The  relationship  between  earthquake  size,  fracture  size, 
and  stress  drop  is  given  by  Eqn  (6).  In  other  words,  assuming  M0  =  3.7xl014  dyne-cm, 
we  vary  stress  drop  and  solve  for  source  volume, 
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Table  1.  Lame  parameters  and  Young’s  moduli  derived  from  seismic  velocity  and  bulk 
density  measurements  of  hand  core  samples  of  several  rock  types  taken  from  TAG  active 
mound  subsurface  [ Ludwig  et  al.,  1998]. 


Rock  Type 

P  (GPa) 

A  (GPa) 

E  (GPa) 

Core  section(s) 
158-957 

Altered  basalt 

42.91- 

30.42 

41.81 

78.45 

M-9R-1, 61-63  cm; 

46.67 

30.42 

44.23 

78.86 

10R- 1,39-41  cm; 

33.49 

35.67 

84.25 

10R- 1,47-49  cm 

33.49 

38.09 

84.80 

Silicified  wallrock  breccia 

72.90 

41.77 

22.78 

98.27 

E-  8R-1,  10-12  cm 

Chloritized  basalt  breccia 

106.80 

44.85 

22.89 

104.86 

E-  15R-1, 30-32  cm 

Nodular  pyrite-silica 

40.14 

45.74 

17.39 

104.08 

P-  9R-1, 4-6  cm 

breccia 

Pyrite-silica  breccia 

56.31 

74.95 

44.16 

177.68 

P-  12R-2,  67-69  cm 

Measured  anhydrite  veins  in  ODP  cores  have  a  median  thickness  of  L= 4  mm  [Humphris 
et  al.,  1996,  Table  5].  If  we  assume  that  veins  obey  a  linear  fracture-length  relationship, 

L  —  C1  (2c),  where  L  is  the  thickness  or  aperture,  c  is  the  half-length,  and  C1  —  l.Ox 
10-3,  a  geologically  representative  aspect  ratio  [Vermilye  and  Scholz,  1995]  then  vein 
half-length  is  2  m.  We  model  the  source  volume  of  the  crack  (Figure  12B)  on  dimensions 
of  anhydrite  veins  (e.g.  a  penny-shaped  crack  4  mm  thick  with  a  4  m  diameter  has  a 
volume  of  F  =  nc2L  —  5xl0-5  m3.  There  is  no  preferred  orientation  observed  in  vein 
dip  and  dips  range  from  0°  to  90°  [Humphris  et  al.,  1996]  so  we  consider  confining  and 
differential  pressures  at  100  m  below  the  seafloor  and  3630  m  below  the  sea  surface. 
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Table  2.  Source  volume  and  stress  drop  required  to  generate  a  -0.95  magnitude 
microearthquake  by  opening  a  cylindrical  crack  in  different  rock  types.  We  see  that 
results  of  calculations  are  similar  for  the  strength  parameters  of  several  different  rock 
types. 


Rock  Type 

Stress 

Source 

Crack 

Stress 

Source 

Crack 

Drop 

Volume  (m3) 

radius  (m) 

Drop 

Volume  (m3) 

radius  (m) 

Altered  basalt 

8.54 

11.08 

0.58 

4.51 

Differential 

8.27 

10.96 

Confining 

0.56 

4.46 

Pressure  of 

9.74 

11.57 

Pressure 

0.66 

4.71 

2.6  MPa 

9.43 

11.45 

of  39 

0.64 

4.66 

Silicified 

wallrock 

breccia 

12.50 

12.58 

MPa 

0.84 

5.12 

Chloritized 
basalt  breccia 

12.73 

12.65 

0.86 

5.15 

Nodular  pyrite- 
silica  breccia 

13.63 

12.94 

0.92 

5.27 

Pyrite-silica 

breccia 

12.22 

12.23 

0.82 

5.08 

Stress  drop  is  perhaps  the  most  difficult  source  parameter  to  constrain  for  any  earthquake. 
Published  values  in  the  literature  range  from  0.1  and  100  MPa  [e.g.  Abercrombie  and 
Leary,  1993;  Prejean  and  Ellsworth,  2001;  Imanishi  et  al.,  2004],  but  large  uncertainties 
are  typical.  To  cover  the  parameter  space  for  the  TAG  active  mound  we  consider  three 
representative  values:  2.6  MPa  =  the  differential  pressure  between  hydrostatic  and 
lithostatic  pressure  at  a  depth  of  125  m,  39  MPa  =  the  confining  pressure  at  a  depth  of 
125  m,  and  100  MPa  =  overpressures  observed  in  salt  weathering  experiments  [Steiger, 
2005a].  Using  these  stress  drop  values  we  can  calculate  the  size/volume  of  the  associated 
fracture  required  to  generate  a  M,  =  -0.95  microearthquake  (Table  1). 


Table  3.  Source  volume  by  stress  drop  required  to  generate  a  -0.95  magnitude 
microearthquake  in  altered  basalt. 


Stress  Drop 

Source  Volume  (m3) 

Crack  radius  (m) 

100  MPa 

0.23 

3.30 

Confining  Pressure  39  MPa 

0.58 

4.51 

Differential  Pressure  2.6  MPa 

8.54 

11.08 
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The  source  volume  required  to  generate  an  earthquake  of  a  given  size  varies  inversely 
with  stress  drop  (Figure  13;  Table  3).  A  stress  drop  equal  to  the  differential  pressure 
requires  unreasonably  large  source  volumes  to  generate  the  microearthquakes  we 
observe,  but  stress  drops  greater  than  or  equal  to  the  confining  pressure  require 
reasonable  source  volumes  of  <  0.58  m\  We  therefore  conclude  that  reaction-driven 
cracking  is  a  plausible  source  mechanism  for  the  microearthquakes  we  observed  if  the 
stress  drop  for  individual  events  is  greater  than  or  equal  to  the  confining  pressure. 


6.3  Implications  for  Anhydrite  Precipitation 

If  our  hypothesis  that  the  microearthquakes  we  observed  are  generated  by  reaction-driven 
fracturing  in  response  to  anhydrite  precipitation  in  the  secondary  circulation  system  is 
correct,  then  we  can  use  our  observations  to  constrain  the  location  and  volume  of 
anhydrite  deposited  during  our  experiment.  The  microearthquakes  cluster  along  two 
linear  trends  that  intersect  near  the  southwest  comer  of  the  active  mound  (Figure  9). 

These  linear  trends  likely  represent  fissures.  Although  the  region  southwest  of  the  mound 
is  not  adequately  covered  by  the  fine-scale  microbathymetric  map  (Figure  2),  fissures 
imaged  in  the  area  by  ARGO  II  cameras  and  mapped  by  Bohnenstiehl  and  Kleinrock 
[2000]  suggest  that  the  trends  may  relate  to  specific  surface  features.  Under  our 
hypothesis  these  two  linear  trends  demarcate  a  zone  of  focused  anhydrite  deposition 
within  the  depth  interval  from  100-200  mbsf.  In  cross-section  we  can  see  that  the  densest 
cluster  of  microearthquakes  is  focused  in  a  narrow  depth  interval  at  ~  125  mbsf,  and  that 
there  are  essentially  no  events  above  this  Tid’ .  These  patterns  imply  that  recharge  for 
secondary  circulation  is  focused  on  two  intersecting  fissures  just  southwest  of  the  mound, 
and  that  temperatures  greater  than  150°C  required  for  anhydrite  precipitation  are  reached 
at  depths  of  less  than  ~125  m. 

We  can  estimate  the  volume  of  anhydrite  precipitated  during  our  experiment  by  summing 
the  volume  change  associated  with  each  microearthquake  that  we  observed.  Because  of 
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the  trade-off  between  stress  drop  and  source  volume  in  Eq.  6  we  constrain  the  volume 
change  by  considering  stress  drops  between  confining  pressure  (39  MPa)  and 
crystallization  pressures  achieved  during  salt  weathering  experiments  (100  MPa).  Under 
these  two  scenarios,  a  Mh  =  -0.95  (M0  =  3.7x10 14  dyne-cm)  microearthquake  generates 
volume  changes  of  31  cm3  and  58  cm3,  respectively.  If  we  assume  a  seismicity  rate  of  243 
microearthquakes  per  day  (as  observed  at  the  beginning  of  the  experiment  before  one  of 
the  instruments  failed)  then  we  obtain  an  annual  volume  change  of  between  27  m3  to  5 1 
m3.  The  total  volume  of  anhydrite  within  the  active  TAG  deposit  has  been  estimated  to  be 
2xl04m3  [Tivey  et  al.,  1998],  such  that  it  would  take  between  391  to  727  years  to 
precipitate  all  the  anhydrite  in  the  deposit  based  on  our  estimated  precipitation  rates. 


We  can  use  these  precipitation  rate  estimates  to  constrain  the  amount  of  heat  and  the 
volume  of  fluid  required  for  anhydrite  deposition.  Assuming  a  precipitation  rate  ( ranh ) 
of  ~0.0143  mol  anhydrite  per  kilogram  of  fluid  [Tivey  et  al.,  1998]  and  a  molar  volume 
(Yanh)  °f  4.6xl0_5m3  mol-1  [Robie  and  Hemingway ,  1995],  first  we  determine  the 
volume  of  entrained  seawater  required  to  match  our  rate  estimates  of  anhydrite  deposition 
(Ranh)  between  27  m3  to  51  m3  per  year.  After  converting  our  estimated  rate  of 
anhydrite  deposition  to  units  of  m3  s'1,  the  rate  of  seawater  entrainment,  rsw,  required  to 
precipitate  the  anhydrite  is 


Uw  R' 


anh\‘  anhv  anh 


v„ 


(7) 


which  comes  to  an  entrainment  rate  of  1.3  to  2.5  kg  s'1  of  seawater,  on  the  order  of  107kg 
of  seawater  per  year.  By  comparison,  the  mass  of  hydrothermal  fluid  feeding  the  Black 
Smoker  Complex  is  on  the  order  of  109  kg  per  year  assuming  that  the  Black  Smoker 
operates  at  about  100  kg  s-1  [ Hannington  et  al.,  1998].  The  estimated  volume  of 
seawater  entrained  is  ~1%  of  the  flux  of  hydrothermal  fluids  through  the  Black  Smoker 
Complex. 
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We  also  use  our  precipitation  rate  estimates  to  constrain  the  heat  flux  required  for 
anhydrite  deposition.  Following  Tivey  et  ah,  [1998],  if  we  entrain  seawater  at  a  rate 
Osw)  °f  1-3  to  2.5  kg  s'1,  then  the  heat  flux  (H)  required  to  heat  seawater  from  2°  to 
160°C  or  even  350°C,  the  precipitation  temperature  of  anhydrite,  is 

H  =  rswcswAT  (8) 

where  csw  is  the  average  heat  capacity  of  seawater  at  400  bars  ( csw  =4100  J  kg  '  K  1 ) 
[Bischoff  and  Rosenbauer ,  1985] .  The  heat  flux  required  to  deposit  our  estimated  volume 
of  anhydrite  is  0.8  to  3.6  MW,  less  than  1%  of  the  ~1  GW  total  heat  flux  of  the  TAG 
active  mound  [ Wickers ,  2005]. 


6.4  Implications  of  reaction-driven  cracking  in  the  TAG  active  mound  subsurface 

The  seismogenic  zone  of  non-double  couple  microearthquakes  is  a  narrow  and  shallow 
area  arcing  along  the  south  and  west  sides  of  the  TAG  active  mound.  The  spatial 
distribution  of  earthquakes  suggests  that  seismicity  illuminates  an  area  of  dilatant 
cracking  driven  by  intense  mineral  precipitation  which  we  infer  to  be  the  result  of 
convective  and  conductive  heating  of  entrained  seawater  (Figure  14B).  Seismicity  is 
offset  from  the  center  of  the  hydrothermal  upflow  zone  under  the  Black  Smoker 
Complex,  presumably  the  highest  thermal  gradient  and  thermal  stresses  in  the  mound. 
The  seismicity  does  not  align  with  dominant  through-going  fault  orientations  (Figure  14), 
suggesting  that  regional  stresses  have  little  influence  on  focusing  secondary  recharge. 

The  occurrence  of  microearthquakes  randomly  in  time  also  suggests  an  interconnected 
hydraulic  system,  which  is  expected  in  the  highly  permeable  shallow  crust. 

Measurements  of  conductive  heat  flow  on  the  west  side  of  the  sulfide  rubble  plateau  that 
surrounds  the  Black  Smoker  Complex  showed  an  area  of  very  low  heat  flow  20-50  m  in 
1994  (Figure  14A)  [Becker  et  al.,  1996]  supporting  a  region  of  local  seawater 
entrainment  west  of  the  black  smokers.  Three-dimensional  modeling  of  on-bottom 
gravity  measurements  show  that  the  sulfide  mound  has  a  sharp  change  in  thickness  near 
the  center  of  the  mound  from  50  m  to  the  north  to  10  m  to  the  south  [Evans,  1996] .  The 
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non-uniform  distribution  of  inferred  secondary  recharge  and  seawater  entrainment  that 
we  observe  may  be  influenced  by  the  asymmetry  of  the  deposit. 
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Figure  1.  A)  Area  map  of  the  Mid-Atlantic  Ridge  at  26°N  and  the  TAG  hydrothermal 
field  [Jones,  2003;  Bird,  2003].  B)  Bathymetry  of  the  TAG  segment  with 
microearthquake  epicenters  (black  dots)  depicting  a  dome-shaped  detachment  fault 
[deMartin  et  al.,  2007] .  Its  extensional  hanging  wall  hosts  relict  and  active  hydrothermal 
mounds  [Tivey  et  al.,  2003].  The  locations  of  relict  and  shimmering  hydrothermal 
mounds  are  marked  by  black  and  white  circles,  respectively,  and  the  TAG  active  mound 
is  at  the  white  star. 
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Figure  2.  Microbathymetry  of  the  TAG  active  mound  [White  et  al.,  1998;  Roman 
and  Singh,  2007] .  The  surface  expression  of  the  ~200  m  diameter  mound  is  indicated 
by  the  outermost  dashed  white  circle.  High-temperature  discharge  (>360°C)  during 
the  time  period  of  the  STAG  experiment  was  focused  at  the  Black  Smoker  Complex 
(BSC),  whereas  low-temperature  (<100°C)  discharge  was  observed  at  discrete 
locations  across  the  upper  terrace  [Sohn,  2007a].  Locations  of  the  5  OBSs  on  the 
periphery  of  the  mound  are  denoted  by  triangles.  STIR- 12  (black  triangle)  failed  77 
days  into  the  experiment.  Local  fault/fissure  orientations  show  great  variation 
[Bohnenstiehl  and  Kleinrock,  2000] . 
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Figure  3.  Cumulative  event  count  in  time  for  events  arriving  on  four  or  more  OBSs  (n  = 
32,078  events).  Event  rate  is  243  events  per  day  prior  to  the  failure  of  STIR- 12  77  days 
into  the  study,  at  which  time  event  rate  decreases  to  128  events  per  day.  Whale  noise  in 
the  area  begins  in  November,  masking  some  microearthquake  events  and  decreasing  the 
observed  event  rate  to  97  events  per  day. 
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Figure  4.  Multitaper  spectral  estimates  (NW=5)  of  (A)  the  microearthquake  catalog  as  a 
point  process  (10  minute  bins),  (B)  pressure  gauge  data,  and  (C)  predicted  solid  earth  tide 
loading  show  a  lack  of  tidal  influence  on  the  microearthquake  catalog.  Major  tidal 
constituents  (e.g.  O,,  K,,  M2,  S2,  MK,,  M4)  are  labeled  in  blue  lines  [ Foreman ,  1977]. 
There  are  no  highly  coherent  periodicities  (with  threshold  C2  >  0.6)  between  the  (D) 
microearthquake  catalog  and  pressure  gauge  data  nor  between  the  (E)  microearthquake 
catalog  and  solid  earth  tide.  Most  coherences  fall  below  the  95%  level  of  significance 
denoted  by  a  red  line  (i.e.  coherences  are  not  significantly  higher  than  zero). 


D)  Coherence:  Earthquake  catalog  vs.  Pressure  gauge  data 
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Figure  5.  Magnitude  histogram  of  6,207  events  with  hypocenter  estimate  root  mean 
square  errors  less  than  20  ms.  Mean  magnitude  is  -0.95  and  mean  log-moment  is  14.57 
(log  dyne  cm). 
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Figure  6.  Waveforms  and  spectra  of  3  microearthquake  events.  Event  waveforms  are 
from  vertical  channels  of  the  5  OBSs  and  have  different  velocity  scales  to  better  convey 
the  diversity  of  waveforms  across  the  network.  Event  power  spectral  density,  which  is 
calculated  using  the  multi-taper  method  with  NW=4,  is  shown  for  each  waveform.  Events 
have  single-phase  arrivals  and  very  different  waveforms  across  the  network  suggestive  of 
non-double  couple  sources  whose  paths  are  affected  by  local  scattering.  See  Figures  8 
and  9  for  estimated  hypocenters  (labeled  EX1 ,  EX2,  and  EX3). 
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NUMBER  OF  EVENTS 


Figure  7.  RMS  histogram  shows  the  17,487  events  arriving  on  four  or  more  OBSs  with 
hypocenter  estimate  root  mean  square  (RMS)  errors  less  than  200  milliseconds.  Epicenter 
and  hypocenter  maps  depict  the  6,207  events  with  hypocenter  estimates  with  RMS  errors 
less  than  20  ms  (dashed  white  line). 
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Figure  8.  A)  Bathymetric  map  of  6,207  epicenters  with  average  l-o  confidence  ellipses 
(+/-  ~400m).  Depth  sections  (no  vertical  exaggeration)  across  (B)  a  west-to-east  plane 
and  (C)  a  south-to-north  plane  show  the  shallow  extent  of  event  hypocenters,  which 
cluster  in  a  cloud  in  the  top  125  m  of  crust  (blue  dashed  lines).  Hypocenter  and  epicenter 
locations  of  example  events  (Figure  6)  are  labeled  EX1,  EX2,  and  EX3. 
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Figure  9.  A)  Bathymetric  map  of  6,207  epicenters  after  collapse  along  their  confidence 
ellipses  [e.g.  Jones  and  Stewart,  1997].  Depth  sections  (no  vertical  exaggeration)  across 
(B)  a  west-to-east  plane  and  (C)  a  south-to-north  plane  show  the  shallow  extent  of  event 
hypocenters,  which  cluster  in  a  cloud  in  the  top  125  m  of  crust  (blue  dashed  lines). 
Hypocenter  and  epicenter  locations  of  example  events  (Figure  6)  are  labeled  EX1 ,  EX2, 
and  EX3. 
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Figure  10.  Density  map  of  6,207  epicenters  prior  to  applying  the  collapsing  method. 
Note  that  the  results  of  the  collapsing  method  (Figure  9)  serve  to  highlight  the  zone  of 
densest  seismicity  on  the  south  and  west  flanks  of  the  mound. 
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Figure  11.  A)  A  photograph  and  (B)  sketch  of  a  massive  to  semi-massive,  brecciated 
pyrite  associated  with  vein-related,  pyrite- anhydrite  breccias  collected  on  ODP  Leg  158, 
Hole975C  (158-957C-11N-2)  (Figure  16  of  Humphris  et  al.,  [1996]).  Angular  clasts  of 
fine-grained,  massive  pyrite  occur  along  the  edges  of  the  anhydrite  veins  and  are 
extensively  fractured.  The  section  is  collected  from  ~33  mbsf  drilled  from  the  upper 
terrace  ~  20  m  southeast  of  the  Black  Smoker  Complex. 
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Figure  12.  A)  Body  force  representation  for  a  mode  I  crack  opening.  Expansion  is 
illustrated  in  the  vertical  direction.  However  no  preferred  orientation  was  clear  from 
orientations  of  veins  in  drill  cores,  which  had  dips  ranging  from  0°  to  90°  [Humphris  et 
al.,  1996] .  B)  Diagram  of  a  crack  opening  under  tension  in  stockwork. 


Figure  13.  A  plot  of  the  stress  drop  and  source  volume  required  to  generate 
microearthquakes  of  the  range  of  moments  observed  (i.e.  1013-1016  dyne-cm)  in  altered 
basalt.  The  dashed  red  lines  mark  confining  pressure  at  39  MPa  and  differential  pressure 
at  2.6  MPa. 
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Figure  14.  A  depth  section  of  2,415  events  located  in  a  200  m  wide  plane  cutting  north  to 
south  through  the  center  of  the  TAG  active  mound  shows  a  seismogenic  zone  focused  to 
the  south  of  the  deposit  and  <  150  m  below  the  seafloor.  Many  hypocenters  lie  within  the 
anhydrite-rich  zone  and  anhydrite-veined  alteration  zones  of  the  TAG  active  mound 
[Humphris  et  al.,  1995;  Alt  and  Teagle,  1998].  A  well-defined  horizontal  band  of 
seismicity  that  likely  marks  a  zone  of  reaction-driven  cracking  where  secondary 
circulation  recharge  fuels  the  most  intense  anhydrite  deposition.  Comment  about  heat 
source 
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Chapter  4:  Statistical  assessment  of  the  correlation  between 
microearthquake  activity  and  exit-fluid  temperatures  at  the  Trans- 
Atlantic  Geotraverse  active  mound,  Mid-Atlantic  Ridge  (26°08’N) 

Abstract 

We  use  statistical  methods  to  assess  the  impact  of  microearthquake  activity  on  exit-fluid 
temperatures  at  the  Trans-Atlantic  Geotraverse  active  mound.  Two  microearthquake 
catalogs  are  derived  from  a  9-month  deployment  of  ocean  bottom  seismometers  at  the 
TAG  segment  of  the  Mid- Atlantic  Ridge  during  which  exit-fluid  temperature  records 
were  measured  at  discrete  discharge  sites  on  the  upper  terrace  of  the  active  mound.  The 
‘segment-scale’  catalog  contains  microearthquakes  associated  with  detachment  and 
antithetic  normal  faulting  events.  The  ‘mound-scale’  catalog  contains  microearthquakes 
local  to  the  TAG  active  mound  and  intrinsic  to  hydrothermal  processes.  We  test  for 
coherent  frequencies  between  the  catalogs  and  each  temperature  record  and  find  coherent 
frequencies  (C2>0.5)  between  the  mound-scale  microearthquake  catalog  and  seven  of  ten 
temperature  probes.  The  microearthquake  catalogs  are  not  perfectly  Poissonian  point 
processes  (i.e.  coefficents  of  variation  >  1),  so  we  use  a  Rotation  Test,  which  preserves 
the  internal  temporal  structure  of  the  microearthquake  catalogs,  to  test  whether  mean  exit 
fluid  venting  temperatures  are  significantly  higher  following  microearthquake  events. 

We  find  that  2  out  of  20  temperature  probe-channels  show  significantly  higher  mean 
temperatures  (95%  level  of  significance)  following  mound-scale  microearthquakes,  but 
results  are  different  on  two  channels  of  the  same  probes  and  therefore  cannot 
conclusively  show  feedbacks  between  subsurface  seismicity  and  venting  at  those  sites. 
Results  of  initial  tests  suggest  that  more  refined  methods  may  be  necessary  to  define 
spatial  and  temporal  relationships  between  subsurface  cracking  and  hydraulic  pathways 
to  surface  vents. 
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1.  Introduction 


Studies  at  deep-sea  hydrothermal  vents  and  geothermal  sites  on  land  have  shown  that 
seismic  activity  can  affect  hydrothermal  flow  by  perturbing  the  thermal  state  and/or 
permeability  structure  of  the  local  crust  [Aki  et  al.,  1982;  Bame  and  Fehler,  1986; 
Ferrazzini  et  al.,  1990;  Sohn  et  al.,  1998;  Johnson  et  al.,  2000;  Bianco  et  al.,  2004;  Davis 
et  al.,  2004;  Foulger  et  al.,  2004;  Wilcock,  2004;  Cuenot  et  al.,  2006] .  While  these 
observations  provide  intriguing  evidence  for  seismic  perturbation  of  hydrothermal  flow, 
the  extent  to  which  they  may  be  generalized  is  not  clear.  We  do  not  yet  understand  how 
earthquake  characteristics  (e.g.,  size,  location,  source  mechanism)  contribute  to 
hydrothermal  perturbations,  nor  do  we  understand  the  role  of  local  hydrogeology  in  the 
perturbation  process(es).  One  way  to  assess  these  issues  is  with  rigorous  statistical 
analyses  of  hydrothermal  flow  records  and  contemporaneous  earthquake  catalogs  to 
examine  relationships  in  a  universal  framework. 

Long-term  flow  records  for  deep-sea  hydrothermal  vents  are  relatively  scarce  owing  to 
the  logistical  difficulties  and  expense  associated  with  making  these  measurements.  In  the 
past  decade,  however,  several  observatory-style  experiments  providing  concurrent 
measurements  of  local  seismicity  and  hydrothermal  exit-fluid  temperatures  have  been 
conducted  [Sohn  et  al.,  1998;  deMartin  et  al.,  2007;  Sohn,  2007a, b;  Van  Ark  et  al.,  2007], 
and  these  data  provide  an  excellent  opportunity  to  assess  the  statistical  correlation 
between  flow  records  and  earthquake  catalogs.  Here  we  focus  on  data  acquired  from  the 
TAG  active  mound  on  the  Mid-Atlantic  Ridge  (MAR,  26°N)  during  the  Seismicity  and 
Fluid  Flow  of  the  TAG  Hydrothermal  Mound  (STAG)  experiment  from  2003-2004.  The 
experiment  consisted  of  concurrent  deployments  of  thirteen  ocean  bottom  seismometers 
(OBSs),  a  seafloor  pressure  gauge,  eight  high-temperature  probes,  and  twelve  low- 
temperature  probes.  Data  acquired  included  concurrent  catalogs  of  microearthquake  data 
on  both  the  scale  of  the  segment  and  the  scale  of  the  hydrothermal  mound,  a  year-long 
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record  of  seafloor  pressure,  and  year-long  records  of  exit  fluid  temperatures  at  low- 
temperature  and  high-temperature  venting  orifices  on  the  surface  of  the  TAG  mound. 


2.  Hydrothermal  Flow  Records  from  the  TAG  Active  Mound 

The  TAG  active  mound  is  a  ~200  m  diameter,  ~50  m  tall  massive  sulfide  deposit  that 
discharges  hydrothermal  fluids  in  excess  of  360°C  from  a  Black  Smoker  Complex 
(B.S.C.)  of  chimneys  at  the  top  of  the  mound,  and  diffuse,  lower  temperature  (<  100°C) 
fluids  from  numerous  discrete  orifices/cracks  on  the  mound’s  upper  terrace  [Campbell  et 
al.,  1988;  Humphris  et  al.,  1995;  Tivey  et  al.,  1995;  Becker  et  al.,  1996].  On  the  basis  of 
geochronological  studies,  high  temperature  venting  at  the  TAG  active  mound  is  thought 
to  have  begun  roughly  20,000  -  50,000  years  ago  [Lalou  et  al.,  1990,  1993,  1998].  The 
geochronological  data  indicate  that  the  mound  has  experienced  discrete  pulses  of  activity 
every  5,000-6,000  years,  with  quiescent  periods  in  between.  Cooling  of  the  mound  during 
quiescence  leads  to  dissolution  of  anhydrite,  mass  wasting,  and  extensive  reworking  of 
the  sulfides  [Humphris  et  al.,  1995] .  The  most  recent  pulse  of  activity  is  believed  to  have 
begun  ~90  years  ago  [Lalou  et  al.,  1998;  You  and  Bickle,  1998]. 

Exit  fluid  temperatures  were  measured  every  8  min  (on  low-temperature  probes)  or  10 
min  (on  high-temperature  probes)  from  17  sites  distributed  across  the  upper  terrace  of  the 
mound  from  June  2003  to  June  2004  [Sohn,  2007a].  High-temperature  records  were 
obtained  using  Deep  Sea  Power  and  Light  SeaLogger  probes  deployed  in  fractures 
discharging  ~360°C  black  smoker  fluids,  and  low-temperature  records  were  obtained 
using  VEMCO  Ltd.  Minilog  probes  deployed  in  cracks  discharging  ~20°C  diffuse  flow 
fluids.  All  (7)  of  the  high-temperature  records  are  clipped  (low,  <  150°C)  for  extended 
periods  of  time  during  the  experiment,  and  are  therefore  unsuitable  for  correlative 
analyses  with  earthquake  catalogs.  The  ten  low-temperature  records,  however,  provide 
continuous  data  for  the  12-month  observation  period  (Figure  1)  with  only  brief  periods  of 
clipping  (high,  >  60°C)  on  a  few  of  the  probes  (LI,  L2,  L4,  and  L7)  (Figure  2).  Each  low- 
temperature  probe  was  equipped  with  two  sensors  placed  a  few  cm  apart  yielding  two 
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separate  records  for  each  station,  thereby  providing  redundancy  in  case  of  sensor  failure 
and  the  capability  to  measure  fine-scale  spatial  discharge  variations.  In  our  analyses,  we 
treat  each  probe  channel  as  a  separate  record,  such  that  we  have  2x10=20  individual 
temperature  records  to  test  against  the  earthquake  catalogs. 

The  diffuse  flow  records  exhibit  a  complex  mix  of  episodic  and  harmonic  variability  with 
little  apparent  temporal  correlation  between  records  from  different  sites  [ Sohn ,  2007a, b]. 
The  temperature  records  form  highly  skewed,  long-tailed,  and  multi-peaked  empirical 
distributions  that  have  been  modeled  as  mixtures  of  Gamma  densities  [Sohn,  2007b]. 
Individual  temperature  records  exhibit  strong  serial  correlations,  with  the  records 
commonly  alternating  between  background  bottom  water  temperatures  indicative  of  no¬ 
flow  or  possibly  recharge  and  higher  temperatures  indicative  of  discharge.  These 
characteristics  led  Sohn  [2007b]  to  model  the  records  as  Hidden  Markov  Processes  with 
three  to  five  flow  states,  and  to  hypothesize  that  the  processes  triggering  flow 
perturbations  occur  at  very  shallow  depths  within  the  mound. 


3.  Microearthquake  Catalogs 

In  June  2003,  a  network  of  13  ocean  bottom  seismometers  (OBSs)  was  deployed  on  the 
TAG  segment  of  the  MAR  (Figure  3)  as  part  of  the  STAG  Experiment.  The  OBSs 
recorded  continuous  data  at  100  Hz  on  4  channels  (3  component,  4.5  Hz  natural 
frequency  geophone  +  hydrophone)  from  June,  2003  to  February,  2004,  corresponding  to 
the  first  ~9  months  of  the  exit-fluid  temperature  records.  The  network  consisted  of  three 
rings  of  instruments  centered  on  the  active  TAG  mound:  a  7  x  7  km  array  of  four  OBSs 
and  a  2  x  2  km  array  of  four  OBSs  (‘segment-scale’  network),  and  a  200  m  diameter 
array  of  five  OBSs  deployed  on  the  periphery  of  the  TAG  active  mound  (‘mound-scale’ 
subnetwork)  (Figure  2). 

Two  microearthquake  catalogs  have  been  derived  from  the  OBS  data.  deMartin  et  al. 
[2007]  generated  a  catalog  of  19,232  regional  microearthquakes  that  were  recorded  by  all 
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of  the  OBSs  in  the  network  (Figure  3B).  These  events,  with  local  magnitudes  ranging 
from  1  <  Ml  <4,  represent  tectonic  extension  on  the  active  detachment  fault  and  the 
deforming  fault  footwall.  More  recently,  Pontbriand  and  Sohn  [submitted]  generated  a 
catalog  of  32,078  microearthquakes  that  were  only  recorded  by  the  mound-scale  sub¬ 
network.  These  events,  with  local  magnitudes  of  -1.4  <  ML  <0.5,  cluster  along  two  linear 
trends  that  intersect  near  the  southwest  corner  of  the  active  mound  (Figure  3C),  and  are 
believed  to  represent  cracking  generated  by  anhydrite  deposition  in  the  secondary 
circulation  system.  The  two  catalogs  are  complementary  in  that  the  regional  events  are 
associated  with  segment-scale  tectonic  extension  whereas  the  local  events  are  associated 
with  mound-scale  fluid  flow  and  mineral  deposition. 

Both  microearthquake  catalogs  exhibit  relatively  steady  seismicity  rates  with  no  apparent 
mainshock-aftershock  behavior.  The  segment-scale  catalog  has  a  rate  of  ~88  events/day, 
corresponding  to  a  mean  inter-event  time  of  17.5  min,  whereas  the  mound-scale  catalog 
has  a  rate  of  ~243  events/day  (until  the  detection  threshold  increases  as  a  result  of 
instrument  failure),  corresponding  to  a  mean  inter-event  time  of  10.3  min  (Figure  4).  The 
tendency  for  earthquakes  to  cluster  in  time  is  described  by  the  coefficient  of  variation 
(Cv)  [Kagan  and  Jackson,  1991],  which  is  the  standard  deviation  of  inter-event  times 
divided  by  the  mean  inter-event  time.  If  the  temporal  behavior  is  completely  random  (i.e., 
Poisson  process)  the  coefficient  of  variation  takes  a  value  of  one.  If  the  events  cluster  in 
time  the  coefficient  of  variation  will  be  greater  than  one,  while  by  contrast  if  the  events 
occur  periodically,  the  coefficient  of  variation  approaches  zero.  Coefficients  of  variation 
for  the  segment-scale  and  mound-scale  catalogs  are  2.42  and  1.25,  respectively, 
indicating  that  both  catalogs  exhibit  statistically  significant  clustering  in  time.  The  inter¬ 
event  times  of  mound-scale  events  are  much  shorter,  as  we  would  expect  from 
magnitude -rate  models  [e.g.  Ogata,  1988],  with  maximum  inter-event  times  an  order  of 
magnitude  shorter  than  those  of  the  segment-scale  catalog  (~4  hours  vs  ~40  hours).  Inter¬ 
event  times  from  neither  catalog  fit  an  exponential  distribution,  providing  additional 
evidence  that  they  cannot  be  modeled  as  a  Poisson  process. 
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Loading  from  ocean  and  earth  tides  has  been  hypothesized  to  trigger  microearthquake 
activity  at  mid-ocean  ridge  hydrothermal  fields  [Tolstoy  et  al.,  2002;  Crone  and  Wilcock, 
2005].  Both  the  segment-scale  and  mound-scale  catalogs  in  this  study  have  coefficients 
of  variation  greater  than  one,  arguing  against  perfectly  periodic  behavior  as  would  be 
expected  from  tidal  triggering,  but  we  use  spectral  methods  to  explicitly  investigate 
periodic  behavior  at  tidal  frequencies.  We  estimate  the  power  spectral  density  of  the 
catalogs  using  the  method  of  Brillinger  [1974]  wherein  the  point  process  is  converted  to  a 
continuous  process  by  counting  the  number  of  events  in  discrete  time  windows.  Each 
catalog  is  binned  on  the  sampling  time  scale  of  the  pressure  gauge  data  (i.e.  10  min  bins). 
Power  spectra  for  the  mound-scale  microearthquake  catalog  (Figure  5A)  does  not  exhibit 
peaks  at  any  of  the  tidal  frequencies,  indicating  that  tidal  loading  does  not  play  a 
significant  role  in  triggering  microearthquakes  associated  with  mound-scale  fluid  flow. 
This  result  is  reinforced  by  the  lack  of  coherency  between  the  catalog  with  concurrently 
measured  ocean  tide  data  (bottom  pressure)  and  predicted  solid  earth  tides  (Figure  5E, 
Figure  5F).  Power  spectra  for  the  segment-scale  microearthquake  catalog  (Figure  5B) 
exhibits  a  peaks  at  the  M2  tidal  constituent,  indicating  that  tidal  loading  at  the  principal 
lunar  semi-diurnal  tidal  periodicity  may  play  a  role  in  triggering  microearthquakes 
associated  with  segment-scale  deformation.  This  result  is  reinforced  by  the  coherency 
between  the  catalog  with  the  predicted  solid  earth  tide  (Figure  5H),  but  is  not  present  in 
the  comparison  between  the  segment-scale  catalog  and  pressure  gauge  data  (Figure  5E). 
Coherence  at  M2  is  robust  to  higher  levels  of  spectral  smoothing. 


4.  Statistical  Tests  and  Initial  Results 

We  use  spectral  statistical  methods  and  a  non-parametric  rotation  test  to  compare 
segment-scale  and  mound-scale  microearthquakes  with  exit-fluid  temperature  records  to 
assess  potential  feedbacks  between  microearthquakes  and  fluid  flow.  Here  we  focus  on 
temporal  correlations,  though  spatiotemporal  models  may  be  possible  in  the  future.  Data 
from  each  temperature  probe  channel  are  treated  individually  against  the  microearthquake 
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catalogs,  such  that  each  observation  site  is  tested  twice  (once  for  each  channel).  In  order 
to  have  ‘significant’  correlation  at  a  given  site  we  require  significant  correlations  for  both 
channels  on  an  exit  fluid  temperature  probe. 


4.1  Spectral  coherence:  Microearthquakes  vs.  exit  fluid  temperatures 
We  use  a  multi-taper  spectral  method  to  compare  the  spectral  content  of  the 
microearthquake  catalogs  and  exit  fluid  temperature  records.  Results  shown  are  robust 
for  both  NW=3  and  NW=5  levels  of  spectral  smoothing,  and  coherences  reported  are 
consistently  high  over  frequency  bands  as  opposed  to  at  single  frequencies.  Earthquake 
data  are  binned  at  the  sample  interval  of  the  exit  fluid  records  (8  minutes).  Periods  with 
coherences  greater  than  0.5  (and  significantly  greater  than  zero  ,  C 2  >  0,  with  95% 
confidence  [Percival  and  Walden,  1993])  are  found  in  Table  1  for  the  mound-scale 
catalog. 
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Table  1.  A  comparison  of  the  mound-scale  microearthquake  catalog  to  each  exit  fluid 
temperature  probe  shows  that  some  coherent  periods  (C 2  >  0.5)  exist  on  both  probe 
channels.  Coherences  reported  are  a  result  of  multitaper  spectral  analysis  with  NW=3, 
but  results  presented  are  robust  at  both  NW=3  and  NW=5  levels  of  smoothing. 


Probe  (vs.  mound-scale 

Period 

Coherence  with 

Coherence  with 

catalog) 

Channel  1  record 

Channel  2  record 

LI 

19.1  min 

0.67 

0.65 

L4 

3.1  hours 

0.71 

0.78 

1.1  hours 

0.74 

0.66 

20.8  min 

0.69 

0.73 

L5 

1 .2  days 

0.50 

0.92 

L8 

23.7  min 

0.61 

0.55 

L9 

55.2  min 

0.62 

0.75 

Lll 

1 .8  hours 

0.76 

0.77 

38.9  min 

0.57 

0.56 

21.3  min 

0.51 

0.51 

L12 

1 .5  hours 

0.81 

0.72 

23.5  min 

0.81 

0.67 

17.54  min 

0.62 

0.58 

Exit  fluid  probes  that  show  consistent  coherent  periods  with  the  mound-scale 
microearthquake  catalog  on  both  channels  are  LI  (Figure  6),  L4  (Figure  7),  L5  (Figure 
8),  L8  (Figure  9),  L9  (Figure  10),  Lll  (Figure  11),  and  L12  (Figure  12).  Interpretation  of 
periodicities  is  difficult,  though  coherent  periods  in  the  range  of  19-24  min  are  seen  on 
LI ,  L4,  L8,  LI  1 ,  and  L12.  Only  L5  shows  a  coherent  period  on  the  length  scale  of  a  day, 
at  1.2  days,  which  is  slightly  longer  than  the  period  of  O,,  the  lunar  diurnal  tidal 
constituent. 


A  comparison  of  the  segment-scale  microearthquake  catalog  to  each  exit  fluid 
temperature  probe  finds  that  no  coherent  periods  (C 2  >  0.5)  exist  on  both  probe 
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channels.  However  one  probe,  LI  1 ,  shows  a  small  level  of  coherence  (C 2  =  0.4)  with 
the  segment-scale  catalog  at  a  period  12.4  hours,  the  lunar  semidiurnal  tidal  period.  This 
is  in  agreement  with  the  observation  of  coherency  between  the  segment-scale  catalog  and 
the  predicted  solid  earth  tide  at  M2  The  exit  fluid  temperature  record  on  probe  LI  1  has 
been  shown  to  have  a  strong  spectral  peak  at  M2  [  Sohn ,  2007a] . 


4.2  N  onparametric  rotation  test  for  behavioral  point  process  data 
Suppose  that  microearthquake  events  are  observed  over  a  period  ( 0,7)  where  in  this  case 
T  is  9  months  from  the  start  of  the  experiment.  We  have  discussed  how  the 
microearthquake  catalogs  contain  clustering,  so  assumptions  of  Poissonian  point 
processes  are  not  valid  when  correlating  the  catalogs  to  exit  fluid  temperature  records.  In 
order  to  make  this  comparison,  we  use  a  rotation  test  that  preserves  the  interevent  times 
of  the  microearthquake  events,  but  “randomizes”  the  start  time  of  the  entire  catalog  with 
respect  to  the  exit  fluid  temperature  catalog.  In  this  test  microearthquake  events  are 
treated  as  having  no  duration.  This  approximation  is  reasonable  because  the  duration  of 
events  (~1  sec)  is  small  in  relation  to  the  interval  over  which  we  are  looking  for  effects  in 
temperature  (8  -64  minutes).  A  point  of  consideration  is  whether  randomly  selected  time 
windows  will  contain  an  earthquake  and  thus  be  un-representative  of  “non-event”  exit 
fluid  temperatures.  If  this  consideration  is  factual,  we  would  expect  to  always  accept  the 
null  hypothesis  that  there  is  no  significant  difference  in  exit  fluid  temperatures  following 
microearthquakes  and  a  randomly  generated  distribution  of  mean  exit  fluid  temperatures. 
However,  we  find  that  this  is  not  the  case  so  time  windows  up  to  64  minutes  appear 
appropriate  for  gauging  the  responsiveness  of  venting  temperatures  to  microearthquakes. 
The  procedure  of  this  test  is  to  randomly  “rotate”  times  of  microearthquakes  to  form  a 
distribution  of  mean  temperatures,  thereby  preserving  the  internal  structure  of  the 
microearthquake  catalog  [Deruiter  and  Solow,  2008] . 

First  we  determine  the  optimal  time  window  following  an  earthquake  event  over  which 
exit  fluid  temperature  is  maximized.  We  do  so  by  finding  the  length  of  the  time  window 
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that  maximizes  the  average  exit  fluid  temperature  following  all  the  microearthquake 
events  in  the  catalog.  Then  we  randomly  rotate  the  times  of  the  microearthquake  events, 
determine  the  optimal  time  delay,  and  calculate  the  mean  temperature  at  these  times  [c.f. 
Deruiter  and  Solow,  2008] .  We  repeat  this  procedure  1000  times  to  approximate  the 
distribution  of  mean  exit  fluid  temperature.  Finally,  we  test  whether  the  mean  exit  fluid 
temperature  following  earthquake  events  is  significantly  higher  than  that  of  the 
bootstrapped  distribution  using  a  non-parametric  rank  test.  This  test  is  repeated  for  each 
exit  fluid  temperature  probe  record. 

Test  hypotheses: 

H0:  Mean  exit  fluid  temperature  following  earthquake  events  is  not  statistically  different 
from  background  temperature  fluctuations. 

H,:  Mean  exit  fluid  temperature  is  significantly  higher  following  microearthquake  events. 

Test  statistic: 

Nonparametric  rank  test  statistic  at  95%  confidence  level  (one-sided  test). 

If  mean  exit  fluid  temperature  following  earthquakes  is  greater  than  the  test  statistic,  we 
reject  H0. 
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Table  3.  Results  of  rotation  test  between  the  mound-scale  microearthquake  catalog 
and  exit  fluid  temperature  records.  If  we  reject  the  null  hypothesis,  then  mean  exit 
fluid  temperature  is  significantly  higher  following  mound-scale  microearthquake  events. 
If  we  do  not  reject  the  null  hypothesis,  then  mean  exit  fluid  temperatures  are  not 
significantly  higher  following  mound-scale  microearthquakes. 


Temp. 

Probe 

Optimal 

time 

window 

(min) 

Mean  temp, 
following 
earthquakes 
(°C) 

Confidence 

level 

Bootstrap 

samples 

Critical 
value  (°C) 

Conclusion 

LI -Chant 

48 

3.718 

95% 

607 

3.720 

Do  not  reject  H0 

Ll-Chan2 

8 

3.970 

“ 

800 

3.974 

Do  not  reject  H0 

L2-Chanl 

64 

8.864 

“ 

129 

8.929 

Do  not  reject  H0 

L2-Chan2 

56 

14.613 

“ 

1000 

14.550 

Reject  H0 

L4-Chanl 

8 

14.488 

“ 

413 

14.522 

Do  not  reject  H0 

L4-Chan2 

8 

7.563 

“ 

1000 

7.512 

Reject  H0 

L5 -Chant 

64 

5.079 

“ 

202 

5.089 

Do  not  reject  H0 

L5-Chan2 

8 

5.929 

“ 

67 

5.935 

Do  not  reject  H0 

L6-Chanl 

32 

4.974 

“ 

134 

4.975 

Do  not  reject  H0 

L6-Chan2 

64 

6.556 

“ 

73 

6.558 

Do  not  reject  H0 

L7 -Chant 

64 

34.799 

“ 

99 

34.821 

Do  not  reject  H0 

L7-Chan2 

8 

40.022 

“ 

178 

40.054 

Do  not  reject  H0 

L8-Chanl 

16 

16.367 

“ 

54 

16.371 

Do  not  reject  H0 

L8-Chan2 

24 

12.890 

“ 

70 

13.121 

Do  not  reject  H0 

L9-Chanl 

40 

4.416 

“ 

116 

4.417 

Do  not  reject  H0 

L9-Chan2 

48 

4.394 

“ 

142 

4.450 

Do  not  reject  H0 

Lll- 

Chanl 

64 

6.870 

53 

6.909 

Do  not  reject  H0 

Lll- 

Chan2 

64 

6.812 

215 

3.477 

Do  not  reject  H0 

L12- 

Chanl 

48 

3.460 

441 

3.618 

Do  not  reject  H0 

LI  2- 
Chan2 

48 

3.418 

210 

3.421 

Do  not  reject  H0 

114 


When  accounting  for  inter-event  time  structure  of  the  mound-scale  microearthquake 
catalog,  only  tests  on  L2-Channel  2  and  L4-Channel  2  (Figure  13)  find  significantly 
higher  mean  exit  fluid  temperatures  following  mound-scale  microearthquake  events.  We 
do  not  see  similar  behavior  on  both  channels  of  these  two  probes. 

There  is  also  a  possible  model  in  which  microearthquakes  increase  permeability  and 
amounts  of  entrained  seawater,  thereby  lowering  exit  fluid  temperatures.  The  rotation 
test  method  is  repeated  to  look  for  evidence  of  significantly  lower  temperatures  following 
microearthquakes.  Results  (Table  4)  show  no  evidence  of  lower  temperatures  on  any  of 
the  exit  temperature  probes  following  local  microearthquakes. 
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Table  4.  Results  of  rotation  test  #2  between  the  mound-scale  microearthquake 
catalog  and  exit  fluid  temperature  records.  If  we  reject  the  null  hypothesis,  then  mean 
exit  fluid  temperature  is  significantly  lower  following  mound-scale  microearthquake 
events.  If  we  do  not  reject  the  null  hypothesis,  then  mean  exit  fluid  temperatures  are  not 
significantly  lower  following  mound-scale  microearthquakes. 


Temp.  Probe 

Optimal 

time 

window 

(min) 

Mean  temp. 

following 

earthquakes 

(°C) 

Confidence 

level 

Critical 

value  (°C) 

Conclusion 

LI -Chant 

16 

3.715 

95% 

3.309 

Do  not  reject  H0 

Ll-Chan2 

64 

3.965 

3.536 

Do  not  reject  H0 

L2-Chanl 

16 

8.860 

8.350 

Do  not  reject  H0 

L2-Chan2 

40 

14.606 

9.305 

Do  not  reject  H0 

L4-Chanl 

64 

14.477 

7.608 

Do  not  reject  H0 

L4-Chan2 

64 

7.558 

“ 

4.438 

Do  not  reject  H0 

L5 -Chant 

24 

5.075 

4.830 

Do  not  reject  H0 

L5-Chan2 

40 

5.923 

5.922 

Do  not  reject  H0 

L6-Chanl 

8 

4.971 

4.762 

Do  not  reject  H0 

L6-Chan2 

8 

6.551 

6.551 

Do  not  reject  H0 

L7 -Chant 

8 

34.776 

34.142 

Do  not  reject  H0 

L7-Chan2 

32 

40.007 

37.366 

Do  not  reject  H0 

L8-Chanl 

8 

16.364 

“ 

16.363 

Do  not  reject  H0 

L8-Chan2 

8 

12.885 

12.879 

Do  not  reject  H0 

L9-Chanl 

8 

4.415 

“ 

4.192 

Do  not  reject  H0 

L9-Chan2 

16 

4.391 

4.115 

Do  not  reject  H0 

Lit -Chant 

24 

6.862 

6.851 

Do  not  reject  H0 

Lll-Chan2 

24 

6.804 

6.803 

Do  not  reject  H0 

L12-Chanl 

16 

3.458 

3.064 

Do  not  reject  H0 

L12-Chan2 

8 

3.416 

66 

3.054 

Do  not  reject  H0 
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In  a  similar  comparison  of  the  segment-scale  microearthquake  catalog  to  each  exit  fluid 
temperature  probe,  each  rotation  test  concluded  that  accept  the  null  hypothesis.  In  other 
words,  no  exit  fluid  temperature  probes  appeared  to  record  significantly  higher 
temperatures  following  segment-scale  earthquakes. 


Discussion 

Given  the  results  of  spectral  analysis  and  coherence  measurements  between  the  mound- 
scale  microearthquake  catalog  and  exit  fluid  temperature  records,  we  suggest  that  exit 
fluid  temperatures  may  experience  some  feedbacks  from  microearthquakes,  which  we 
suggest  are  caused  by  the  opening  of  fluid  flow  pathways  in  response  to  reaction-driven 
cracking  of  veins  of  anhydrite  deposition.  Interpretation  of  coherent  periodicities  is 
difficult,  though  coherent  periods  in  the  range  of  19-24  min  are  seen  on  both  channels  of 
five  temperature  probes  (LI ,  L4,  L8,  LI  1 ,  and  L12).  Similarly  short  periods  of  22  min, 
26.5  min,  33.4  min,  37.7  min,  53.2  min,  and  90  min  are  evident  in  pressure  gauge  data 
recorded  near  the  TAG  active  mound.  Periods  in  pressure  gauge  data  are  interpreted  to 
represent  vertical  motion  of  the  seafloor  in  response  to  hydrothermal  flow  -  similar  to 
ground  surface  deformation  at  terrestrial  geysers  [c.f.  Sohn  et  al.,  2009] .  Perhaps  a  related 
process,  such  as  a  pressure  cycle  associated  with  flux  of  fluid  or  heat  into  the  subsurface 
of  the  hydrothermal  mound,  may  influence  the  mound-scale  microearthquakes,  exit  fluid 
temperatures,  and  pressure  gauge  data,  but  much  more  detailed  analysis  is  required  to 
better  understand  correlations  in  the  dataset.  Although  two  probe-channels  at  different 
sites  find  significantly  higher  mean  exit  fluid  temperatures  following  mound-scale 
microearthquakes,  the  rotation  tests  fail  to  find  any  sites  where  both  probe  channels  test 
‘positive’  against  the  catalog,  suggesting  that  the  mound-scale  catalog  as  a  whole  did  not 
have  a  measurable  impact  on  the  discharge  temperatures  that  we  measured.  Considering 
the  complex  hydraulic  structure  and  wide  spatial  distribution  of  mound-scale 
microearthquakes  with  respect  to  the  temperature  probes,  finding  a  response  in 
temperature  on  only  two  probes  for  the  rotation  test  is  a  significant  discovery  of  signal  in 
the  noise.  More  refined  comparisons  of  the  spatiotemporal  behavior  of  exit  fluid 
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temperatures  and  shallow  mound-scale  microearthquakes  is  warranted  (e.g.  comparison 
to  state  changes  in  HMMs). 

Tidal  signals  are  apparent  in  the  variability  of  exit  fluid  temperatures  and  in  the  segment- 
scale  microearthquake  catalog  (as  M2).  The  mound-scale  microearthquake  catalogs  lack 
coherence  with  major  tidal  constituents,  however.  Thus  tidal  loading  appears  to 
contribute  to  regional  fault  failure  on  a  lunar  semidiurnal  period,  but  does  not  affect  the 
timing  of  dilatant  fracture  opening.  No  other  notable  correlations  are  found  between  the 
segment-scale  microearthquakes  and  exit  fluid  venting  temperatures  at  the  TAG  active 
mound. 
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Figure  1.  Low-temperature  exit  fluid  temperature  records  displayed  clockwise  around  the  mound  starting 
with  probe  L12  in  the  ODP  6  depression.  Two  records  are  shown  for  each  probe  (left).  Time  series  data  and 
(right)  data  histograms.  The  y-axis  of  the  histograms  represents  normalized  frequency,  which  is  the  number 
of  observations  in  each  data  bin  divided  by  the  total  number  of  observations  (data  are  binned  in  digitization 
bins  of  the  data  loggers).  Many  of  the  records  contained  large  fractions  of  measurements  in  the  first  few 
data  bins  corresponding  to  background  water  column  temperature  levels  (>2.7°C),  and  in  such  cases  a 
small  y  axis  interval  of  0-0. 1  is  retained  to  bring  out  the  details  of  the  histogram  shape,  but  the  fraction  of 
observations  in  the  background  temperature  bins  is  explicitly  listed  (e.g.,  the  two  probes  at  site  L12 
returned  background  temperature  measurements  45%  and  56%  of  the  time,  respectively. 

Time-series  Histograms 


Figure  2.  Bathymetry  of  the  area  surrounding  the  TAG  active  mound  [Roman  and  Singh , 
2007] .  The  active  hydrothermal  deposit  (white  dashed  circle  labeled  MOUND)  is  the  200 
m  diameter  feature  with  two  terraces  surrounded  by  5  OBSs  of  the  mound-scale  network 
(triangles).  STIR- 12  (black  triangle)  failed  77  days  into  the  experiment.  The  main 
chimney,  or  Black  Smoker  Complex  (B.S.C.)  sits  atop  the  upper  terrace  within  the  small 
dashed  white  circle.  Ten  low  temperature  probes  (black  circles)  are  deployed  on  the 
upper  terrace  of  the  sulfide  deposit. 
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Figure  3.  (A)  Area  map  of  the  Mid- Atlantic  Ridge  at  26°N  and  the  TAG  hydrothermal 
field.  (B)  Bathymetry  of  the  TAG  segment  with  microearthquake  epicenters  (black  dots) 
detected  on  the  segment-scale  network  depicting  a  dome-shaped  detachment  fault  and 
antithetic  normal  faults  [deMartin  et  al.,  2007] .  The  extensional  hanging  wall  hosts  relict 
and  active  hydrothermal  mounds  [Tivey  et  al.,  2003] .  The  locations  of  relict  and 
shimmering  hydrothermal  mounds  are  marked  by  black  and  white  circles,  respectively, 
and  the  TAG  active  mound  is  at  the  white  star.  (C)  Bathymetric  map  of  6,207  epicenters 
detected  on  the  mound-scale  network  locate  in  a  seismogenic  zone  focused  to  the 
southwest  of  the  deposit  and  <  150  m  below  the  seafloor.  Many  hypocenters  lie  within 
the  anhydrite-rich  zone  and  anhydrite- veined  alteration  zones  of  the  TAG  active  mound 
[Humphris  et  al.,  1995;  Alt  and  Teagle,  1998]. 
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Figure  4.  Cumulative  event  count  in  time  for  mound-scale  events  (blue,  n  =  32,078 
events)  and  segment-scale  events  (red,  n  =  19,827  events).  Mound-scale  event  rate  is  243 
events  per  day  prior  to  the  failure  of  STIR-12  77  days  into  the  study,  at  which  time  event 
rate  decreases  to  128  events  per  day.  Average  segment-scale  event  rate  is  80  events  per 
day.  Whale  noise  in  the  area  begins  in  November,  masking  some  microearthquake  events 
and  decreasing  the  observed  mound-scale  event  rate  to  97  events  per  day. 
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Figure  5.  Normalized  multitaper  spectral  estimates  (NW=5)  of  the  (A)  mound-scale  catalog,  (B)  segment- 
scale  catalog,  (C)  pressure  gauge  data,  and  (D)  solid  earth  tidal  prediction  at  tidal  periods.  Coherences  of 
the  (E)  mound-scale  catalog  vs  pressure  gauge  data,  (F)  mound-scale  catalog  vs  solid  earth  tides  show  little 
coherent  frequencies  at  tidal  periods.  Major  tidal  constituents  (e.g.  Ol5  K1;  M2,  S2,  MK,,  M4)  are  labeled  in 
blue  lines  [Foreman,  1977].  Coherence  at  tidal  periodicity  is  evident  between  the  segment-scale 
microearthquakes  and  the  solid  earth  tides  at  M2,  the  principal  lunar  semi-diurnal  tide. 
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Figure  6.  The  mound-scale  microearthquake  catalog  as  a  point  process  (8  minute  bins) 
and  exit  fluid  temperature  record  for  both  channels  of  probe  LI  are  shown  below.  The 
coherence  plots  (at  NW=3)  between  each  probe  channel  and  the  mound-scale  catalog 
show  the  19  min  period  at  which  significant  coherence  exists  between  both  probes  and 
the  catalog  (orange  line).  Results  shown  are  robust  to  spectral  smoothing  level.  The  zero 
significance  level  for  coherence  ( C 2  >  0  with  95%  confidence)  is  depicted  as  a  grey  line. 
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Figure  7.  The  mound-scale  microearthquake  catalog  as  a  point  process  (8  minute  bins) 
and  exit  fluid  temperature  record  for  both  channels  of  probe  L4  are  shown  below.  The 
coherence  plots  (at  NW=3)  between  each  probe  channel  and  the  mound-scale  catalog 
show  3.1  hour,  1.1  hour,  and  20.8  min  periods  at  which  significant  coherences  exist 
between  both  probes  and  the  catalog  (orange  lines).  Results  shown  are  robust  to  spectral 
smoothing  level.  The  zero  significance  level  for  coherence  ( C 2  >  0  with  95% 
confidence)  is  depicted  as  a  grey  line. 
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Figure  8.  The  mound-scale  microearthquake  catalog  as  a  point  process  (8  minute  bins) 
and  exit  fluid  temperature  record  for  both  channels  of  probe  L5  are  shown  below.  The 
coherence  plots  (at  NW=3)  between  each  probe  channel  and  the  mound-scale  catalog 
show  the  1 .2  day  period  at  which  significant  coherence  exists  between  both  probes  and 
the  catalog  (orange  line).  Probe  L5  is  the  only  site  where  a  longer  coherent  period  is 
measured,  and  the  1.2  day  period  is  slightly  longer  than  the  lunar  diurnal  tidal  constituent 
(O , ) .  Results  shown  are  robust  to  spectral  smoothing  level.  The  zero  significance  level 
for  coherence  (C2  >  0  with  95%  confidence)  is  depicted  as  a  grey  line. 
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Figure  9.  The  mound-scale  microearthquake  catalog  as  a  point  process  (8  minute  bins) 
and  exit  fluid  temperature  record  for  both  channels  of  probe  L8  are  shown  below.  The 
coherence  plots  (at  NW=3)  between  each  probe  channel  and  the  mound-scale  catalog 
show  the  23.7  min  period  at  which  significant  coherence  exists  between  both  probes  and 
the  catalog  (orange  line).  Results  shown  are  robust  to  spectral  smoothing  level.  The  zero 
significance  level  for  coherence  ( C 2  >  0  with  95%  confidence)  is  depicted  as  a  grey  line. 
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Figure  10.  The  mound-scale  microearthquake  catalog  as  a  point  process  (8  minute  bins) 
and  exit  fluid  temperature  record  for  both  channels  of  probe  L9  are  shown  below.  The 
coherence  plots  (at  NW=3)  between  each  probe  channel  and  the  mound-scale  catalog 
show  the  period  of  55.2  min  at  which  significant  coherence  exists  between  both  probes 
and  the  catalog  (orange  line).  Results  shown  are  robust  to  spectral  smoothing  level.  The 
zero  significance  level  for  coherence  (C2  >  0  with  95%  confidence)  is  depicted  as  a  grey 
line. 
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Figure  11.  The  mound-scale  microearthquake  catalog  as  a  point  process  (8  minute  bins) 
and  exit  fluid  temperature  record  for  both  channels  of  probe  Lll  are  shown  below.  The 
coherence  plots  (at  NW=3)  between  each  probe  channel  and  the  mound-scale  catalog 
show  periods  of  1.8  hours,  38.9  min,  and  21.3  min  at  which  significant  coherences  exist 
between  both  probes  and  the  catalog  (orange  lines).  Results  shown  are  robust  to  spectral 
smoothing  level.  The  zero  significance  level  for  coherence  (C2  >  0  with  95% 
confidence)  is  depicted  as  a  grey  line. 
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Figure  12.  The  mound-scale  microearthquake  catalog  as  a  point  process  (8  minute  bins) 
and  exit  fluid  temperature  record  for  both  channels  of  probe  LI 2  are  shown  below.  The 
coherence  plots  (at  NW=3)  between  each  probe  channel  and  the  mound-scale  catalog 
show  periods  of  1.5  hours,  23.5  min,  and  17.5  min  at  which  significant  coherences  exist 
between  both  probes  and  the  catalog  (orange  lines).  Results  shown  are  robust  to  spectral 
smoothing  level.  The  zero  significance  level  for  coherence  ( C 2  >  0  with  95% 
confidence)  is  depicted  as  a  grey  line. 
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Figure  13.  Rotation  test  distributions  (n=1000)  for  (A)  L2-  Channel  2  and  (B)  L4- 
Channel  2  compared  to  the  mound-scale  microearthquake  catalog.  Both  tests  reject  the 
null  hypothesis  and  conclude  that  mean  exit  fluid  temperature  (°C)  is  significantly  higher 
following  mound-scale  microearthquake  events. 
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